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ABSTRACT 


The  purpose  of  this  program  was  to  predict  the  performance  of  eight  marine  composites  and  one 
textile  wall  covering  in  the  ISO  9705  room-comer  test  using  the  fire  growth  model  developed  by 
Quintiere  [Quintiere,  1993].  Four  sets  of  material  properties  were  developed  using  data  obtained  from  the 
Cone  Calorimeter  and  the  Lateral  Ignition  and  Flame  spread  Test  (LIFT)  apparatus.  These  properties 
were  input  into  a  FORTRAN  version  of  Quintiere’ s  model  previously  published  by  Dillon  [Dillon,  1998], 
The  predicted  heat  release  rate  was  compared  with  the  actual  heat  release  rate  measured  in  ISO  9705  tests. 
These  initial  simulations  did  not  provide  acceptable  predictions  of  the  heat  release  rate  and  the  time  to 
flashover  (characterized  by  a  measured  heat  release  rate  of  1,000  kW). 

Quintiere’s  model  was  converted  to  Microsoft®  Excel  97,  Visual  Basic  for  Applications  (VBA), 
and  a  validation  was  conducted  to  ensure  that  identical  results  were  obtained  as  with  the  FORTRAN 
version.  All  modifications  to  Quintiere’s  fire  growth  model  were  conducted  on  the  VBA  version  from 
that  point  on. 

In  an  attempt  to  improve  the  heat  release  rate  predictions  by  the  fire  growth  model,  the  material 
properties  were  re-examined.  It  was  determined  that  the  convection  coefficient,  1^,  in  the  Cone 
Calorimeter  is  a  function  of  incident  heat  flux.  To  account  for  this  heat  flux  dependency,  a  new  method 
had  to  be  developed  to  obtain  ignition  properties  from  the  ignition  times  measured  at  different  heat  flux 
levels.  The  heat  of  gasification  values  appeared  to  be  excessive,  and  were  believed  to  be  the  main  cause 
of  the  poor  flashover  time  predictions  with  Quintiere’s  original  model.  Therefore,  the  heat  release  rate 
algorithm  based  on  the  total  pyrolysis  area,  heat  of  combustion,  and  heat  of  gasification  was  abandoned  in 
favor  of  the  use  of  an  exponentially  decaying  heat  release  rate  curve  derived  from  Cone  Calorimeter  data 
at  50  kW/m2. 

The  fire  growth  model  was  also  revised  to  reflect  the  new  properties.  The  equations  used  in  the 
original  FORTRAN  version  were  replaced  by  equations  previously  used  by  Janssens  [Janssens,  1995]  in 
his  QBasic  version  of  the  Quintiere  model.  This  approach  provides  a  more  simplified  method  with  the 
pyrolysis  and  burned  out  areas  represented  by  rectangular  areas  as  opposed  to  complex  trapezoids.  The 
main  modification  in  Janssens’  version  of  the  model  is  an  improved  algorithm  to  provide  a  more  realistic 
characterization  of  the  ignition  burner  flame.  The  heat  flux  to  the  material  in  contact  with  the  burner 
flame  is  determined  based  on  the  heat  output  of  the  burner  and  the  temperature  of  the  material,  as  opposed 
to  simply  using  a  constant  flux.  Additional  modifications  were  made  as  described  below.  The  heat 
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release  rate  from  the  burning  material  was  simulated  using  the  aforementioned  exponentially  decaying 
curve.  As  the  flame  front  progresses,  the  pyrolysis  area  increases.  At  every  incremental  time  step,  a  new 
area  may  ignite  and  start  burning.  The  modified  model  tracks  and  sums  the  heat  release  rate  from  each 
incremental  area  based  on  the  exponentially  decaying  heat  release  rate  curve  to  determine  the  total  heat 
release  rate  from  the  material.  This  method  automatically  accounts  for  burnout,  i.e.,  an  incremental  area 
bums  out  when  its  heat  release  rate  reaches  the  end  of  the  exponential  curve.  An  algorithm  was  added  to 
calculate  the  emissivity  of  the  upper  gas  layer  as  a  function  of  the  specific  extinction  area  measured  in  the 
Cone  Calorimeter  (Quintiere  assumes  the  emissivity  to  be  1.0).  A  routine  was  also  added  to  estimate  the 
smoke  production  rate  in  the  room  based  on  the  specific  extinction  area  measured  in  the  Cone 
Calorimeter. 

The  results  of  the  prediction  using  the  revised  fire  growth  model  are  quite  reasonable.  For 
materials  with  flashover  times  less  than  600  (Material  Nos.  3, 4,  and  8),  the  model  predictions  are  very 
close  to  the  experimental  data.  For  materials  with  flashover  times  between  600  and  1200  (Material 
Nos.  5  and  9),  the  predicted  flashover  times  fall  within  the  same  300-kW  exposure  period.  For  the 
remaining  four  materials,  the  model  consistently  predicts  that  flashover  does  not  occur.  The  model  also 
correctly  predicts  that  the  heat  release  rate  criteria  for  fire  restricting  materials  are  not  exceeded. 

However,  two  of  these  four  materials  (Material  Nos.  1  and  6)  failed  marginally  on  smoke  production  in 
the  tests,  while  the  model  predicts  that  all  four  materials  would  meet  the  smoke  requirements. 

Based  on  the  predictions  achieved,  the  fire  growth  model  presented  in  this  report  represents  a 
simple  method  for  predicting  material  performance  in  the  ISO  9705  Room-Comer  Test  based  on  a  limited 
amount  of  test  requirements  and  input  parameters.  A  single  Cone  Calorimeter  run  at  50  kW/m2, 
supplemented  with  surface  temperature  measurements,  is  all  that  is  needed.  However,  it  is  recommended 
to  perform  at  least  two  or  three  runs  to  improve  the  confidence  in  the  input  data  (and,  hence,  the 
predictions).  More  work  is  needed  to  improve  the  smoke  predictions,  and  to  validate  the  conjecture  that 
the  surface  temperature  at  ignition  can  be  determined  experimentally,  in  lieu  of  being  inferred  from  the 
analysis  of  ignition  data  obtained  at  different  heat  flux  levels. 
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1.0  INTRODUCTION 


This  report  presents  the  results  of  a  program  that  was  conducted  at  Southwest  Research  Institute 
(SwRI)  for  Hughes  Associates  Inc.  in  Baltimore,  MD,  under  subcontract  No.  E00 150-4,  and  covers  part 
of  the  work  to  be  performed  by  Hughes  for  the  U.S.  Coast  Guard  under  Contract  No. 
DTCG39-97-D-E00 1 50. 

The  performance  of  eight  marine  composite  materials  and  one  textile  wall  covering,  previously 
tested  in  accordance  with  ISO  9705  “Fire  tests — full-scale  room  test  for  surface  products”  (also  referred 
to  as  the  room-corner  test),  was  predicted  using  a  fire  growth  model  developed  by  Quintiere  [Quintiere, 
1993],  The  nine  materials  considered  are  presented  in  Table  A-l. 


Table  A-l.  Marine  Composite  Materials. 


Material  No. 

Generic  Name 

1 

FR  Phenolic 

2 

Fire-Restricting  Material 

3 

FR  Polyester 

4 

FR  Vinylester 

5 

FR  Epoxy 

6 

Coated  FR  Epoxy 

7 

Wall  Covering  Material 

8 

Polyester 

9 

FR  Modified  Acrylic 

Originally,  a  FORTRAN  version  of  Quintiere’s  model  was  used  to  predict  the  performance  of  the 
materials.  However,  the  calculations  were  in  poor  agreement  with  the  experimental  data.  Therefore,  the 
model  was  revised  to  include  new  methods  for  determining  ignition  and  flame  spread  properties,  as  well 
as  new  algorithms  for  the  ignition  burner  exposure,  heat  release  and  smoke  production  rate  predictions, 
and  the  emissivity  of  the  upper  layer. 
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2.0  ORIGINAL  MATERIAL  PROPERTIES  AND  MODEL  SIMULATIONS 


2.1  Ignition  Properties 

2.1.1  LIFT  Ignition  Data 

The  ignition  times  obtained  in  the  Lateral  Ignition  and  Flamespread  Test  (LIFT)  apparatus  are 
reported  in  Appendix  E2  of  USCG  R&D  Center  Report  No.  CG-D-22-98  [Janssens,  Garabedian  &  Gray, 
1998].  The  ignition  temperature,  Tig,  and  the  thermal  inertia,  kpc,  were  calculated  from  these  data  using 
the  procedures  outlined  in  American  Society  for  Testing  and  Materials  (ASTM)  Standard  E  1321 
“Standard  Test  Method  for  Determining  Material  Ignition  and  Flamespread  Properties.”  The  results  are 
presented  in  Section  7  of  the  aforementioned  report.  The  ignition  property  values  required  for  input  into 
the  fire  growth  model  are  repeated  in  Table  A-2  below.  Note  that  no  ignition  properties  could  be  obtained 
for  Material  Nos.  1  and  2  because  the  samples  did  not  ignite  at  the  highest  heat  flux  that  can  be  obtained 
in  the  LIFT  apparatus. 

Table  A-2.  Ignition  Properties  Derived  from  LIFT  Ignition  Data. 


Material  No. 

Tig 

(°C) 

kpc 

(kW2s/m4K2) 

1 

— 

— 

2 

— 

— 

3 

375 

1.65 

4 

370 

1.89 

5 

453 

1.73 

6 

643 

8.00 

7 

647 

0.27 

8 

337 

0.74 

9 

385 

1.72 

2.1.2  Cone  Calorimeter  Ignition  Data 

The  nine  materials  were  tested  in  the  Cone  Calorimeter  in  accordance  with  ISO  5660-1:1993  “Fire 
Tests — Reaction  to  fire-Heat  release  rate  of  building  products,”  with  smoke  production  measurements 
made  in  accordance  with  the  draft  version  of  ISO  5660-2.  A  portion  of  the  Cone  Calorimeter  ignition 
data  is  reported  on  the  standard  data  sheets  presented  in  Appendix  C  of  USCG  R&D  Center  Report  No. 
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CG-D-22-98.  Ignition  times  were  obtained  at  additional  heat  flux  levels,  but  these  times  are  not  presented 
in  the  USCG  report.  The  complete  set  of  Cone  Calorimeter  time  to  ignition  data,  including  the  minimum 
heat  flux  necessary  for  ignition,  4o,ig  >  is  presented  in  Attachment  A  of  this  report.  Material  properties 

were  obtained  using  the  same  calculation  procedure  described  in  ASTM  E  1321  as  used  for  the  LIFT  data. 
However,  the  convection  coefficient  of  he  =  15  W/m2-K,  that  is  applicable  to  the  LIFT  apparatus,  was  not 
used.  A  coefficient  applicable  to  the  Cone  Calorimeter  for  specimens  exposed  in  the  horizontal 
orientation  with  the  retainer  frame  was  determined  and  used  for  the  calculations.  This  coefficient  was 
obtained  on  the  basis  of  a  heat  transfer  analysis  of  a  12.7-mm  thick  calcium  silicate  board  exposed  at 
different  heat  flux  levels,  and  instrumented  with  fine  thermocouples  on  the  exposed  and  unexposed 
surfaces.  This  analysis  is  described  in  detail  in  Attachment  B  of  this  report.  Note  that  the  convection 
coefficient  appears  to  be  a  function  of  the  heat  flux.  Therefore,  the  convection  coefficient  evaluated  at  the 
critical  heat  flux  was  used  to  calculate  the  surface  temperature  at  ignition,  Tig.  This  value  was  also  used 
to  calculate  the  thermal  inertia,  kpc,  from  the  slope  b,  which  is  conservative.  The  resulting  material 
ignition  properties  are  listed  in  Table  A-3. 


Table  A-3.  Ignition  Properties  Derived  from  ASTM  E  1354  Cone  Calorimeter  Ignition  Data. 


Material  No. 

T.g 

(°C) 

kpc 

(kW/m-K)2s 

1 

603 

2.23 

2 

632 

1.61 

3 

398 

1.29 

5 

408 

1.83 

6 

518 

1.02 

7 

453 

1.91 

8 

389 

0.83 

9 

408 

2.84 

2.2  Flame  Spread  Properties 


2.2. 1  LIFT  Flame  Spread  Data 

Flamespread  data  were  also  obtained  according  to  the  procedures  described  in  ASTM  Standard 
E  1321.  The  results  are  reported  in  Appendix  E3  of  the  USCG  report.  The  corresponding  flame  spread 
properties  (flame  heating  parameter,  <|>,  and  minimum  temperature  for  lateral  spread,  Ts.  mm)  are  repeated 
in  Table  A-4.  Note  that  for  five  of  the  nine  materials,  lateral  flame  spread  data  could  not  be  obtained 
from  the  LIFT  apparatus. 

Table  A-4.  Lateral  Flame  Spread  Properties  Derived  from  LIFT  Data. 


Material  No. 

<t> 

(kW2/m3) 

Ts.mln 

(°C) 

1 

— 

— 

2 

— 

— 

3 

19.5 

325 

4 

6.3 

428 

5 

— 

’  — 

6 

— 

— 

7 

— 

— 

8 

19.1 

234 

9 

32.9 

307 

2. 2. 2  IMO  Flame  Spread  Data 

A  second  set  of  lateral  flame  spread  properties  was  obtained  using  the  same  calculation  procedure 
as  for  the  LIFT  data,  but  based  on  flame  spread  data  obtained  according  to  International  Maritime 
Organization  (IMO)  Resolution  A.653(16)  “Recommendation  on  Improved  Fire  Test  Procedures  for 
Surface  Flammability  of  Bulkhead,  Ceiling,  and  Deck  Finish  Materials.”  The  flamespread  data  obtained 
from  the  IMO  surface  flammability  test  are  presented  in  Table  A-5.  As  in  the  LIFT,  flame  spread  data 
could  not  be  obtained  for  five  of  the  materials. 
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Table  A-5.  Lateral  Flamespread  Properties  Derived  from  IMO  Resolution  A.653(16)  Test  Data. 


Material  No. 

4> 

(kW2/m3) 

Ts,min 

(°C) 

1 

— 

— 

2 

— 

— 

3  . 

4.8 

406 

4 

18.4 

419 

5 

— 

— 

6 

— 

—  ■ 

7 

— 

— 

8  , 

18.9 

177 

9 

13.6 

284 

2.3  Heat  Release  Rate  Properties 

Representative  values  of  the  effective  heat  of  combustion,  AHc,  and  total  heat  released  per  unit 
area,  Q",  were  obtained  for  each  material.  These  values  were  calculated  as  the  average  of  all  values 
measured  for  that  material  when  tested  in  the  Cone  Calorimeter.  In  addition,  the  heat  of  gasification,  L, 
was  calculated  according  to  a  procedure  recommended  by  Quintiere  [Quintiere,  1993].  The  heat  of 
gasification  values  were  obtained  by  plotting  the  maximum  30-sec  average  heat  release  rate,  HRR3o.max, 
with  respect  to  the  external  heat  flux  from  the  Cone  heater,  q" .  An  effective  heat  of  gasification  was 
calculated  based  on  the  slope  of  a  liner  fit  through  the  data: 


slope 


(Al) 


The  HRR30i  max  values  were  obtained  from  Section  5.4  of  the  USCG  report  and  are  presented  in 
Attachment  C  of  this  report.  The  graphs  of  HRR3o,max  plotted  with  respect  to  q*  and  the  slopes  of  the 

best-fit  lines  used  to  obtain  the  heat  of  gasification  can  also  be  found  in  Attachment  C.  The  resulting  heat 
release  rate  properties  are  presented  in  Table  A-6. 
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Table  A-6.  Heat  Release  Rate  Properties  Derived  from  Cone  Calorimeter  Data. 


Material  No. 

AHC 

(MJ/kg) 

Q" 

(MJ/m2) 

L 

(MJ/kg) 

1 

7.7 

15.9 

9.4 

2 

9.7 

8.1 

16.6 

3 

11.3 

35.6 

17.6 

4 

13.4 

49.0 

13.3 

• 

5 

8.2 

12.6 

19.5 

6 

7.6 

11.7 

15.2 

7 

9.1 

5.8 

8.4 

8 

21.6 

61.5 

6.8 

9 

12.3 

41.1 

13.5 

2.4  Model  Simulations 

2.4.1  FORTRAN  Version 

Simulations  of  the  ISO  9705  room-comer  tests  were  performed  for  the  nine  materials,  using  the 
most  current  FORTRAN  version  of  Quintiere’s  model  that  was  available.  This  version  has  been 
previously  published  by  Dillon  [Dillon,  1998].  Detailed  descriptions  of  the  Quintiere’s  model  and  the 
computer-based  FORTRAN  version  are  available  in  the  literature  [Quintiere,  1993;  Dillon,  1998; 
Janssens,  1995;  Haynes,  1996;  Quintiere,  1995;  Quintiere,  Haynes  &  Rhodes,  1995]  and  will  not  be 
repeated  in  this  report. 

The  FORTRAN  model  requires  the  input  of  the  seven  material  properties  (Tig,  kpc,  (j>,  TS;min,  AHC, 
Q",  and  L),  a  description  of  the  compartment  geometry,  characteristics  of  the  ignition  burner,  as  well  as 
several  modeling  parameters.  A  maximum  of  four  sets  of  property  data  were  developed  using 
combinations  of  the  different  sets  of  data  presented  in  Sections  2.1  through  2.3.  For  the  materials  for 
which  no  flame  spread  data  could  be  obtained  (Material  Nos.  1, 2,  5,  6,  and  7)  a  flame  heating  parameter 
of  <|>  =  0  kW2/m3,  and  a  minimum  surface  temperature  for  spread  of  Ts>min  =  20  °C  were  chosen. 
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Predictions  of  the  heat  release  rate  from  the  materials  were  performed  for  each  set  of  property 
data.  Graphical  comparisons  of  the  measured  and  predicted  heat  release  rate  curves  are  presented  in 
Attachment  D. 

2.4.2  VBA  Version 

After  predicting  the  performance  of  the  materials  using  the  FORTRAN  version  of  Quintiere’s 
model,  the  source  code  was  converted  to  Microsoft®  Excel  97,  Visual  Basic  for  Applications  (VBA).  The 
use  of  VB  A  provides  a  user-friendly  interface  with  which  to  run  the  model.  The  source  code  of  the  model 
is  run  as  a  macro  with  the  material  properties,  room  dimensions,  burner  characteristics,  model  parameters, 
etc.  entered  into  a  windows-based  user  form. 

A  validation  was  conducted  to  ensure  that  identical  results  were  obtained  with  both  the 
FORTRAN  and  VBA  versions.  All  additional  modifications  to  the  model  from  this  point  on  were 
performed  using  the  VBA  version. 


3.0  MODIFIED  MATERIAL  PROPERTIES  AND  MODEL  SIMULATION 

3.1  Ignition  Properties 

Procedures  to  obtain  material  properties  from  piloted  ignition  data  at  different  heat  flux  levels 
commonly  assume  that  the  surface  heat  losses  partly  involve  Newtonian  cooling  which  is  characterized  by 
a  constant  convection  coefficient.  It  is  shown  in  Attachment  B  that  the  convection  coefficient  in  the  Cone 
Calorimeter,  for  specimens  in  the  horizontal  orientation  tested  with  the  retainer  frame,  can  be  expressed 
as  a  linear  function  of  the  external  heat  flux  from  the  Cone  heater: 

hc  s h0  +hjq"  (A2) 

where  h0  =11.8  W/m2-K  and  hi  =  0.00034  1/K  at  heat  flux  levels  below  50  kW/m2,  and  ho  =  25.5  W/m2-K 
and  hi  =  0.000065  1/K  at  heat  flux  levels  equal  to  or  greater  than  50  kW/m2. 

Consider  a  semi-infinite  solid  with  constant  properties  k,  p,  and  c  exposed  to  a  constant  radiant 
heat  flux,  q",  with  radiative  and  convective  heat  losses  from  the  surface: 

0T  ,  02T 
pc —  =  k — r 

51  (A3.a) 


A-ll 


with 
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Equation  (A5.a)  suggests  that  (l/tig)° 5  be  plotted  as  a  function  of  q" .  The  intercept  of  a  straight 
line  fitted  through  the  data  points  is  equal  to  Co/Cu  which  can  be  used  to  determine  Tig.  Once  the  surface 
temperature  at  ignition  is  known,  kpc  can  be  calculated  from  the  slope  of  the  linear  fit.  Because  ho  and  hi 
have  different  values  for  heat  fluxes  below  and  above  50  kW/m2,  the  slope  of  the  linear  fit  is  slightly 
smaller  at  heat  fluxes  below  50  kW/m2.  This  is  illustrated  for  Material  Nos.l  and  2  on  pages  AE-1  and 
AE-2  of  Attachment  E,  respectively.  These  two  materials  did  not  ignite  at  heat  fluxes  below  50  kW/m2, 
and  the  ignition  properties  were  estimated  from  a  correlation  of  ignition  data  at  heat  fluxes  at  50  kW/m2 
and  higher  (solid  line  in  the  graphs).  Based  on  the  resulting  Tigand  kpc  values,  the  location  of  the 
Equation  (A5.a)  was  calculated  for  heat  fluxes  below  50  kW/m2  (dashed  line  in  the  graphs). 

A  different  approach  was  used  for  the  remaining  materials,  all  of  which  ignited  at  heat 
flux  levels  below  50  kW/m2.  Because  the  samples  generally  do  not  exhibit  semi-infinite  solid 
behavior  at  low  heat  flux  levels,  the  data  points  do  not  fall  on  a  straight  line,  and  Equation  (A5.a) 
cannot  be  used  to  determine  Tig  by  extrapolation  [Janssens,  1991].  T;g  can  be  determined 
experimentally,  by  instrumenting  some  of  the  ignition  test  specimens  with  fine  surface 
thermocouples.  This  approach  will  be  recommended  in  Section  4.0  for  future  work,  because  it 
significantly  reduces  the  number  of  small-scale  tests  that  are  needed  to  obtain  the  material 
properties  needed  to  estimate  ISO  9705  (or  other)  room-comer  test  performance.  In  the  absence 
of  experimental  ignition  temperatures,  it  was  assumed  that  the  critical  heat  flux  for  ignition  was  5 
kW/m2  below  the  lowest  heat  flux  level  at  which  ignition  was  observed  within  a  20-min  period 
of  exposure.  This  critical  flux  value  was  used  to  estimate  Tig,  and  kpc  was  then  calculated  from 
the  slope  of  a  line  connecting  the  estimated  critical  heat  flux  on  the  abscissa,  and  the  average  of 
the  data  points  at  50  kW/m2.  The  correlations  are  shown  in  Attachment  E,  pages  AE-3  through 
AE-9.  Ignition  data  for  Material  Nos.  3  through  9,  obtained  at  lower  heat  flux  levels,  are 
presented  in  the  graphs  of  Attachment  E,  however  the  values  are  not  used  in  the  determination  of 
kpc.  The  resulting  ignition  properties  are  summarized  in  Table  A-7. 
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Table  A-7.  Modified  Cone  Calorimeter  Ignition  Properties. 


Material  No. 

Tlg 

(°C) 

kAc 

(kW2s/m4K2) 

1 

450 

1.5252 

2 

380 

2.9935 

3 

337 

1.1547 

4 

337 

1.3283 

5 

393 

1.1169 

6 

483 

0.1615 

7 

412 

0.1973 

8 

323 

0.6061 

9 

349 

1.4393 

3.2  Flame  Spread  Properties 

The  ASTM  E  132 1  LIFT  data  analysis  procedure  specifies  that  l/VV(x)  be  plotted  <  s  a  function 
of  q"  (x),  and  a  straight  line  be  fitted  through  the  data  points.  The  flame  heating  parameter  (j>  is  calculated 
from  the  slope  C  as  follows: 


kpc 

Fh 


(A6) 


Since  the  ignition  properties  have  changed,  the  flame  heating  parameter  was  recalculated.  The 
results  are  presented  in  Table  A-8. 


Table  A-8.  Modified  Flame  Spread  Properties. 


Material  No. 

big 

(W/m-K) 

kAc 

(kW2-s/m4-K2) 

C 

N 

(kW2/m3) 

3 

37.8 

1.1547 

0.22 

16.7 

4 

37.8 

1.3283 

0.42 

5.3 

8 

36.4 

0.6061 

0.16 

17.9 

9 

39.1 

1.4393 

0.26 

13.9 
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3.3  Heat  Release  Rate  Properties 


The  simulation  data  presented  in  Attachment  D  indicate  that  the  original  model  generally 
underestimates  fire  growth.  This  is  attributed,  at  least  in  part,  to  high  heat  of  gasification  values.  To 
eliminate  this  problem,  it  was  decided  to  use  actual  Cone  Calorimeter  heat  release  rate  curves,  instead  of 
heat  release  properties  derived  from  Cone  data.  The  Cone  Calorimeter  data  show  that  heat  flux  effects  are 
not  significant  for  most  materials,  and  the  experimental  data  at  50  kW/m2  were  selected  for  (an  initial) 
analysis.  Quintiere  [Quintiere,  1993]  uses  an  incident  heat  flux  of  60  kW/m2  for  the  ISO  9705  ignition 
burner  flames  and  30  kW/m2  for  a  vertical  wall  flame.  However,  experiments  by  Dillon  [Dillon,  1998] 
and  calculations  by  Janssens  [Janssens,  1995]  indicate  that  a  heat  flux  of  45  to  50  kW/m2  may  be  more 
appropriate  for  the  ISO  9705  burner  with  an  output  of  100  kW.  Therefore,  the  selection  of  Cone 
Calorimeter  data  at  50  kW/m2  is  reasonable  for  this  analysis.  Because  there  were  no  heat  release  data  at 
50  kW/m2  for  Material  Nos.  2  and  5,  the  Cone  Calorimeter  results  at  75  kW/m2  were  used  instead.  The 
average  heat  release  curve  for  all  runs  conducted  at  50  (or  75  kW/m2)  was  approximated  by  an 
exponentially  decaying  function,  as  shown  below: 

0<t-t,8<30  :  Q'  =  HRR,0„„ 

30<t-tig<t„  :  Q-  =  HHR,0,mJ.e-M“'^,0)  (A7) 

t-t„>t,  :  Q"  =  0 

The  decay  parameter  X  is  determined  such  that  the  area  under  the  curve  is  equal  to  the  average 
total  heat  release  rate  measured  in  the  Cone  Calorimeter  experiments.  The  resulting  fitting  parameters  are 
given  in  Table  A-9  and  the  exponential  heat  release  rate  curves  are  presented  in  Attachment  F.  The  idea 
to  use  an  exponentially  decaying  function  was  based  on  earlier  work  by  Magnusson  and  Sundstrom 
[Magnusson  &  Sundstrom,  1985].  The  maximum  30-sec  sliding  average  heat  release  rate  was  used 
instead  of  the  peak,  because  the  former  has  been  proposed  as  one  of  the  criteria  to  qualify  fire-restricting 
materials  on  the  basis  of  Cone  Calorimeter  data. 

The  Cone  Calorimeter  heat  release  rate  curves  presented  in  Attachment  F  represent  an  average  of 
the  measured  heat  release  rates  from  the  tests  conducted.  To  average  the  data,  the  data  was  synchronized 
such  that  the  ignition  times  corresponded  to  the  to  the  average  ignition  time  and  an  average  value  was 
taken  at  each  time  step. 
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Table  A-9.  Heat  Release  Rate  Fitting  Parameters. 


Material  No. 

HRR30,  max 

tb 

8 

(kW/m2) 

(s) 

(1/s) 

1 

30 

241 

0.00630 

2 

32 

265 

0.00179 

3 

109 

720 

0.00194 

4 

124 

978 

0.00150 

5 

80 

431 

0.00430 

6 

25 

50 

0.14000 

7 

39 

143 

0.00000 

8 

350 

355 

0.00474 

9 

118 

1008 

0.00243 

Bold  Italic:  Values  at  75  kW/m2 


3.4  Smoke  Release  Rate  Properties 

A  method  was  developed  to  obtain  a  rough  estimate  of  the  smoke  production  rate  in  the  ISO  9705 
room-comer  test.  The  heat  release  rate  of  the  wall  material  is  divided  by  the  heat  of  combustion  (based 
on  Cone  Calorimeter  data  obtained  at  50  or  75  kW/m2).  The  resulting  mass  loss  rate  is  multiplied  with 
the  specific  extinction  area,  a,  measured  in  the  Cone  calorimeter  to  determine  the  smoke  production  rate. 
The  heat  of  combustion  and  specific  extinction  area  data  are  provided  in  Table  A- 10. 

3.5  Modifications  to  the  Model 

The  VBA  input  data  forms  and  associated  code  were  modified  to  include  the  heat  release  rate  and 
smoke  production  rate  calculations  described  above.  The  geometry  of  the  burning  area  was  simplified 
based  on  equations  previously  used  by  Janssens  [Janssens,  1995]  in  his  QBasic  version  of  the  Quintiere 
model.  The  pyrolysis  and  burned  out  areas  are  represented  by  rectangular  areas  as  opposed  to  complex 
trapezoids.  This  approach  provides  a  simple  means  for  calculating  and  tracking  flamespread,  heat  release, 
and  burnout. 
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Table  A-10.  Heat  of  Combustion  and  Specific  Extinction  Area  Data. 


Material  No. 

)Hc 

(kJ/g) 

(m2/kg) 

1 

4.9 

84 

2 

11.3 

20 

3 

H-4 

594 

4 

12.9 

742 

5 

8.8 

219 

6 

11.4 

391 

7 

6 

48 

8 

21.3 

746 

9 

13 

68 

Bold  Italic:  Values  at  75  kW/nr 


Algorithms  were  also  included  that  better  describe  the  geometry  and  thermal  environment  created 
by  the  ignition  source.  The  heat  flux  to  the  material  in  contact  with  the  burner  flame  is  determined  based 
on  the  heat  output  of  the  burner  and  the  temperature  of  the  material,  as  opposed  to  simply  using  a  constant 
flux. 


Equations  to  estimate  heat  release  rate  on  the  basis  of  the  heat  of  combustion,  heat  of  gasification, 
and  net  heat  flux  were  replaced  with  calculations  on  the  basis  of  the  exponentially  decaying  curve.  The 
total  heat  release  rate  from  the  wall  material  at  time  t  is  given  by 


HRR(t)  =  £(Ai-AM)Q'(t-i) 

i=l 


where  A;  and  Aj.i,are  the  burning  areas  at  time  i  and  i-1  seconds,  respectively. 


(A8) 


As  the  flame  front  progresses,  the  pyrolysis  area  increases.  At  every  incremental  time  step,  a  new 
area  may  ignite  and  start  burning.  The  modified  model  tracks  and  sums  the  heat  release  rate  from  each 
incremental  area  based  on  the  exponentially  decaying  heat  release  rate  curve  to  determine  the  total  heat 
release  rate  from  the  material.  This  method  automatically  accounts  for  burnout,  i.e.,  an  incremental  area 
bums  out  when  its  heat  release  rate  reaches  the  end  of  the  exponential  curve. 


Routines  were  also  added  to  estimate  the  smoke  production  rate,  SPR,  and  the  corresponding 
emissivity  of  the  hot  gas  layer,  eg.  The  smoke  production  rate  is  calculated  as  the  heat  release  rate  divided 
by  the  heat  of  combustion  and  multiplied  by  the  specific  extinction  area  measured  in  the  Cone 
Calorimeter. 


The  emissivity  of  the  upper  gas  layer  is  also  calculated  as  a  function  of  the  specific  extinction 
area.  Quintiere’s  original  model  assumes  eg  =  1,  which  appears  to  be  overly  conservative  for  low  smoke 
producing  materials.  The  emissivity  of  a  hot  and  smoky  gas  layer  can  be  estimated  from  [Seader  & 
Einhorn,  1976] 

eg  =  1  -  exp(-  (0.33  +  0.47C,  )(H  -  Z; ))  (A9) 

where  Cs  is  the  concentration  of  soot  particles  (g/m3),  H  is  the  room  height  (m),  and  Z;  is  the  layer 
interface  height  (m).  The  soot  concentration  can  be  estimated  from 


SPR  aHRR(t) 
kmV0  '  kmV0AHc 


(A10) 


where  km  =  7.6  m2/g  based  on  data  by  Seader  and  Einhorn  for  flaming  wood  and  plastics  fires  [Seader  & 
Einhorn,  1976],  The  layer  depth  was  set  equal  to  1  m,  based  on  observations  in  the  ISO  9705  tests.  To 
obtain  a  conservative  (high)  estimate  of  Sg,  a  volumetric  vent  flow  of  0.5  m3/s  was  chosen.  The  resulting 
eg  estimates  vary  between  0.4  and  1.0. 


3.6  Modified  Model  Simulations 


The  results  of  the  simulations  with  the  modified  model  are  presented  in  Attachment  G.  A 
comparison  is  provided  between  the  calculated  and  measured  heat  release  rate,  upper  gas  layer 
temperature,  smoke  production  rate,  and  heat  flux  to  the  floor.  Comparisons  with  the  experimental 
temperature  data  are  based  on  the  readings  of  the  ceiling  thermocouple  opposite  the  burner.  The  floor 
flux  was  calculated  on  the  basis  of  radiation  from  the  upper  gas  layer,  i.e.,  the  heat  flux  from  the  flames 
was  ignored.  As  a  result,  the  model  underestimates  the  floor  flux  for  materials  with  low  heat  release  rate. 


The  heat  release  predictions  are  quite  reasonable.  For  the  materials  with  flashover  times  less  than 
600  s  (Nos.  3, 4,  8),  the  model  predictions  are  very  close  to  the  measured  heat  release  rate.  For  materials 
with  flashover  times  between  600  and  1200  s  (Nos.  5  and  9),  the  predicted  flashover  times  fall  within  the 
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same  300-kW  exposure  period.  For  the  remaining  materials,  the  model  correctly  predicts  that  flashover 
does  not  occur.  Measured  and  predicted  flashover  times  (times  to  1  MW  total  heat  release  rate,  except 
where  noted)  are  compared  in  Table  A-l  1.  The  measured  and  predicted  average  heat  release  and  smoke 
production  rates  used  to  qualify  fire  restricting  materials  are  presented  in  Table  A-12.  For  the  materials 
that  did  not  flash  over,  the  model  correctly  predicts  that  the  heat  release  criteria  for  fire  restricting 
materials  are  not  exceeded.  However,  two  of  the  four  materials  (No.  1  and  No.  6)  failed  marginally  on 
smoke,  while  the  model  predicts  that  all  four  materials  would  meet  the  smoke  requirements. 

Table  A-ll.  Comparison  of  Measured  and  Predicted  Flashover  Times. 


Material  No. 

Experimental  (s) 

Model  (s) 

3 

342 

345 

4 

306 

305 

5 

978 

666 

8 

102 

56 

9 

672 

611 

Bold  Italic:  Time  to  750  kW 


3.7  Conclusions 

The  modified  model  does  a  reasonable  job  in  predicting  flashover  times  and  heat  release  rates. 
However,  the  model  underestimates  smoke  production,  and  more  work  is  needed  to  address  this  problem. 


Table  A-12.  Measured  and  Predicted  Average  Heat  Release  and  Smoke  Production  Rates. 


Material 

No. 

Experimental  Data 

Model  Predictions 

HRRavg 

(kW) 

SPRavg 

(m2/s) 

SPR60s 

(m2/s) 

HRRavg 

(kW) 

HRR30S 

(kW) 

SPRavg 

(m2/s) 

SPR60s 

(m2/s) 

1 

62 

159 

1.5 

■a 

36 

120 

0.6 

2.0 

2 

31 

129 

0.2 

0.5 

47 

134 

0.1 

0.2 

6 

134 

1.5 

3.5 

8 

193 

0.3 

4.5 

7 

17 

131 

0.1 

0.2 

48 

194 

0.4 

1.5 

Criteria 

<  100 

<500 

<  1.4 

<8.3 

<  100 

<500 

<  1.4 

<8.3 

Bold  Italic:  Marginal 


failure 
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4.0  SUMMARY  AND  RECOMMENDATIONS 


A  modified  version  of  Quintiere’s  fire  growth  model  was  developed  to  predict  performance  of 
wall  and  ceiling  linings  in  the  ISO  9705  room-comer  test.  Good  quantitative  agreement  was  found 
between  the  predicted  and  measured  heat  release  rates  for  a  set  of  eight  marine  composite  materials  and 
one  textile  wall  covering. 

The  results  of  this  analysis  indicate  that  a  minimal  amount  of  small-scale  data  are  required  to  run 
the  model.  A  single  Cone  Calorimeter  test  at  50  kW/m2  is  all  that  is  needed.  However,  it  is 
recommended  to  perform  at  least  two  or  three  runs  to  improve  the  confidence  in  the  input  data  (and, 
hence,  the  predictions).  The  choice  of  using  Cone  data  at  50  kW/m2  is  based  on  a  maximum  heat  flux  to 
the  walls  of  approximately  45  to  50  kW/m2  from  the  flames  of  the  ISO  burner  at  100  kW.  However,  the 
heat  flux,  from  the  burner  with  an  output  of  300  kW  can  range  between  50  kW/m2  to  approximately  80  to 
90  kW/m2  near  the  wall-ceiling  interface  and  approximately  100  kW/m2  to  the  ceiling  just  above  the 
burner  (Dillon,  1998).  Therefore,  after  10  minutes  it  would  seem  reasonable  to  use  Cone  Calorimeter 
data  at  a  higher  incident  heat  flux.  Only  Material  Nos.  5  and  9  produced  flashover  conditions  after  10 
minutes,  but  the  model  conservatively  predicted  flashover  faster  than  the  experiments.  This  indicates  that 
the  simple  model  presented  in  this  report  may  be  reasonable  for  predicting  flashover  for  a  variety  of 
materials.  However,  further  analysis  of  this  of  scenario  is  clearly  warranted. 

Although  the  flamespread  data  were  developed,  it  was  found  that  lateral  flamespread  was 
insignificant  to  the  predictions  of  material  performance  in  the  ISO  9705  test.  Therefore,  it  is  suggested 
that  LIFT  and  IMO  surface  flammability  test  data  are  unnecessary  for  performing  predictions  by  the 
method  presented  in  this  report.  This  can  be  advantageous  due  to  the  high  cost  of  lateral  flamespread 
tests  in  comparison  with  Cone  Calorimeter  tests. 

A  first  attempt  was  made  at  estimating  the  smoke  production  rate  on  the  basis  of  the  specific 
extinction  area  measured  in  the  Cone  Calorimeter.  The  smoke  predictions  are  in  qualitative  agreement 
with  the  measurement,  but  in  some  cases  they  err  on  the  unconservative  side.  Therefore,  additional  work 
is  needed  to  improve  the  smoke  predictions. 

It  has  been  suggested  that  the  critical  heat  flux  for  ignition  can  be  determined  on  the  basis  of 
supplementary  surface  temperature  measurements.  However,  since  the  tests  conducted  were  not  provided 
with  surface  mounted  thermocouples,  the  materials  should  be  re-tested  to  verify  this  conjecture. 
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Original  Cone  Calorimeter  Ignition  Properties 
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Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  3 


q'ig  =  18  kW/m2 


•  ft 

qe 

(kW/m2) 

*ig 

(s) 

<ig 

t-  1/2 
>> 

100 

14 

0.180 

3.742 

95 

16 

0.189 

4.000 

90 

17 

0.200 

4.123 

75 

25 

0.240 

5.000 

75 

28 

0.240 

5.292 

50 

71 

0.360 

8.426 

50 

59 

0.360 

7.681 

25 

248 

0.720 

15.748 

25 

250 

0.720 

15.811 

18 

397 

1.000 

19.925 

17 

NI 

1.059 

N/A 

15 

NI 

1.200 

N/A 

MATERIAL  #3:  IGNITION  DATA 
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Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  4 


=  18  kW/mJ 


< 

(kW/m2) 

*ig 

(*) 

<«g 

C 

t  1/2 
>) 

100 

13 

0.180 

3.606 

95 

14 

0.189 

3.742 

90 

15 

0.200 

3.873 

75 

35 

0.240 

5.916 

75 

32 

0.240 

5.657 

50 

79 

0.360 

. 

8.888 

50 

70 

0.360 

8.367 

25 

265 

0.720 

16.279 

25 

355 

0.720 

18.841 

20 

667 

0.900 

25.826 

18 

735 

1.000 

27.111 

17 

NI 

1.059 

N/A 

15 

NI 

1.200 

N/A 

MATERIAL  #4:  IGNITION  DATA 


Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  5 


1.2 


q;ig=19kW/m’ 


• » 

qe 

(kW/m2) 

t|g 

(s) 

<<* 

c 

t  1/2 
>) 

100 

36 

0.190 

6.000 

95 

38 

0.200 

6.164 

90  . 

39 

0.211 

6.245 

85 

49 

0.224 

7.000 

75 

58 

0.253 

7.616 

75 

60 

0.253 

7.746 

50 

135 

0.380 

11.619 

50 

110 

0.380 

10.488 

35 

134 

0.543 

11.576 

30 

205 

0.633 

14.318 

25 

349 

0.760 

18.682 

23 

930 

0.826 

30.496 

20 

407 

0.950 

20.174 

18 

NI 

1.056 

N/A 

15 

NI 

1.267 

N/A 

MATERIAL  #5:  IGNITION  DATA 


Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  6 


q'ig  =  33  kW/m2 


c 

(kW/m2) 

*ig 

(») 

Cig 

Qe 

t  1/2 
>) 

100 

22 

0.330 

4.690 

100 

20 

0.330 

4.472 

95 

19 

0.347 

4.359 

90 

21 

0.367 

4.583 

85 

27 

0.388 

5.196 

75 

30 

0.440 

5.477 

57 

57 

0.579 

7.550 

50 

68 

0.660 

8.246 

45 

71 

0.733 

8.426 

40 

84 

0.825 

9.165 

38 

122 

0.868 

11.045 

36 

205 

0.917 

14.318 

34 

218 

0.971 

14.765 

32 

NI 

1.031 

N/A 

25 

NI 

1.320 

N/A 

MATERAL#6:  IGNmON  DATA 
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Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  7 


“  24  kW/m2 


.  H 

qe 

(kW/m2) 

(s) 

q'ig 

Qe 

t  1/2 

<v> 

75 

14 

0.320 

3.742 

75 

13 

0.320 

3.606 

70 

14 

0.343 

3.742 

65 

16 

0.369 

4.000 

60 

18 

0.400 

4.243 

50 

26 

0.480 

5.099 

50 

29 

0.480 

5.385 

25 

1165 

0.960 

34.132 

25 

690 

0.960 

26.268 

23 

NI 

1.043 

N/A 

22 

NI 

1.091 

N/A 

20 

NI 

1.200 

N/A 

MATERIAL  #7:  IGNITION  DATA 


Note:  X  indicates  data  points  that  were  neglected  in  the  determination  of  the  best-fit  line. 


Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  8 


q'  ,lg  =  17  kw/m2 


.  H 

qe 

(kW/m2) 

*ig 

(s) 

q' 

t  1/2 
>) 

85 

11 

0.200 

3.317 

80 

11 

0.213 

3.317 

75 

16 

0.227 

4.000 

75 

15 

0.227 

3.873 

70 

13 

0.243 

3.606 

65 

14 

0.262 

3.742 

50 

33 

0.340 

5.745 

50 

26 

0.340 

5.099 

25 

100 

0.680 

10.000 

25 

145 

0.680 

12.042 

20 

161 

0.850 

12.689 

18 

381 

0.944 

19.519 

17 

P80 

1.000 

29.665 

16 

NI 

1.063 

N/A 

15 

NI 

1.133 

N/A 

MATERIAL  #8:  IGNITION  DATA 


Note:  x  indicates  a  data  point  that  was  neglected  in  the  determination  of  the  best-fit  line. 


Original  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  9 


1.2 


<L,ig  =  19  kW/m2 


•  n 

<le 

(kW/m2) 

(s) 

t-  1/2 
<?“> 

100 

20 

0.190 

4.472 

95 

24 

0.200 

4.899 

90 

24 

0.211 

4.899 

85 

43 

0.224 

6.557 

75 

65 

0.253 

8.062 

75 

55 

0.253 

7.416 

75 

65 

0.253 

8.062 

50 

85 

0.380 

9.220 

50 

105 

0.380 

10.247 

50 

90 

0.380 

9.487 

25 

540 

0.760 

23.238 

25 

420 

0.760 

20.494 

25 

425 

0.760 

20.616 

20 

998 

0.950 

31.591 

19 

1104 

1.000 

33.226 

18 

NI 

1.056 

N/A 

MATERIAL  #9:  IGNITION  DATA 
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Cone  Calorimeter  Convective  Heat  Transfer  Coefficient  Determination 

(Consisting  of  6  pages) 
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Cone  Calorimeter  Convection  Heat  Transfer  Coefficient  Determination 


The  convective  heat  transfer  coefficient,  fr.,  for  materials  tested  horizontally  in  the  Cone 
Calorimeter  with  the  sample  retainer  frame  was  determined  by  exposing  a  piece  of  calcium  silicate 
(Ca-Si)  board  to  the  radiant  flux  of  the  cone  heater  and  measuring  the  temperatures  of  the  exposed  and 
unexposed  surfaces.  The  temperatures  of  the  two  surfaces  were  used  to  calculate  the  net  heat  flux  into  the 
calcium  silicate  board  by  a  finite  difference  method  using  the  material  properties  presented  in  Table 
AB-1 .  The  net  heat  flux  was  then  used  to  calculate  the  convective  heat  transfer  coefficient  based  on  the 
following  expression: 


h  = 


sC-sq  (Ts4-t)-q;et 
(TS-TJ 


103 


(AB-1) 


where  he  =  convective  heat  transfer  coefficient,  W/m2-K 
q*  =  external  heat  flux  from  the  cone  heater,  kW/m2 


£  =  surface  emissivity  of  the  test  specimen 

a  =  Stefan-Boltzmann  constant,  5.67  x  10'n  kW/m2-K4 

Ts  =  surface  temperature,  K 

To,  =  ambient  temperature,  K 

q'e,  =  net  heat  flux  into  the  specimen,  kW/m2 


Table  AB-1.  Properties  of  Marinite  I  Calcium  Silicate  Board  (BNZ  Materials,  Inc.). 


Thermal  Conductivity,  k 

0.12  W/m-K 

Density,  p 

737  kg/m3 

Specific  Heat,  c 

1.17  kJ/kg  °C  @  93°C 

1.25  kJ/kg  °C  @  205°C 
1.34  kJ/kg°C  @  316°C 
1.42  kJ/kg-°C  @  425°C 

The  test  specimen  was  a  100  x  100  mm  (4x4  in.)  x  12.7  mm  (1/2  in.)  thick  piece  of  “Marinite  F’ 
Ca-Si  board  obtained  from  BNZ  Materials,  Inc.  Prior  to  testing,  the  specimen  was  heated  to  200°C  in  a 
muffle  furnace  for  approximately  12  hours  to  remove  any  trapped  moisture  and  stored  with  drierite  to 
maintain  a  moisture  free  specimen.  A  10-mil  diameter,  Type  K,  thermocouple  (TC)  was  used  to  measure 
the  exposed  surface  temperature.  The  TC  was  fixed  to  the  board  by  drilling  small  holes  through  the  board 
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Cone  Calorimeter  Convection  Heat  Transfer  Coefficient  Determination 


approximately  13  mm  on  either  side  of  the  center  (along  the  diagonal)  and  running  the  TC  wires  through 
the  holes  with  the  TC  bead  securely  against  the  sample  surface.  A  10-mil  diameter  Inconel  sheathed 
thermocouple  was  attached  to  the  bottom  surface  to  measure  the  unexposed  surface  temperature.  The  test 
specimen  was  backed  with  refractory  blanket  and  placed  in  a  standard  Cone  Calorimeter  sample  holder 
with  the  sample  retainer  frame.  The  specimen  was  exposed  to  radiant  heat  fluxes  of  10  to  100  kW/m2  in  a 
horizontal  orientation  and  the  surface  temperatures  were  measured  and  recorded  every  0.5  sec  for 
approximately  25  to  40  min. 

In  order  to  calculate  the  convective  heat  transfer  coefficient  using  Equation  (AB-1),  the 
surface  emissivity,  £,  of  the  Ca-Si  board  was  required.  The  emissivity  was  not  available  in  the 
material  literature,  so  it  was  determined  based  on  the  measured  surface  temperature.  After 
heating  the  surface  of  the  calcium  silicate  board  to  approximately  700°C  in  the  Cone 
Calorimeter,  the  specimen  was  removed  and  allowed  to  cool,  in  the  sample  holder,  in  a  vertical 
orientation.  The  surface  temperature  at  the  center  of  the  board  was  measured  for  over  35 
minutes  using  the  surface  TC  as  well  as  an  infrared  pyrometer  set  to  an  emmissivity  of  £  =  1 .0. 
The  surface  emissivity  was  calculated  based  on  the  two  measured  temperatures.  The  surface 
temperatures  of  the  Ca-Si  board  and  the  calculated  surface  emissivity  values  are  plotted  with 
respect  to  time  in  Figure  AB-1.  The  figure  indicates  that  the  average  surface  emissivity  for 
Marinite  I  calcium  silicate  board  is  0.88  over  a  range  of  40  to  700°C. 

The  measured  temperatures  of  the  exposed  and  unexposed  surfaces  of  the  calcium  silicate  board 
were  used  to  calculate  the  net  heat  flux  and  the  associated  convective  heat  transfer  coefficient.  The 
experimentally  determined  surface  emissivity  was  used  along  with  q'et  in  Equation  (AB-1),  to  calculate 

he  for  horizontal  specimens  as  a  function  of  time  over  a  range  of  external  heat  flux  levels:  10  to  100 
kW/m2.  An  example  of  the  result  of  such  a  calculation  is  presented  graphically  in  Figure  AB-2.  Average 
he  values  were  determined  for  all  the  heat  flux  exposures  and  are  presented  numerically  in  Table  B-2  and 
graphically  with  respect  to  the  external  heat  flux  in  Figure  AB-3.  Based  on  the  data  presented  in  Figure 
AB-3,  an  exponential  expression  for  the  convective  heat  transfer  coefficient  was  developed: 

hc  =  6.56-q'035 
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Cone  Calorimeter  Convection  Heat  Transfer  Coefficient  Determination 


Figure  AB-1.  Emissivity,  e,  of  Marinite  I,  Calcium  Silicate  Board  (BNZ  Materials,  Inc.). 

The  critical  heat  flux  for  ignition,  q'  ig ,  was  used  in  Equation  (AB-2)  to  calculate  the  correct  1^ 
value  in  the  determination  of  the  temperature  for  material  ignition,  Tjg,  in  Section  2.1.2. 


0  300  600  900  1200  1500  1800 

Time  (s) 

Figure  AB-2.  Example  of  the  Convective  Heat  Transfer  Coefficient,  hc,  in  the  Cone 
Calorimeter  Calculated  with  Respect  to  Time:  q'  =  10  kW/m2. 
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Table  AB-2.  Convective  Heat  Transfer  Coefficients  for  Horizontal  Specimens  with  the 
Sample  Retainer  Frame  in  SwRI’s  Cone  Calorimeter. 


Figure  AB-3:  Convective  Heat  Transfer  Coefficient,  h„  Plotted  with  Respect  to  the  External  Flux 
for  Horizontal  Specimens  Tested  in  the  Cone  Calorimeter  with  the  Sample  Retainer 
Frame. 
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Cone  Calorimeter  Convection  Heat  Transfer  Coefficient  Determination 


In  order  to  calculate  the  modified  ignition  properties  in  Section  3.1,  the  exponential  expression 
determined  in  Figure  AB-3  was  simplified  to  produce  two  linear  expressions  for  he  based  on  the  external 
heat  flux  (see  Figure  AB-4): 


hc  =  0.34  -qj  +  l  1.8 
hc  =  0.065 -q'+ 25.5 


q'  <  50  kW/m2 

q'  >  50  kW/m2  (AB-3) 


Figure  AB-4:  Convective  Heat  Transfer  Coefficient,  hc,  with  Respect  to  the  External  Flux  for  Horizontal 
Specimens  Tested  in  the  Cone  Calorimeter  with  the  Sample  Retainer  Frame. 
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HRRao  ,  max  (kW/m1) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  1 


.  tt 

qe 

(kW/m2) 

HRR30,  max 
(kW/m2) 

50 

30 

75 

76 

75 

63 

100 

60 

100 
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MATERIAL#!:  HEAT  OF  GASIFICATION 


0  10  20  30  40  50  60  70  80  90  100  110  120 

External  Heat  Flux  (kW/m2) 


HRR30,  max  (kW/m2) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  2 


■  ff 

qe 

(kW/mJ) 

HRRjo,  max 
(kW/m2) 

75 

36 

75 

27 

100 

43 

100 

49 

MATERIAL  #2:  HEAT  OF  GAS  IFICATION 


0  10  20  30  40  50  60  70  80  90  100  110  120 

External  Heat  Flux  (kW/m2) 


HRR?o,  ma*  (kW/m2) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  3 


•  ff 

qe 

(kW/m2) 

HRRjg^  max 

(kW/m2) 

25 

93 

25 

94 

50 

113 

50 

104 

75 

125 

75 

126 

MATERIAL #3:  HEAT  OF  GASIFICATION 


0  10  20  30  40  50  60  70  80  90  100  110  120 

External  Heat  Flux  (kW/m2) 


HRRw  ,  max  (kW/m1) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  4 


.  n 

qe 

(kW/m2) 

HRR30,  mai 
(kW/mJ) 

25 

89 

25 

104 

50 

121 

50 

126 

75 

154 

75 

140 

MATERIAL #4:  HEAT  OF  GASIFICATION 


0  10  20  30  40  50  60  70  80  90  100  110  120 

External  Heat  Flux  (kW/m2) 
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MATERIAL  NO.  5 


■  H 

qe 

(kW/m2) 

HRR30,  max 

(kW/m2) 

50 

54 

50 

75 

75 

85 

75 

74 

100 

84 

MATERIAL #5:  HEAT  OF  GASIFICATION 


External  Heat  Flux  (kW/m2) 


HRRjo  ,  max  (kW/m2) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  6 


.  If 

qe 

(kW/m2) 

HRR30,  max 

(kW/m2) 

50 

27 

75 

77 

75 

15 

100 

53 

100 

54 

MATERIAL #6:  HEAT  OF  GASIFICATION 


0  10  20  30  40  50  60  70  80  90  100  110  120 

External  Heat  Flux  (kW/m2) 


I 


I 


I 


TTWK 


qe 

(kW/m2) 

HRR30,  max 

(kW/m2) 

25 

38 

25 

34 

50 

58 

50 

65 

75 

46 

75 

134 

Slope  =  01.080 


HRRio.  max  (kW/m1) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 


MATERIAL  NO.  8 


.  /f 

qe 

(kW/m2) 

HRR30,  max 
(kW/m2) 

25 

245 

25 

288 

50 

370 

50 

330 

75 

451 

75 

399 

MATERIAL #8:  HEAT  OF  GASIFICATION 
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HRR30,  max  (kW/m1) 


Original  Cone  Calorimeter  Heat  of  Gasification  Data 
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•  n 

qe 

(kW/m2) 

HRR30,  max 
(kW/m2) 

25 

87 

25 

91 

25 

103 

50 

106 

50 

120 

50 

129 

75 

127 

75 

152 

75 

139 

MATERIAL  #9:  HEAT  OF  GASIFICATION 
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Original  Model  Predictions 


Predictions  of  the  performance  of  the  nine  USCG  materials  in  the  ISO  9705  Room-Comer  test 
were  made  using  the  material  properties  presented  in  Section  2.0  input  into  Quintiere’s  fire  growth  model. 
Predictions  were  made  using  the  ignition  data  from  the  Cone  Calorimeter  and  the  Lateral  Ignition  and 
Flamespread  Test  (LEFT)  apparatus  and  flamespread  data  from  the  LIFT  apparatus  and  the  International 
Maritime  Organization  (IMO)  Surface  Flammability  Test. 

The  predictions  for  the  USCG  materials  are  presented  below.  The  legend  of  the  graphical 
presentations  of  the  material  performance  indicates  the  material  properties  that  were  used  for  each 
prediction.  The  first  term  indicates  the  ignition  properties  and  the  second  term  indicates  the  flame  spread 
properties.  For  example  “Cone  &  IMO  (Cl)”  indicates  that  the  ignition  data  from  the  Cone  Calorimeter 
and  flamespread  data  from  the  IMO  Surface  Flammability  Test  were  used.  For  materials  where  no 
flamespread  data  were  available,  only  the  set  of  ignition  data  used  is  presented.  For  example  “Cone” 
simply  indicates  that  there  were  no  flame  spread  data  available  and  the  Cone  Calorimeter  ignition  data 
were  used.  For  such  a  material  the  flame  heating  parameter,  O,  was  input  as  0  kW2/m3  and  the  minimum 
temperature  for  lateral  flame  spread,  Ts  min,  was  input  as  20°C  (293  K),  in  order  to  indicate  that  no 
flamespread  should  be  calculated  by  the  prediction  model. 
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Heat  Release  Rate  (kW)  Heat  Release  Rate  (kW) 


Original  Model  Predictions 


ISO  9705  ROOM-CORNER  TEST  PREDICTIONS:  MATERIAL  #3 
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ISO  9705  ROOM-CORNER  TEST  PREDICTIONS:  MATERIAL  #4 
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ISO  9705  ROOM-CORNER  TEST  PREDICTIONS:  MATERIAL  #5 
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ISO  9705  ROOM-CORNER  TEST  PREDICTIONS:  MATERIAL  #7 
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Modified  Cone  Calorimeter  Ignition  Properties 

MATERIAL  NO.  1 


Heat  Flux 
(kW/m2) 

*ig 

(S) 

(1/ti/5 

(s^) 

(kpc/tjg)05 

(k\V/m2K) 

Fit 

75 

77 

0.1140 

0.1378 

0.1116 

75 

78 

0.1132 

0.1378 

0.1116 

65 

140 

0.0845 

0.1086 

0.0879 

60 

210 

0.0690 

0.0940 

0.0761 

55 

211 

0.0688 

0.0794 

0.0643 

50 

324 

0.0556 

0.0648 

0.0524 

50 

380 

0.0513 

0.0648 

0.0524 

27.84 

0.0000 

0.0000 

0.0000 

e  =  0.92 

Intercept  (>  50)  =  27.84  kW/m2 
he  =  27.3  W/m2-K 
hj„=  59.6  W/m2-K 
Tig=  450  °C 

kpc  =  1.5252  kW2-s/m4K2 
Intercept  (<  50)  =  24.48  kW/m2 
Minimum  Flux  =  47.0  kW/m2 

MATERIAL  NO.  1 


Heat  Flux  (kW/m2) 


Modified  Cone  Calorimeter  Ignition  Properties 
MATERIAL  NO.  2 


Heat  Flux 
(kW/m2) 

(s) 

<l/tj/5 

(kpc/tjg)®'5 

(kW/m2-K) 

Fit 

75 

72 

0.1179 

0.1915 

0.1107 

75 

83 

0.1098 

0.1915 

0.1107 

70 

107 

0.0967 

0.1740 

0.1006 

65 

138 

0.0851 

0.1565 

0.0904 

60 

163 

0.0783 

0.1389 

55 

192 

0.0722 

0.1214 

54 

198 

0.0711 

0.1179 

20.40 

0.0000 

0.0000 

0.0000 

£  = 

0.92 

Intercept  (>  50)  = 

20.4  kW/m2 

K  = 

26.8  W/m2-K 

hj8  — 

52.1  W/m2-K 

TiB  = 

380  °C 

kpc  = 

2.9935  kW2-s/m4-K2 

Intercept  (<  50)  = 

16.74  kW/m2 

Minimum  Flux  = 

53.0  kW/m2 

MATERIAL  NO.  2 


s 


o 
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Modified  Cone  Calorimeter  Ignition  Properties 

MATERIAL  NO.  3 


Heat  Flux 
(kW/m2) 

^»g 

(s) 

(1/ti/5 

(kpc/tig)0-5 

(kW/m2K) 

Fit 

50 

71 

0.1187 

0.1337 

0.1244 

50 

59 

0.1302 

0.1337 

0.1244 

25 

248 

0.0635 

25 

250 

0.0632 

18 

397 

0.0502 

13 

0.0000 

0.0000 

0.0000 

e  =  0.92 

Critical  Flux  =  13.0  kW/m2 
he  =  16.2  W/m2-K 
hi,  =  37.8  W/m2K 
Tj(.=  337  °C 

kpc  =  1.1547  kW2-s/m4-K2 
Minimum  Flux  =  17.0  kW/m2 


MATERIAL  NO.  3 


Heat  Flux  (kW/m2) 
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Modified  Cone  Calorimeter  Ignition  Properties 
MATERIAL  NO.  4 


Heat  Flux 

*'g 

(kpc/tjg)0'5 

Fit 

(kW/m2) 

(s) 

(s^5) 

(kW/m2-K) 

50 

79 

0.1125 

0.1337 

0.1160 

50 

70 

0.1195 

0.1337 

0.1160 

25 

265 

0.0614 

25 

355 

0.0531 

20 

667 

0.0387 

18 

397 

0.0502 

13 

0.0000 

0.0000 

0.0000 

e  =  0.92 

Critical  Flux  =  13.0  kW/m2 
h«=  16.2  W/m2-K 
hie=  37.8  W/m2-K 
TiR=  337  °C 

kpc=  1.3283  kW2-s/m4-K2 
Minimum  Flux  =  17.0  kW/m2 


0.09 


MATERIAL  NO.  4 


Heat  Flux  (kW/m2) 


Modified  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  5 


Heat  Flux 
(kW/m2) 

*ig 

(s) 

(1/ti/-5 

<o 

(kpc/tjg)0-5 

(kW/m2K) 

Fit 

50 

135 

0.0861 

0.0959 

0.0907 

50 

110 

0.0953 

0.0959 

0.0907 

35 

134 

0.0864 

30 

205 

0.0698 

25 

349 

0.0535 

23 

930 

0.0328 

18 

0.0000 

0.0000 

0.0000 

e  =  0.92 

Critical  Flux  =  18.0  kW/m2 
h<.=  17.9  W/m2-K 
hiR=  44.4  W/m2K 
Tig=  393  °C 

kpc  =  1.1169  kW2s/m4-K2 
Minimum  Flux  =  22.0  kW/m2 


MATERIAL  NO.  5 


=?  0.09 


Heat  Flux  (kW/m2) 
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Modified  Cone  Calorimeter  Ignition  Properties 

MATERIAL  NO.  6 


Heat  Flux 
(kW/m2) 


50  68  0.1213 

45  71  0.1187 

40  84  0.1091 

38  122  0.0905 

36  205  0.0698 

34  218  0.0677 

29  0.0000 


e  =  0.92 

Critical  Flux  =  29.0  kW/m2 
hc=  21.7  W/m2-K 
hig  =  57.6  W/m2-K 
Tjg=  483  °C 

kpc=  0.1615  kW2-s/m4K2 
Minimum  Flux  =  32.0  kW/m2 


(kpc/t|g)0'5 

(kW/m2K) 

Fit 

0.0487 

0.1213 

0.0000 

0.0000 

MATERIAL  NO.  6 


Heat  Flux  (kW/m2) 


Modified  Cone  Calorimeter  Ignition  Properties 
MATERIAL  NO.  7 


Heat  Flux 
(kW/m2) 


(i/y0-5 

(sJs) 

(kpc/tjg)0'5 

(kW/m2K) 

Fit 

0.1961 

0.0848 

0.1909 

0.1857 

0.0848 

0.1909 

0.0293 

0.0381 

0.0000 

0.0000 

0.0000 

9  0.15 

Vi 


M 

4=  0.10 


e  =  0.92 

Critical  Flux  =  20.0  kW/m2 
1^  =  18.6  W/m2-K 
his=  46.9  W/m2-K 
Ti*=  412  °C 

kpc=  0.1973  kW2-s/m4-K2 
Minimum  Flux  =  23.0  kW/m2 


MATERIAL  NO.  7 


Heat  Flux  (kW/m2) 


Modified  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  8 


Heat  Flux 

*ig 

(kpc/tjg)0'5 

Fit 

(kW/m2) 

(s) 

m 

(kW/m2K) 

50 

33 

0.1741 

0.1441 

0.1851 

50 

26 

0.1961 

0.1441 

0.1851 

25 

100 

0.1000 

25 

145 

0.0830 

20 

161 

0.0788 

18 

381 

0.0512 

17 

880 

0.0337 

12 

0.0000 

0.0000 

0.0000 

E  =  0.92 

Critical  Flux  =  12.0  kW/m2 
hc=  15.9  W/m2K 
hjg=  36.4  W/m2K 
Ti(,=  323  °C 

kpc=  0.6061  kW2-s/m4-K2 
Minimum  Flux  =  16.0  kW/m2 


Heat  Flux  (kW/m2) 


S)  c.n(3,J/l) 


Modified  Cone  Calorimeter  Ignition  Properties 


MATERIAL  NO.  9 


Heat  Flux 
(kW/m2) 

tig 

(s) 

(1/tj/5 

(kpc/tig)0'5 

(kW/m2-K) 

50 

85 

0.1085 

0.1246 

50 

105 

0.0976 

0.1246 

50 

90 

0.1054 

0.1246 

25 

540 

0.0430 

25 

420 

0.0488 

25 

425 

0.0485 

20 

998 

0.0317 

19 

1104 

0.0301 

14 

0.0000 

0.0000 

e  =  0.92 

Critical  Flux  =  14.0  kW/m2 
1^=  16.6  W/m2K 
hj„=  39.1  W/m2K 
TiB  =  349  °C 

kpc  =  1.4393  kW2-s/m4-K2 
Minimum  Flux  =  18.0  kW/m2 


MATERIAL  NO.  9 


0 


10 


Heat  Flux  (kW/m2) 
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Modified  Heat  Release  Rate  Properties 

HEAT  RELEASE  RATE  -  50  kW/m2 


MATERIAL#! 


Time  (sec) 


01 

+ r 

tb 

HRR3omax 

THR 

X 

|  324  s 

241  s 

30  kW/m2 

4.40  MJ/m2 

0.00630 

MATERIAL  #2 


60  120 

180 

240  300  360  420 

Time  (sec) 

480 

4 

HRR30max 

THR 

X 

NI 

NI 

_ 

1 

1 

1 

0.00  MJ/m2 

0 

540 


600 
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Heat  Release  Rate  (kW/m!)  Heat  Release  Rate  (kW/m1) 


Modified  Heat  Release  Rate  Properties 


HEAT  RELEASE  RATE  -  50  kW/m2 

MATERIAL  #3 


Time  (sec) 


tig 

tb 

HRRjOma* 

THR 

X 

65  s 

720  s 

109  kW/m2 

44.72  MJ/m2 

0.00194  1/s 

MATERIAL  #4 


Time  (sec) 


HQ 

tb 

HRR30max 

THR 

X 

75  s 

978  s 

124  kW/m2 

66.44  MJ/m2 

0.0015  1/s 
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Heat  Release  Rate  (kVV/m2)  Heat  Release  Rate  (kW/m2) 


Heat  Release  Rate  (kW/m1)  Heat  Release  Rate  (kW/m2) 


Modified  Heat  Release  Rate  Properties 

HEAT  RELEASE  RATE  -  50  kW/mJ 

MATERIAL  #7 


tb 

HRR30max 

THR 

X 

28  s 

143  s 

39  jkW/m2 

5.58  MJ/m2 

IE-15  1/s 

MATERIAL  #8 


Time  (sec) 


tig 

tb 

HRR3omax 

THR 

30  s 

355  s 

350  kW/m2 

68.52  MJ/m2 

0.00474  1/s 

'-5 


Heat  Release  Rate  (kW/m1) 


Modified  Heat  Release  Rate  Properties 


HEAT  RELEASE  RATE  -  50  kW/m2 


MATERIAL  #9 


Time  (sec) 


*ig 

tb 

HRR30mm 

THR 

93  s 

1008  s 

118  kW/m2 

47.59  MJ/m2 

0.00243  1/s 

'-6 


Heat  Release  Rate  (kVV/m!)  Heat  Release  Rate  (kW/m2) 


Modified  Heat  Release  Rate  Properties 


HEAT  RELEASE  RATE  -  75  kW/m2 


MATERIAL  #2 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/m1)  Smoke  Production  Rate  (m2/s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 


Time  (s) 

HEAT  FLUX  TO  THE  FLOOR:  MATERIAL#! 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/m1)  Smoke  Production  Rate  (mVs) 


5 


SMOKE  PRODUCTION  RATE:  MATERIAL  #2 


Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/mJ)  Smoke  Production  Rate  (m2/s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


1200 


HEAT  RELEASE  RATE:  MATERIAL  #4 


Time  (s) 

UPPER  GAS  LAYER  TEMPERATURE:  MATERIAL  #4 


Time  (s) 


Heat  Flux  (kW/m2)  Smoke  Production  Rate  (m2/s) 


SMOKE  PRODUCTION  RATE:  MATERIAL  #4 


0  120  240  360  480  600  720  840  960  1080  1200  1320 


Time  (s) 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/m2)  Smoke  Production  Rate  (mVs) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/m5)  ,  Smoke  Production  Rate  (mJ/s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 


HEAT  FLUX  TO  THE  FLOOR:  MATERIAL  #6 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  (kW/tn1)  Smoke  Production  Rate  (mVs) 


Temperature  (°C)  Heat  Release  Rate  (kW) 


1200 


HEAT  RELEASE  RATE:  MATERIAL  #8 


Time  (s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 
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Time  (s) 


Heat  Flux  (kW/mJ)  Smoke  Production  Rate  (mJ/s) 


SMOKE  PRODUCTION  RATE:  MATERIAL  #8 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 


0  120  240  360  480  600  720  840  960  1080  1200  1320 

Time  (s) 
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Temperature  (°C)  Heat  Release  Rate  (kW) 


Heat  Flux  {kW/m1)  Smoke  Production  Rate  (m2/s) 


SMOKE  PRODUCTION  RATE:  MATERIAL  #9 


Time  (s) 
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Symbols 


NOMENCLATURE 


CP 

J/kg-K  or  kJ/kg-K 

Specific  heat  at  constant  pressure 

D 

m 

Characteristic  dimension  of  ignition  burner 

AHch 

kJ/g 

Chemical  heat  of  combustion 

AHeff 

kJ/g 

Effective  heat  of  combustion 

AH  net, c 

kJ/g 

Net  heat  of  complete  combustion 

AH  rad 

kJ/g 

Radiative  heat  of  combustion 

k 

W/m-K  or  kW/m-K 

Material  thermal  conductivity 

fcrit 

W/m2  or  kW/m2 

Critical  heat  flux  for  ignition 

Qo 

t 

Non-dimensional  heat  release  rate 

s 

m 

Mean  beam  length 

Tts 

K 

Material  ignition  temperature 

Ts.min 

K 

Minimum  surface  temperature  for  lateral  flame  spread 

m 

Flame  height  (used  in  flame  height  correlation) 

Greek 

Symbols 

U flame 

t 

Flame  Absorptivity 

Xa 

t 

Combustion  efficiency 

Xr 

t 

Radiative  fraction  of  heat  release  rate 

8 

m 

Material  thickness 

Aflame 

t 

Emissivity  of  flames  from  the  burning  material 

% 

t 

Emissivity  of  the  upper  gas  layer 

£k 

t 

Material  surface  emissivity  of  element  k 

t 

Total  emissivity  of  upper  gas  layer 

0 

kW2/m3 

Lateral  flame  spread  parameter 

K 

1/m 

Mean  absorption  coefficient  (0.47) 

Ks 

1/m 

Absorption  coefficient  of  the  soot  in  the  upper  layer 

P 

kg/m3 

Material  density 

* flame 

t 

Flame  transmissivity 

Notes: 

+  Dimensionless  variable 
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1.0  INTRODUCTION 


This  report  contains  an  evaluation  of  an  enhanced  version  of  the  flamespread  algorithm  of  Mitler 
(Mitler,  1993,  1994)  as  implemented  in  CFAST  (Peacock,  1993,  1997)  by  WPI FPE  MS  student  Mark 
Wright.  It  is  based  on  his  thesis  (Wright,  1999)  which  discusses  the  material  contained  in  this  report  in 
greater  detail.  The  changes  that  were  made  to  the  original  algorithm  of  Mitler  (Mitler,  1993,  1994)  are 
discussed  and  the  processes  through  which  material  input  properties  are  derived  from  bench-scale 
experimental  data  are  explained.  The  flame  spread  algorithm  has  been  evaluated  against  the  data 
provided  by  Janssens  (Janssens,  1998).  A  sensitivity  analysis  was  conducted  to  better  understand  the 
dynamics  of  the  model  with  the  intention  of  identifying  areas  that  may  require  higher  precision.  This 
sensitivity  analysis  also  focuses  on  quantifying  the  effects  of  cone  calorimeter  (ISO  5660),  and  LIFT 
apparatus  (IMO  Resolution  A.653,  ASTM  E  1321-97-a)  data  interpretation  methods  on  the  predictions  of 
the  model.  Additionally  the  model’s  strengths  and  weakness  are  identified,  and  suggestions  are  made  as 
to  how  research  should  proceed  to  improve  upon  the  model. 

2.0  MODIFICATIONS  TO  ORIGINAL  FLAME  SPREAD  ALGORITHM 

The  original  flame  spread  algorithm  (Mitler,  1993,  1994)  was  designed  to  simulate  flame  spread 
along  a  side-wall  with  the  burner  located  at  a  distance  from  the  comer  such  that  the  comer  effects  were 
negligible.  The  algorithm  has  been  modified  to  simulate  the  ISO  9705  test  configuration.  Several  aspects 
of  the  algorithm  required  modification  to  extend  the  model  to  the  comer  configuration.  These 
modifications  are  discussed  in  the  following  sections. 

2.1  Change  in  Geometry  for  Corner  Configuration 

The  change  in  geometry  of  flame  spread  on  the  walls  themselves  requires  little  modification  when 
altering  the  algorithm  for  the  comer  configuration.  The  area  of  spread  is  “folded”  into  the  comer  along  its 
centerline.  The  change  in  geometry  of  flame  spread  across  the  ceiling  is  more  difficult  to  define.  The 
two  comer  walls  constrain  the  flame  spread  across  the  ceiling,  creating  a  more  bounded  geometry.  It  is 
necessary  to  define  the  shape  of  flame  spread  across  the  ceiling  in  order  to  properly  track  the  pyrolysis 
and  burnout  front  positions  on  the  ceiling  and  to  calculate  the  enhanced  radiative  exchange  between  the 
walls  and  ceiling  in  the  comer. 

A  quarter-circular  geometry  was  used  to  represent  the  shape  of  the  pyrolysis  and  burnout  fronts 
on  the  ceiling.  This  was  done  primarily  because  the  geometry  of  spread  is  intuitively  realistic  and  the 
radiation  calculations  are  slightly  easier  with  the  quarter-circular  geometry  of  spread  as  opposed  to  a 
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triangular  or  square  geometry.  With  this  geometry,  the  pyrolysis  and  burnout  fronts  spread  radially  from 
the  comer,  with  growth  of  the  pyrolysis  area  simplified  as  a  radial  progression. 

2.2  Corner  Flame  Height  Correlation 

The  height  of  flames  from  a  burner  in  a  compartment  comer  is  different  from  those  against  the 
center  of  a  compartment  wall  (Drysdale,  1999).  Therefore,  it  was  necessary  to  develop  an  alternate 
expression  for  flame  height  in  the  comer.  Several  correlations  for  square  burners  in  a  comer  were 

examined  using  the  standard  non-dimensional  group  Q*D  (Drysdale,  1999).  The  first  was  developed  by 
Hasemi  and  Tokunaga  (Hasemi,  1984)  for  square  burners  in  an  open  comer  configuration.  The  next 
correlation  (identified  by  Dillon  (Dillon,  1998))  was  developed  by  Kokkala.  The  last  correlation  for  a 
square  burner  flame  height  in  a  compartment  comer  (also  identified  by  Dillon  (Dillon,  1998))  was 
originally  developed  by  Heskestad  and  revised  by  Kokkala. 

Values  of  Zf/D  versus  QD  were  calculated  for  each  of  the  three  individual  correlations.  These 
three  curves  were  then  collapsed  into  one  by  taking  the  arithmetic  mean  of  all  three  correlations  for  Zf/D 
versus  QD  in  order  to  arrive  at  an  approximation  of  flame  height  from  a  square  burner  in  a  comer.  A 
power  law  curve-fit  was  backed  out  of  this  arithmetic  mean  and  the  resulting  expression,  as  given  by 
Equation  (B2-1),  replaces  Mitler’s  original  flame  height  correlation. 

-^  =  -1.7  +  5.8(g/))  (B2-1) 

2.3  Enhanced  Radiative  Exchange  in  the  Corner 

The  incident  radiation  distribution  to  the  wall  and  ceiling  in  the  comer  configuration  is  different 
as  compared  to  the  side-wall  configuration  because  the  geometry  allows  for  radiative  exchange  between 
the  two  comer  walls  and  also  between  the  comer  walls  and  the  ceiling.  The  total  incident  flux  is 
enhanced,  and  this  flux  in  turn  controls  the  mass  loss  rate  for  a  pyrolyzing  region  or  adds  to  the  degree  of 
preheating  for  a  non-pyrolyzing  region. 

A  subroutine  was  added  that  calculates  the  radiative  exchange  among  the  walls,  ceiling,  and  the 
hot  upper  gas  layer  in  the  comer.  The  radiative  fluxes  from  the  burner  and  wall  flames  are  also  added  as 
external  fluxes.  A  factor  to  attenuate  radiation  as  it  passes  through  a  flame  is  included.  It  is  assumed  that 
both  the  walls  and  the  upper  layer  gas  are  gray  bodies  with  uniform  properties.  The  emissivity  of  the  gas 
layer  without  soot  (sg)  is  calculated  in  the  usual  way,  using  the  partial  pressures  of  H20  and  C02,  beam 
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length,  and  the  Hottel  (Tien,  1995)  charts.  A  subroutine  was  borrowed  directly  from  the  source  code  of 
the  WPI/Fire  model  (Satterfield,  1990)  for  this  purpose.  The  equation  used  to  calculate  the  emissivity  of 
the  hot  upper  gas  layer  is  from  Tien,  et  al.  (Tien,  1995)  and  is  given  as  Equation  (B2-2): 

£,  =(l- e~KsS  )  +  £ge~KsS  (B2-2) 

The  walls  are  discretized  into  strips  with  a  width  of  approximately  0.12m  because  this  element 
size  was  found  to  provide  adequate  resolution  while  allowing  acceptable  execution  times.  The  geometry 
of  the  comer  was  broken  into  five  separate  zones  and  the  upper  gas  layer,  as  shown  in  Figure  B2-1 .  View 
factor  algebra  and  reciprocity  were  used  to  generate  a  radiation  network.  The  resulting  matrix  is  solved 
by  LU  factorization  and  the  overall  radiative  exchange  is  calculated  between  wall  elements  on  opposing 
sides  of  the  comer,  as  well  as  between  wall  and  ceiling  elements.  The  main  output  of  this  subroutine  is 
the  incident  radiant  flux  at  each  element. 


Figure  B2-1.  -  Geometry  of  the  elements  used  in  radiation  network. 
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2.4  Experimental  Corner  Heat  Flux  Map 

The  original  flame  spread  algorithm  (Mitler,  1993, 1994)  simulated  the  heat  flux  from  a  line 
burner  to  a  sidewall  by  adding  separate  theoretical  expressions  for  convective  and  radiative  heat  transfer 
to  calculate  a  total  flux.  A  different  approach  was  used  in  the  modified  algorithm.  Experimental 
information  on  the  flux  distribution  in  the  ISO  9705  configuration  was  recorded  by  Dillon  (Dillon,  1998). 
His  graphs  provide  two-dimensional  distributions  for  the  total  (i.e.,  convective  and  radiative)  flux; 
however,  the  model  only  requires  a  one-dimensional  distribution  of  the  flux.  For  this  reason,  the  flux 
distribution  maps  were  translated  to  a  one-dimensional  analogue. 

The  original  flux  distributions  caused  the  model  to  behave  somewhat  unrealistically.  Ignition 
always  occurred  on  the  wall  just  below  the  ceiling  at  the  location  of  highest  flux,  but  general  observation 
on  comer  fire  experiments  indicates  that  sustained  ignition  occurs  farther  down  the  wall.  As  a  result  of 
this  discrepancy,  the  flux  distributions  were  altered  to  ensure  ignition  at  the  area  of  elevated  flux  on  the 
wall  behind  the  burner  flame.  For  the  100  kW  burner  distribution,  the  flux  was  decreased  to  30  kW/m2  at 
a  height  of  approximately  1.3-m  and  set  constant  until  the  30-kW/m2  boundary  was  reached  on  the 
ceiling.  A  similar  procedure  was  followed  for  the  300  kW  burner  distribution,  except  the  cutoff  was 
60kW/m2  at  a  height  of  approximately  1 .7-m. 

2.5  Multiple  Pyrolysis  Zones 

Provisions  were  made  to  track  multiple  pyrolysis  zones  in  the  modified  flame  spread  algorithm. 
The  original  algorithm  allowed  only  one  pyrolysis  zone  at  any  given  time  in  the  simulation.  All  elements 
below  a  node  at  which  bum-out  occurred  were  marked  as  non-burning.  This  often  led  to  the  situation 
where  a  significant  number  of  elements  were  extinguished  even  though  there  was  still  available  fuel.  The 
shift  in  the  lower  pyrolysis  bound  at  bum-out  was  removed  to  correct  this  problem.  Bum-out  is  now 
controlled  by  the  remaining  fuel  at  each  node.  Once  multiple  pyrolysis  zones  exist  and  they  are 
sufficiently  far  apart,  the  validity  of  the  flame  height  and  heat  flux  calculations  is  questionable  because 
the  algorithm  was  written  to  simulate  only  one  pyrolysis  zone. 

2.6  Insertion  Into  CFAST 

The  modified  flame  spread  algorithm  was  incorporated  into  the  zone  model  CFAST  (Peacock, 
1993,  1997).  Gas  temperatures,  wall  temperatures,  interface  height,  etc.  are  calculated  internally  by  the 
zone  model.  CFAST  has  the  ability  to  simulate  the  burning  of  multiple  objects  within  a  compartment. 
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The  time  at  which  ignition  occurs  can  be  specified,  or  CFAST  can  dynamically  determine  when  ignition 
occurs  based  on  the  localized  attainment  of  an  ignition  criterion,  such  as  Tjg  or  qcrit .  The  burning  of  an 

object  is  governed  by  its  thermophysical  properties  as  specified  in  CFAST’ s  integrated  thermal  and  object 
databases.  CFAST  tracks  the  combustion  of  a  burning  object  and  calculates  its  effect  on  the  overall 
compartmental  conditions.  The  new  flame  spread  algorithm  was  inserted  into  CFAST  (Peacock,  1993, 
1997)  as  a  new  object  type,  designated  as  object  type  five. 

3.0  DEVELOPMENT  OF  INPUT  PARAMETERS 

The  data  used  to  derive  the  input  parameters  required  by  the  flame  spread  algorithm  were 
provided  by  Southwest  Research  Institute  (San  Antonio,  TX).  These  data  were  the  basis  of  a  report 
(Janssens,  1998)  that  described  the  testing  of  seven  fire-retardant  (FR)  woven  glass  composite  materials,  a 
textile  wall  covering  material,  and  one  non-fire-retardant  composite  material.  Each  material  was  tested  in 
the  cone  calorimeter  (ISO  5660),  LIFT  apparatus  (IMO  Resolution  A.653,  ASTM  E  1321-97-a),  and  full- 
scale  room  scenario  (ISO  9705).  The  material  ignition  and  lateral  flame  spread  properties  were  derived 
using  the  procedure  described  in  ASTM  E  1321.  Out  of  the  nine  materials  tested,  four  produced  complete 
sets  of  material  data.  The  other  five  were  missing  values  for  the  minimum  temperature  for  lateral  flame 
spread  (T5, „„,,),  and  the  lateral  flame  spread  parameter  («P).  This  section  describes  the  input  parameters 
required  by  the  flamespread  algorithm  and  explains  how  they  are  derived  from  experimental  data. 

3.1  Specific  Heat  and  Thermal  Conductivity 

The  available  experimental  data  provided  a  value  for  each  material’s  thermal  inertia  ( kpcp ); 
however,  the  input  routine  of  the  flame  spread  algorithm  requires  the  individual  specification  of  thermal 
conductivity  (k),  density  (p),  and  specific  heat  (cp).  The  experimental  data  provided  a  separate  value  for 
material  density  ( p ).  A  study  by  Grenier  (Grenier,  1998)  on  similar  composite  materials  was  consulted  to 
obtain  a  value  for  specific  heat  (cp)  of  1000  J/kg-K.  The  value  for  thermal  conductivity  was  backed  out 
from  kpcp,  using  the  known  values  for  p  and  cp. 

3.2  Net  Heat  of  Complete  Combustion  and  Combustion  Efficiency 

The  algorithm  is  set  up  to  use  an  effective  heat  of  combustion  by  multiplying  a  net  heat  of 
complete  combustion  by  a  combustion  efficiency  as  follows: 

AHeff={XA){^nel,c)  (B3-1) 
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The  primaiy  material  that  is  consumed  during  combustion  of  composite  materials  is  the  resin  that  binds 
the  glass  together.  For  this  reason,  the  net  heat  of  complete  combustion  for  the  resin  material  alone  is 
used  with  an  appropriate  combustion  efficiency  to  arrive  at  the  effective  heat  of  combustion  derived  from 
the  cone  calorimeter  tests.  The  net  heat  of  complete  combustion  for  these  materials  was  taken  from 
Tewarson  (Tewarson,  1995).  If  no  data  could  be  found  on  a  particular  material,  a  combustion  efficiency 
of  0.5  was  assumed  and  an  appropriate  net  heat  of  complete  combustion  was  used  to  arrive  at  the 
experimentally  determined  effective  heat  of  combustion. 

3.3  Radiative  Fraction  of  Heat  Release  Rate 

The  radiative  fraction  of  heat  release  rate  (#?)  was  derived  in  a  similar  manner  as  the  combustion 
efficiency,  using  the  radiative  heat  of  combustion  and  the  chemical  heat  of  combustion: 


Xr  = 


A#™, 


(B3-2) 


The  chemical  heat  of  combustion  of  the  resin  material  alone  was  used  in  this  case  instead  of  the  value  for 
the  overall  composite  material  provided  by  Janssens  (Janssens,  1998).  If  these  values  could  not  be  found, 
a  radiative  fraction  of  0.3  was  assumed  because  this  is  a  nominal  value  for  most  fuels  (Drysdale,  1999). 


3.4  Emissivity 


The  emissivity  of  these  materials  was  not  determined  during  testing.  Based  on  Grenier’s  study  of 
similar  materials  (Grenier,  1998),  an  emissivity  of  0.9  was  assumed  to  be  reasonable  for  these 
simulations. 

3.5  Flame  Transmissivity 

A  value  for  flame  transmissivity  was  also  not  determined  during  testing  of  the  composite 
materials,  and  little  guidance  on  this  value  was  found  in  the  literature.  The  only  basis  found  was  from 
Karlsson  (Karlsson,  1992),  where  he  assumed  a  flame  emissivity  of  0.5.  If  the  assumption  is  made  that 
Efhme  =  0Cfiame,  and  af\ame  =  (1  -  Tfiame),  then  Tflame  is  0.5.  This  is  the  value  currently  used  for  these  materials. 
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3.6  Lateral  Flame  Spread  Properties 

Two  sets  of  the  values  governing  lateral  flame  spread  (Ts  min  and  0)  were  available  for  each 
material.  This  is  due  to  the  fact  that  data  on  Ts  min  and  0  were  reduced  following  two  different 
experimental  procedures,  specifically  IMO  Resolution  A.653,  and  ASTM  E  1321-97-a.  The  average 
values  from  the  two  data  reduction  methods  were  used  in  the  simulations. 

3.7  Cone  Calorimeter  Curve  at  Multiple  Irradiances 

Several  cone  calorimeter  tests  were  performed  on  each  material.  Generally,  two  tests  were 
completed  at  each  of  three  irradiance  levels.  The  ignition  times  of  the  two  curves  were  aligned  and  the 
arithmetic  mean  of  heat  release  per  unit  area  versus  time  was  found  at  each  discrete  time  for  which  data 
were  available.  This  created  an  average  heat  release  rate  curve  for  each  material  at  each  irradiance  level. 
In  some  cases,  the  bum-time  of  one  cone  calorimeter  test  was  significantly  longer  than  that  of  the  other, 
leading  to  a  possible  spike  (i.e.,  a  curve  that  is  discontinuous  in  the  first  derivative)  when  the  shorter 
curve  ends.  In  order  to  resolve  this  issue,  the  newly  created  average  curve  was  extended  via  extrapolation 
from  the  time  corresponding  to  the  end  of  the  shorter  cone  curve.  More  specifically,  the  tail  end  of  the 
longer  cone  curve  was  shifted  vertically  upward  or  vertically  downward  to  avoid  discontinuity  at  the  point 
where  the  shorter  cone  curve  ends,  as  shown  in  Figure  B3-1  at  730  seconds. 

Once  the  corrected  average  curve  had  been  developed,  a  number  of  points  (typically  on  the  order 
of  five  to  twenty  points,  depending  on  the  complexity  of  the  curve)  were  selected  that  captures  the  overall 
behavior  of  the  cone  calorimeter  curve.  A  simple  five-point  moving  average  was  used  to  help  select 
points  when  the  corrected  average  curve  exhibited  a  significant  degree  of  noise.  These  are  the  points  that 
are  provided  to  the  flamespread  algorithm  as  input.  An  example  is  shown  in  Figure  B3-1 .  These  are 
provided  in  Section  7  for  all  nine  materials. 
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Heat  Release  Rate  In  Cone  Calorimeter 
Material  f4  ©75  kW/m2  Irradiance 


4.0  SIMULATION  RESULTS  AND  MODEL  VALIDATION 

This  section  presents  selected  results  generated  by  the  modified  flame  spread  algorithm  and 
compares  them  to  experimental  data.  Heat  release  rate,  upper  layer  gas  temperature,  smoke  production 
rate,  and  heat  flux  to  the  floor  of  the  compartment  were  selected  as  validation  criteria  and  are  discussed  in 
this  section.  Although  graphical  data  from  individual  runs  are  presented,  the  results  are  discussed  from  a 
holistic  point  of  view  because  the  qualitative  relationships  among  the  four  validation  criteria  are 
consistent  for  all  materials.  The  baseline  results  for  all  nine  materials  are  given  in  Sections  8  through  1 1 . 
A  summary  of  the  simulation  results  compared  to  experimental  data  is  given  in  Table  B4-1  and  the  times 
to  flashover  are  given  in  Table  B4-2.  A  total  heat  release  rate  of  1000  kW  was  used  as  the  criterion  for 
flashover  except  for  material  #5  where  750  kW  was  used  because  flames  were  seen  to  come  from  the  door 
during  full-scale  testing  at  a  heat  release  rate  of  approximately  750  kW.  In  Table  B4-1,  “Net  Avg.”  refers 
to  the  average  net  HRR  (i.e.,  burner  HRR  has  been  subtracted)  over  the  entire  test,  “Net  30s  Avg  Max”  is 
the  maximum  net  30-second  (or  60-second)  sliding  average,  and  “Net  Max”  is  the  highest  instantaneous 
net  heat  release  rate  that  was  predicted  during  the  simulation. 
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Table  B4-1.  -  Comparison  of  predicted  HRR  and  SPR  to  experimental  data. 


Heat  Release  Rate  (kW)  j 

|  Smoke  Production  Rate  (m2/s) 

Experimental 

Simulation  \ 

I  Experimental 

Simulation 

Material 

Net 

Net 

Net 

Net 

Net 

Net 

Net 

Net 

Net 

Net 

Avg 

30s  Avg 
Max 

Avg 

30s  Avg 
Max 

Max 

Avg 

60s  Avg 
Max 

Avg 

60s  Avg 
Max 

Max 

i 

62 

62 

48 

192 

194 

1.50 

5.41 

2.51 

7.45 

7.71 

2 

31 

31 

28 

109 

118 

0.15 

0.47 

1.77 

4.33 

4.66 

3 

191 

191 

119 

366 

388 

10.00 

21.70 

8!66 

23.34 

25.05 

4 

190 

190 

152 

660 

693 

9.08 

32.10 

10.17 

38.23 

41.20 

5 

115 

115 

44 

218 

223 

6.39 

26.40 

5.80 

22.69 

24.97 

6 

28 

28 

67 

164 

166 

1.45 

3.46 

8.27 

18.91 

19.25 

7 

17 

17 

11 

77 

82 

0.10 

0.16 

0.91 

1.41 

1.41 

8 

170 

170 

131 

361 

901 

2.28 

4.10 

3.86 

5.30 

32.23 

9 

109 

109 

109 

488 

512 

0.42 

3.81 

1.38 

3.30 

3.51 

Table  B4-2.  -  Time  to  flashover. 


Material 

Experimental 

Model 

3 

5.8  minutes 

N/A 

4 

5.1  minutes 

N/A 

5 

16.4  minutes* 

N/A 

8 

1.8  minutes 

1.4  minutes 

9 

1 1 .3  minutes 

N/A 

*  Flashover  assumed  to  be  total  HRR  of  750kW 


4.1  Interpretation  of  the  Pyrolysis/Burnout  Front  Graph 

The  pyrolysis/bumout  front  graphs  presented  in  the  following  discussion  are  somewhat  difficult 
to  understand  at  first  glance;  therefore,  an  example  is  provided  as  Figure  B4-1  for  clarity.  The  arrows 
indicate  the  general  direction  of  the  fronts,  and  the  hatched  area  represents  burning.  The  thick  dashed  line 
represents  the  boundary  between  the  wall  and  the  ceiling,  with  all  points  above  this  line  being  on  the 
ceiling.  The  height  scale  does  not  change  for  the  ceiling,  but  the  length  scale  on  the  ceiling  is  a  radius 
instead  of  a  vertical  height.  To  obtain  a  radius,  the  virtual  height  of  the  room  (ceiling  height  minus  the 
height  of  the  burner)  is  subtracted  from  the  value  on  the  graph.  This  graph  can  also  be  used  to  determine 
the  burning  time  of  a  particular  element  by  choosing  a  height  and  moving  horizontally  to  estimate  the 
time  that  elapses  between  the  passage  of  the  pyrolysis  and  burnout  fronts. 
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4.2  Heat  Release  Rate 

The  simulated  heat  release  rate  of  the  marine  grade  polyester  composite  (material  #8)  follows  the 
experimental  curve  quite  well,  as  shown  by  its  baseline  curve  in  Figure  B4-2.  The  fire  growth  on  this 
material  is  quite  rapid,  and  flashover  conditions  (heat  release  rate  of  1MW)  are  reached  roughly  40 
seconds  after  ignition.  The  entire  wall  is  burning  and  the  ceiling  ignites  about  15  seconds  after  ignition. 
The  rapid  growth  after  ceiling  ignition  is  a  function  of  the  burning  nodes  quickly  reaching  the  peak 
burning  rate  on  their  cone  calorimeter  curve  combined  with  rapid  growth  across  the  ceiling  elements. 
Unlike  the  wall  elements,  the  ceiling  elements  increase  in  area  the  further  they  are  from  the  comer.  This 
characteristic  of  the  algorithm  is  important  to  remember  when  examining  spread  across  the  ceiling  and  the 
resulting  rate  of  heat  release. 

The  peak  heat  release  rate  from  the  cone  calorimeter  is  quite  high  for  this  material  (~350kW/m2), 
and  is  probably  the  most  dominant  factor  in  its  rapid  fire  growth.  As  will  be  discussed  in  Section  5,  some 
slight  perturbations  from  the  baseline  curve  occur  when  varying  the  parameters  in  the  sensitivity  matrix, 
but  the  overall  curvature  of  the  simulation  for  material  #8  remains  similar  to  the  baseline  graph  shown  in 
Figure  B4-2. 


B-16 


Heat  Release  Rate  vs.  Time 
Material  f8  -  Polyester 
Baseline  va.  Experimental 


Time  (seconds) 

Figure  B4-2.  -  Baseline  heat  release  rate  for  material  #8. 


The  behavior  of  the  fire-retardant  materials  is  significantly  different  from  the  marine  grade 
polyester  composite.  A  total  of  seven  fire-retardant  materials  were  tested,  and  these  materials  can  be 
divided  into  two  separate  sub-categories.  The  first  category  consists  of  the  materials  that  exhibited  little 
burning  during  full  scale  testing.  The  materials  that  did  now  show  significant  burning  are  material  #1  (FR 
phenolic),  material  #2  (fire  restricting  material),  and  material  #6  (coated  FR  epoxy).  The  second  sub¬ 
category  consists  of  the  fire-retardant  materials  that  caused  flashover  during  full-scale  testing.  These  are 
material  #3  (FR  polyester),  material  #4  (FR  vinylester),  material  #5  (FR  epoxy),  and  material  #9  (FR 
acrylic).  It  should  be  noted  that  1000  kW  was  used  as  the  criterion  for  flashover  except  with  material  #5. 
This  was  necessary  because  flames  were  seen  to  come  from  the  door  during  full-scale  testing  of  material 
#5  at  a  time  corresponding  to  a  heat  release  rate  of  approximately  750  kW;  therefore,  this  value  is  the 
criterion  for  flashover  for  material  #5.  Material  #7  was  predicted  but  its  results  are  not  considered  to  be 
valid  because  the  material  was  seen  to  fall  on  the  floor  during  full-scale  testing. 

The  model  captured  the  general  behavior  of  the  materials  that  did  not  exhibit  significant  burning. 
Figure  B4-3  shows  a  comparison  between  the  experimental  data  and  the  model  prediction  for  material  #1 . 
This  behavior  is  qualitatively  comparable  to  the  remaining  cases  that  did  not  show  significant  burning. 
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Heat  Release  Rate  ve.  Time 
Material  fl  -  FR  Phenolic 
Baseline  ve.  Experimental 


Figure  B4-3.  -  Baseline  heat  release  rate  for  material  #1. 


The  simulated  behavior  of  the  fire-retardant  composite  materials  that  went  to  flashover  during 
full-scale  testing  is  significantly  different  from  the  previously  discussed  cases.  The  remainder  of  this 
section  examines  the  behavior  of  these  materials  in  detail.  The  simulations  for  materials  #3  and  #4  do  not 
match  the  experimental  curves  well  except  near  ignition,  as  can  be  seen  from  Figures  4-4  and  4-5.  The 
total  experimental  heat  release  rate  of  material  #5  was  seen  to  approach  800  kW  after  the  ignition  burner 
strength  was  increased  to  300  kW.  The  model,  as  shown  in  Figure  B4-6,  did  not  capture  this  behavior 
after  600  seconds.  The  predicted  behavior  of  material  #9  follows  the  experimental  results  fairly  well 
before  600  seconds  but  the  model  does  not  track  the  experimental  behavior  after  the  burner  strength 
increase,  as  shown  in  Figure  B4-7. 


Comparison  of  Experimental  HRRto  Model  Results 
Material  (6  -  FR  Epoxy 
Bacellne  vt.  Experimental 


Figure  B4-6.  -  Baseline  heat  release  rate  for  material  #5. 


Heat  Release  Rate  vs.  Time 
Material  *9  ■  FR  Acrylic 
Baseline  vs.  Experimental 


Figure  B4-7.-  Baseline  heat  release  rate  for  material  #9. 


The  results  from  material  #4  shown  in  Figure  B4-5  are  used  to  explain  the  simulated  behavior 
generally  exhibited  by  these  fire-retardant  materials.  After  ignition,  a  period  of  sustained  growth  occurs 
as  flame  spreads  over  the  area  of  high  heat  flux  on  the  wall,  followed  shortly  by  the  spread  of  flame  over 
the  area  of  constant  heat  flux  specified  by  the  heat  flux  map  (115  seconds  to  180  seconds).  Figure  B4-8 
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shows  the  heat  flux  distribution  at  three  discrete  times  for  material  #4.  The  area  of  high  flux  is  followed 
by  an  area  of  constant  flux  that  ends  at  a  height  of  approximately  2.6-m,  corresponding  to  a  radius  of  0.35 
m  on  the  ceiling.  Figure  B4-9  shows  the  pyrolysis  and  burnout  front  locations  versus  time  for  material 
#4.  The  initial  rapid  growth  occurs  from  1 15  seconds  to  180  seconds,  ending  at  approximately  2.6  m. 
Thus,  the  initial  growth  is  directly  related  to  the  initial  heat  flux  distribution. 

The  section  of  the  curve  from  180  to  600  seconds  in  Figure  B4-5  is  representative  of  slow  flame 
spread  across  the  ceiling  coupled  with  nodes  beginning  to  bum  with  a  heat  release  rate  determined  by  the 
decay  portion  of  the  cone  calorimeter  curve.  This  results  in  a  slower  rate  of  growth.  The  burnout 
simulated  by  the  algorithm  also  tends  to  occur  more  rapidly  on  the  ceiling  as  opposed  to  the  wall,  i.e.,  the 
ceiling  releases  a  greater  amount  of  energy  in  a  shorter  period  of  time. 


Total  Heat  Flux  va.  Distance 
Material  #4  -  FR  VInylester 
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Pyrolycis  and  Burnout  Fronts  v«.  Time 
Material  <4  -  FR  Vlnyleeter 


The  section  of  the  heat  release  rate  curve  after  600  seconds  in  Figure  B4-5  shows  the  behavior  of 
the  material  after  the  ignition  burner  strength  is  increased  from  100  kW  to  300  kW.  It  is  typical  for  the 
remaining  non-burning  nodes  to  ignite  in  a  short  span  of  time,  as  occurs  from  600  to  675  seconds  in  the 
sample  figure.  If  initial  burnout  does  not  happen  before  600  seconds,  it  generally  occurs  during  this 
segment  of  the  curve.  Note  that  the  slope  of  heat  release  rate  curve  in  this  segment  is  also  quite  high. 

This  is  similar  to  the  high  slope  shown  for  material  #8  in  Figure  B4-2  after  the  entire  wall  has  ignited.  In 
Figure  B4-5,  the  rapid  growth  after  the  burner  output  increase  is  a  direct  result  of  the  increasing  size  of 
the  elements  across  the  ceiling  as  previously  described.  The  segment  of  the  curve  in  Figure  B4-5  from 
675  to  1200  seconds  is  indicative  of  additional  burnout  and  continued  burning  of  nodes  with  their  heat 
release  rates  being  determined  by  the  decay  portion  of  the  cone  calorimeter  curve. 

4.3  Upper  Layer  Gas  Temperature 

The  upper  and  lower  layer  gas  temperatures  (as  well  as  wall  temperatures)  are  calculated 
internally  by  CFAST.  The  experimental  data  provided  by  Janssens  (Janssens,  1998)  include  temperatures 
from  several  thermocouples  positioned  throughout  the  test  room.  It  should  be  noted  that  these 
thermocouple  readings  might  not  represent  a  true  gas  temperature  due  to  radiative  losses  or  gains, 
depending  on  their  location  in  the  compartment.  CFAST  calculates  the  mean  temperature  of  a  well- 
mixed,  homogeneous  gas  layer,  whereas  these  thermocouples  record  temperature  at  several  discrete 
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locations  in  the  test  room.  These  thermocouples  can  be  used  to  estimate  the  gas  temperatures  that  were 
present  during  full-scale  testing,  but  the  reading  of  a  single  thermocouple  is  not  directly  analogous  to  the 
temperature  of  a  well-mixed  homogeneous  gas  layer  as  is  reported  by  CFAST.  The  model  predictions 
were  compared  to  the  temperatures  of  a  thermocouple  at  the  center  of  the  room  and  a  thermocouple  in  the 
doorway. 

It  was  generally  found  that  when  the  heat  release  rate  is  over-predicted,  the  same  is  true  of  the 
upper  layer  gas  temperature.  A  comparison  of  experimental  to  predicted  data  for  material  #8  is  shown  in 
Figure  B4-10.  Although  the  upper  layer  gas  temperature  calculated  by  CFAST  is  higher  than  the 
thermocouple  readings,  this  is  consistent  with  its  heat  release  rate  curve  which  is  also  over-predicted  (see 
Figure  B4-2).  It  is  also  generally  true  that  an  under-predicted  heat  release  rate  corresponds  to  an  under¬ 
predicted  upper  layer  gas  temperature  compared  to  the  thermocouple  readings. 


Comparison  of  Thermocouple  Raadlng  and  Prsdlctad  Upper  Layer  Gae  Temperature 

Material  (8  -  Polyester 
Baseline  vs.  Experimental 


Tim*  {seconds} 

Figure  B4-10.  -  Upper  layer  gas  temperature  for  material  #8. 


It  was  found  that  the  upper  layer  gas  temperature  calculated  by  CFAST  was  similar  to  the 
thermocouple  temperatures  when  the  heat  release  rate  predicted  by  the  model  was  also  similar  to  that 
observed  experimentally.  The  results  for  materials  #1  and  #9  are  shown  in  Figures  4-11  and  4-12, 
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Temperature  (K)  Temperature  (K) 


respectively.  These  upper  layer  gas  temperatures  should  be  examined  concurrently  with  Figures  4-3  and 
4-7,  the  predicted  heat  release  rate  curves  for  materials  #1  and  #9. 


Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 

Material  *1  -  FR  Phenolic 
Baseline  vs.  Experimental 


Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 

Material  #9-FR  Acrylic 
Baseline  vs.  Experimental 
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4.4  Smoke  Production  Rate 


The  smoke  production  rate  was  calculated  from  the  upper  layer  optical  density  and  vent  flow  as 
computed  by  CFAST.  CFAST  allows  the  user  to  specify  the  yields  of  CO,  C02,  soot,  and  several  other 
chemical  species  for  each  individual  object  in  the  compartment.  The  yields  given  by  Tewarson 
(Tewarson,  1995)  for  the  resin  material  only  (as  opposed  to  a  glass-reinforced  resin)  were  used  as  input  to 
the  model.  No  yield  data  were  found  in  the  literature  for  vinylester  or  acrylic.  For  this  reason,  cone 
calorimeter  data  were  used  to  calculate  a  soot  yield  for  these  two  materials.  CFAST  requires  that  all 
species  yields  be  specified  as  a  ratio  to  the  yield  of  C02.  The  C02  yield  for  PMMA  was  substituted  for 
that  of  acrylic,  and  the  C02  yield  for  polyester  was  substituted  for  vinylester. 

The  main  CFAST  output  file  contains  a  time  history  of  species  concentrations  and  optical  density 
in  the  upper  and  lower  layers.  Also  given  are  the  temperature  of  each  layer  and  mass-based  vent  flows. 
The  mass  flow  out  of  the  upper  layer  was  converted  to  a  volumetric  flowrate  using  the  ideal  gas  law  and 
upper  layer  gas  temperature.  The  product  of  this  volumetric  flowrate  and  the  upper  layer  optical  density 
was  converted  to  a  smoke  production  rate  using  a  conversion  factor  of  2.3. 

As  was  seen  with  the  upper  layer  gas  temperatures,  the  smoke  production  rate  is  consistent  with 
the  heat  release  rate.  The  shape  of  the  smoke  production  rate  curve  follows  the  heat  release  rate  curve. 

An  over-prediction  in  the  heat  release  rate  generally  corresponds  to  an  over-prediction  in  the  smoke 
production  rate.  Conversely,  the  smoke  production  rate  is  generally  under-predicted  when  the  heat 
release  rate  is  under-predicted.  The  smoke  production  rates  for  materials  #1  and  #8  are  shown  in  Figures 
4-13  and  4-14.  Their  corresponding  heat  release  rate  curves  can  be  found  in  Figures  4-3  and  4-2. 
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Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  fl  -  FR  Phenolic 
Baeellne  ve.  Experimental 


Figure  B4-13.  -  Smoke  production  rate  material  #1. 


Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  *8  -  Polyester 
Baseline  vs.  Experimental 
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Figure  B4-14.  -  Smoke  production  rate  material  #8. 


4.5  Heat  Flux  to  the  Floor 

CFAST  allows  the  user  to  specify  targets  at  arbitrary  locations  and  orientations  throughout  the 
compartment.  The  flux  at  these  targets  is  then  calculated  internally  by  CFAST.  The  simulations  were  run 
with  a  target  placed  at  the  center  of  the  floor  to  estimate  this  heat  flux.  The  results  of  the  simulation  show 
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that  the  direct  radiation  from  the  ignition  burner  flame  and  wall  flames  are  not  included.  The  heat  flux  at 
this  target  is  primarily  due  to  radiative  heat  transfer  from  the  hot  upper  gas  layer.  For  this  reason,  the  heat 
flux  at  the  center  of  the  floor  was  consistently  underestimated  as  compared  to  experimental  data.  Future 
versions  of  the  model  will  include  a  module  to  account  for  direct  ignition  burner  and  wall  flame  radiation 
to  a  target  on  the  floor.  A  typical  time  history  of  heat  flux  to  the  floor  is  shown  in  Figure  B4-15.  It  can 
be  seen  that  the  simulation  results  are  consistent  with  the  upper  layer  gas  temperature.  Figure  B4-1 1 . 


Heat  Flux  to  the  Floor  v«.  Time 
Material  «1  •  FR  Phenolic 
Baaaline  va.  Experimental 


Figure  B4-15.  -  Heat  flux  to  the  floor  for  material  #1. 


5.0  SENSITIVITY  ANALYSIS 

A  self-consistent  set  of  input  data  was  assembled  for  each  material.  A  sensitivity  matrix  was 
developed  and  used  to  quantify  the  effects  of  altering  various  input  parameters.  The  key  baseline  input 
parameters  can  be  found  on  the  first  page  of  Section  7,  and  the  complete  set  of  input  parameters 
(including  algorithm  options)  for  the  nine  materials  immediately  follows.  With  the  intention  of  keeping 
this  report  brief,  the  sensitivity  analysis  results  from  four  of  the  nine  materials  will  be  discussed  here.  It 
should  be  noted  that  the  general  results  discussed  in  this  report  extend  to  all  materials  for  which 
simulations  were  run  because  similar  phenomena  were  observed  for  all  nine  materials.  The  complete  set 
of  baseline  results  can  be  found  in  Sections  8  through  11. 
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5.1  Selection  of  Parameters  for  Variation  in  Sensitivity  Analysis 

The  thermal  inertia  (kpcp)  of  the  lining  material  affects  ignition  time,  material  surface 
temperature,  and  fire  growth  rate;  therefore,  it  is  important  to  examine  the  effect  of  varying  kpcp.  A  value 
for  density  ip)  was  available,  and  the  value  for  specific  heat  ( cp )  was  held  constant  at  1000  J/kg-K  as 
previously  mentioned.  For  these  reasons  thermal  conductivity  ( k )  was  varied,  effectively  altering  the 
thermal  inertia.  A  factor  of  two  was  selected  because  work  conducted  by  Jacoby  (Jacoby,  1998)  on 
similar  composite  materials  showed  that  the  level  of  uncertainty  in  calculation  of  ignition  parameters  was 
roughly  within  a  factor  of  two. 

The  mass  loss  rate,  and  thus  the  heat  release  rate,  is  calculated  directly  from  the  cone  calorimeter 
curve  that  is  provided  to  the  model  as  input.  For  this  reason,  examining  the  effect  of  using  cone 
calorimeter  curves  at  multiple  irradiance  levels  is  important.  Average  cone  calorimeter  curves  from  three 
irradiance  levels  were  used  for  all  materials  apart  from  material  #2,  for  which  there  were  data  at  only  two 
irradiance  levels.  It  was  decided  that  the  cone  calorimeter  curves  at  the  median  irradiance  level  would  be 
considered  to  be  the  baseline  case,  although  this  is  somewhat  of  an  arbitrary  choice.  In  the  case  of 
material  #2,  the  lower  irradiance  (75  kW/m2)  was  chosm  to  constitute  the  baseline.  The  complete  set  of 
average  cone  calorimeter  curves  and  the  individual  points  that  were  given  to  the  model  as  input  can  be 
found  in  Section  7. 

The  material  emissivity  (e)  and  the  flame  transmissivity  itfiame)  were  selected  for  variation 
because  they  are  used  in  the  radiation  calculations.  These  variables  are  used  as  part  of  the  radiation 
network  used  to  compute  the  enhanced  radiative  exchange  in  the  comer.  Their  variation  is  important 
because  few  sources  were  found  that  suggest  baseline  values  to  use  for  emissivity  or  flame  transmissivity. 

The  ignition  temperature  (Tig)  has  several  effects  on  the  simulation  results  and  was  therefore 
selected  for  variation.  The  time  at  which  the  wall  lining  material  first  ignites  is  dependent  upon  its 
ignition  temperature.  After  ignition,  the  rate  of  fire  growth  is  affected  because  a  non-burning  element  is 
marked  as  burning  when  it  has  been  heated  to  its  ignition  temperature.  Additionally,  the  degree  to  which 
a  burning  element  radiates  is  proportional  to  the  fourth  power  of  its  ignition  temperature  (i.e.,  its 
maximum  surface  temperature),  thus  affecting  the  overall  level  of  radiative  exchange  in  the  comer.  The 
variation  of  ignition  temperature  was  selected  as  ±  25%  to  represent  a  reasonable  variation  due  to 
experimental  error  and  uncertainties  in  cone  or  LIFT  data  reduction.  The  minimum  temperature  at  which 
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lateral  flame  spread  occurs  (Ts  min)  was  altered  proportionally  with  ignition  temperature.  This  was  done 
because  the  decreased  ignition  temperature  would  generally  be  lower  than  the  unaltered  value  of  Tsmm. 

It  was  decided  to  vary  the  inputs  affecting  lateral  flame  spread  to  determine  the  significance  of 
lateral  flame  spread  in  the  simulations.  These  variables  were  the  minimum  temperature  at  which  lateral 
flame  spread  occurs  (Tsmi„)  and  the  lateral  flame  spread  parameter  (0).  The  variation  of  each  parameter 
was  selected  as  ±  25%  to  represent  a  reasonable  variation  due  to  experimental  error  and  uncertainties  in 
cone  or  LIFT  data  reduction.  The  variation  of  &by  ±  25%  is  straightforward;  however,  some  difficulties 
arise  when  varying  Ts  mi„  by  +  25%  due  to  its  close  proximity  with  the  ignition  temperature.  If  Ts  min  + 
25%  was  greater  than  the  Tig,  then  TStmin  was  set  to  a  value  infinitesimally  smaller  than  Tig.  This  prevented 
the  possibility  of  division  by  zero  but  effectively  prevented  lateral  flame  spread. 

The  resulting  sensitivity  matrix  is  shown  in  Table  B5-1.  The  leftmost  column  describes  the 
manner  in  which  the  input  data  were  modified.  The  next  column  refers  to  the  particular  variation,  later 
referred  to  as  “cases.”  For  example,  case  3-  refers  to  the  baseline  simulation  of  material  #3,  and  case  6g 
refers  to  the  simulation  of  material  #6  with  the  material  emissivity  set  to  0.95. 


Table  B5-1.  -  Sensitivity  analysis  matrix. 


USCG  MATERIAL  NUMBER 

VARIATION 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Baseline 

- 

k  *  0.5 

a 

0.637 

3.354 

0.391 

0.611 

0.479 

0.267 

0.543 

0.299 

0.756 

CM 

4c 

b 

2.548 

13.416 

1.564 

2.442 

1.916 

1.068 

2.170 

1.194 

3.022 

Cone  Curve  @  Lowest  Irradiance 

c 

50 

N/A 

25 

25 

50 

50 

25 

25 

25 

Cone  Curve  @  Highest  Irradiance 

e 

100 

100 

75 

75 

100 

100 

75 

75 

75 

£  =  0.75 

t 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

e  =  0.95 

g 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

f/a  m  =  0.75 

h 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

0.75 

Tfem  =  0.25 

i 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

T/g  +  25%  {Tsmjn  +  25%) 

i 

1095(986} 

1131  {1018} 

839(799} 

839 {799} 

851 {766} 

989  {890} 

908  {817} 

828  {599} 

851  {711} 

T*  -  25%  {Ts,^  -25%) 

k 

657(591} 

679  {611} 

503  {479} 

503  {479} 

511 {460} 

593(534} 

545  {490} 

497  {359} 

511  {427} 

7V™,  +  25% 

1 

876 

905 

671 

671 

681 

791 

726 

599 

681 

rs,™  -25% 

m 

591 

611 

479 

479 

460 

534 

490 

359 

427 

0  +  25% 

n 

19.49 

19.49 

15.16 

15.39 

19.49 

19.49 

19.49 

23.78 

29.03 

0  -25% 

0 

11.69 

11.69 

9.10 

9.24 

11.69 

11.69 

11.69 

14.27 

17.42 
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5.2  Variation  of  Thermal  Conductivity 


The  main  effects  of  changing  thermal  conductivity  for  all  materials  was  a  delay  or  advance  in  the 
ignition  time  and  a  change  in  the  rate  of  fire  growth,  as  indicated  by  Figure  B5-1.  A  decrease  in  thermal 
conductivity  causes  ignition  to  occur  at  an  earlier  time  and  also  induces  more  rapid  fire  growth  after 
ignition.  This  is  because  the  diffusion  of  heat  through  the  wall  is  inhibited,  leading  to  a  more  rapid 
surface  temperature  increase.  This  is  due  to  the  fact  that  the  cone  calorimeter  curve  exhibits  rapid  decay 
after  its  peak  heat  release  rate  has  been  reached.  An  increase  in  thermal  conductivity  has  the  opposite 
effect,  as  heat  readily  diffuses  through  the  wall  (case  3b  in  Figure  B5-1).  The  surface  temperature  rises 
more  slowly,  leading  to  delayed  ignition  and  slower  fire  growth. 

A  secondary  effect  of  changing  the  thermal  conductivity  also  becomes  apparent  upon 
examination  of  Figure  B5-1 .  When  ignition  is  delayed  and  fire  growth  is  slower,  a  greater  amount  of 
unbumed  material  remains  at  the  time  when  the  ignition  burner  strength  is  increased  to  300kW.  This 
promotes  a  longer  burning  time  after  the  increase  in  burner  strength. 


Haat  Relaaee  Rata  va.  Time 
Material  #3  -  FR  Polyeater 
Variation  of  Thermal  Conductivity 


Tim*  (ft«cond«) 

Figure  B5-1.  -  Variation  of  thermal  conductivity  for  material  #3. 
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5.3  Cone  Calorimeter  Curve  at  Multiple  Irradiances 

From  Figure  B5-2  it  is  evident  that  the  predicted  results  are  highly  sensitive  to  the  usage  of 
different  cone  calorimeter  curves.  Except  for  the  non-fire-retardant  polyester  based  material,  all  of  the 
low-irradiance  cone  calorimeter  curves  produced  significant  effects  similar  to  the  one  shown  in  Figure 
B5-2.  Rapid  burnout  occurs  with  lower  irradiance  levels.  Figure  B5-3  confirms  this  observation  by 
showing  that  the  pyrolysis  fronts  are  indeed  followed  shortly  by  a  burnout  front.  This  rapid  burnout  is  a 
direct  result  of  the  cone  calorimeter  curve  for  material  #4  at  an  irradiance  of  25kW/m2,  as  shown  in 
Figure  B5-4.  Very  little  area  exists  under  this  curve;  therefore,  a  small  amount  of  fuel  is  available  at  each 
node,  inducing  rapid  burnout. 

Also  important  in  this  simulation  is  the  relative  bum  time  of  the  wall  in  comparison  to  the  ceiling. 
The  material  on  the  ceiling  is  consumed  extremely  rapidly,  with  some  nodes  burning  for  less  than  one 
minute.  Similar  results  are  exhibited  with  materials  #3  and  #9,  although  the  effect  is  less  evident  with 
material  #9.  Material  #8  did  not  show  significant  variation  among  the  different  cone  calorimeter  curves. 
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Pyrolysis  and  Burnout  Fronts  vs.  Time 
Matsrtsl  #4  -  FR  Vlnylester 
Slmultlon  4c.  Cone  Irrad.  ”  25kW/m2 
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H«at  RticaM  Rat*  in  Con*  Caloflmatai 
MatarUI  #4  @J5kW/m2  Irradlanc. 


5.4  Variation  of  Emissivity  and  Flame  Transmissivity 

Variation  of  material  emissivity  and  flame  transmissivity  did  not  cause  a  significant  change  in  the 
shape  of  the  heat  release  rate  curves,  but  rather  small  deviations  from  the  baseline  case.  Figures  5-5  and 
5-6  show  the  variation  of  these  parameters  for  material  #9,  where  it  becomes  obvious  that  altering  these 
variables  had  little  impact.  The  heat  release  rate  curves  for  all  other  materials  exhibited  similar 
perturbations  from  the  baseline  case  when  material  emissivity  and  flame  transmissivity  were  varied. 


B-33 


Heat  Raleace  Rata  va.  Tima 
Material  *3  -  FR  Polyester 
Variation  of  Ignition  Temperature 


5.5  V ariation  of  Ignition  T emper  ature 

A  change  in  the  material  ignition  temperature  produced  similar  effects  as  with  the  variation  of 
thermal  conductivity.  A  lower  ignition  temperature  produced  faster  ignition,  followed  by  more  rapid  fire 
growth  and  burnout.  A  higher  ignition  temperature  produced  delayed  ignition,  slower  fire  growth,  and 
longer  bum  times.  This  is  shown  in  Figure  B5-7. 

5.6  Variation  of  Lateral  Flame  Spread  Inputs 

Variation  of  the  minimum  temperature  for  lateral  flame  spread  and  the  lateral  flame  spread 
parameter  showed  no  effects  in  the  simulations.  This  is  due  to  the  fact  that  the  upper  wall  temperature  as 
calculated  by  CFAST  never  exceeded  the  minimum  temperature  for  lateral  flame  spread,  even  for  the 
cases  where  Tsmin  was  reduced  by  25%.  For  this  reason,  no  graphs  are  presented  that  show  the  variation 
of  Ts  min  and  &  as  these  cases  are  no  different  from  the  baseline  examples. 

5.7  Flame  spread  behavior  on  fire-retardant  materials 

The  overall  prediction  of  the  behavior  of  the  fire-retardant  materials  that  reached  flashover  before 
the  increase  in  ignition  burner  output  (e.g.,  materials  #3  and  #4)  was  not  accurate.  It  is  hypothesized  that 


B-35 


the  material  properties  derived  from  the  bench-scale  cone  calorimeter  and  LIFT  testing  cannot  be  directly 
extrapolated  to  predict  full-scale  behavior  due  to  the  fire-retardants  present  in  these  materials. 

A  possible  explanation  for  this  is  that  when  materials  are  tested  in  the  cone  calorimeter,  the  fire- 
retardants  are  able  to  successfully  inhibit  the  combustion  process.  The  result  is  an  overall  reduced  rate  of 
heat  release,  and  this  behavior  is  translated  to  the  model  prediction  when  the  resulting  cone  calorimeter 
curve  is  provided  to  the  flame  spread  algorithm  as  input.  It  is  hypothesized  that  in  the  full-scale  ISO  9705 
tests,  the  fire-retardants  function  well  after  the  burner  is  first  ignited  by  hindering  combustion  of  the 
lining  material.  As  time  passes,  the  combustion  reactions  of  the  lining  material  flames  as  well  as  the 
burner  flames  consume  a  significant  portion  of  the  fire-retardants.  At  this  point,  the  lining  material  is  able 
to  support  more  rapid  combustion. 

Figure  B5-8  shows  the  experimental  results  from  materials  #3  and  #4,  indicating  potential  times 
where  the  fire-retardants  may  be  active.  Note  that  there  is  overlap  in  the  suggested  times,  as  there  is  a  fair 
degree  of  uncertainty  in  the  behavior  of  the  fire-retardants.  Material  #9  is  not  shown  on  this  graph 
because  the  fire-retardants  appear  to  have  a  more  constant  effect  until  the  time  when  the  heat  release  rate 
of  the  ignition  burner  is  increased  and  are  overcome,  as  indicated  in  Figure  B4-7. 

This  hypothesis  was  tested  by  conducting  two  additional  simulations.  A  cone  calorimeter  curve 
from  a  similar  non-fire-retardant  material  was  substituted  for  the  fire-retardant  cone  calorimeter  curve  of 
material  #3  and  #4  to  determine  if  the  effects  of  the  retardants  could  be  removed  from  the  simulation. 

The  cone  calorimeter  curve  for  material  #8  was  substituted  for  the  original  material  #3  curve  because 
material  #8  is  a  non-fire-retardant  polyester  based  composite  and  material  #3  is  a  fire-retardant  polyester 
based  composite.  The  substitute  for  the  material  #4  cone  calorimeter  curve  came  from  a  composite 
material  found  in  Jacoby’s  (Jacoby,  1998)  work  that  was  a  non-fire-retardant  vinylester  based  composite 
material  (8mm  thick,  with  a  ceramic  fiberboard  backing).  The  remaining  material  properties  were  not 
changed.  A  comparison  of  the  fire-retardant  and  non-fire-retardant  cone  calorimeter  curves  used  for 
material  #3  and  #4  are  shown  in  Figures  5-9  and  5-10.  The  resulting  simulations  for  material  #3  and  #4 
are  presented  in  Figures  5-1 1  and  5-12. 

The  results  from  the  change  in  cone  calorimeter  curve  for  material  #3  do  not  offer  much  support 
for  the  hypothesis.  Using  the  cone  curve  from  material  #8  causes  material  #3  to  reach  flashover 
conditions  very  quickly  in  a  manner  similar  to  that  seen  with  the  original  baseline  for  material  #8. 
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Heat  Release  Rate  ve.  Time 

Material  #3  (FR  polyester)  and  Material  (4  (FR  Vinyleeter) 
Experimental  Curves,  Indicating  Fire-retardant  Activity 


Time  (seconds) 

Figure  B5-8.  -  Possible  fire-retardant  activity  for  material  #3  and  #4. 


Heat  Release  Rate  in  Cone  Calorimeter 
Old  and  New  Curves  for  Material  93 


Figure  B5-9.  -  Change  of  cone  calorimeter  curve  for  material  #3. 


B-37 


Heat  Release  Rate  ve.  Time 
Material  f4  •  New  Non-FR  Cone  Calorimeter  Curve 
Comparison  to  Old  Baseline 


0  200  400  BOO  BOO  1000  1200 

Tim*  (Mcondt) 

Figure  B5-12.  -  Results  from  non-FR  cone  calorimeter  curve  for  material  #4. 

Switching  the  cone  curve  for  material  #4  shows  better  results.  The  heat  release  rate  continues  to 
grow  after  ignition,  but  the  slop  e  of  the  continued  growth  is  still  shallower  than  that  of  the  experimental 
curve.  This  simulation  was  conducting  using  the  material  properties  ( kpcp ,  Tig,  k,  etc,)  of  the  fire-retardant 
material  and  the  cone  calorimeter  curve  of  the  non-fire-retardant  material.  It  is  possible  that  the  effect  of 
the  fire-retardants  combined  with  experimental  error  or  uncertaintedrom  cone  calorimeter  or  LIFT  data 
reduction  has  produced  slightly  inaccurate  material  properties. 

For  this  reason,  the  thermal  conductivity  was  decreased  by  a  factor  of  two  and  an  additional 
simulation  was  run.  As  shown  in  Figure  B5-12,  the  slope  of  this  new  curve  near  flashover  matches  the 
experimental  curve  more  closely.  Additionally,  the  predicted  time  to  flashover  matches  that  observed 
experimentally.  The  purpose  of  doing  this  is  to  make  the  case  that  the  fire-retardants  contained  in  these 
materials  may  be  adversely  affecting  the  accuracy  of  the  simulation.  This  is  most  likely  the  result  of  a 
combination  of  mechanisms,  including  the  skewing  of  the  experimentally  derived  material  properties  and 
the  generation  of  a  cone  calorimeter  curve  that  cannot  be  directly  extrapolated  to  predict  full-scale 
behavior  due  to  the  presence  of  fire-retardants.  This  also  illustrates  the  importance  of  conducting  a 
sensitivity  analysis  on  material  properties. 
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5.8  Additional  variations 

In  addition  to  the  parameters  discussed  in  the  previous  sections,  the  compartment  size  and 
ventilation  conditions  were  altered.  The  temperature  of  the  wall  lining  material  and  ambient  temperature 
were  also  varied.  These  effects  are  outlined  in  the  following  sections. 

5.8.1  Compartment  size  and  ventilation 

The  size  of  the  compartment  was  varied  by  increasing  its  length  and  width  by  a  factor  of  two  and 
three,  effectively  multiplying  the  compartment  volume  by  four  and  nine,  respectively.  The  height  was 
held  constant  so  the  flame  spread  on  the  ceiling  is  comparable  to  the  baseline  case.  The  specimen  size 
was  also  held  constant  so  that  the  total  amount  of  fuel  available  for  combustion  does  not  change. 
Simulations  were  conducted  for  all  nine  materials.  The  results  of  these  simulations  show  that  doubling 
the  length  and  width  of  the  compartment  has  almost  no  effect  on  the  heat  release  rate.  A  moderate  effect 
was  observed  when  the  length  and  width  of  the  compartment  were  tripled.  The  results  for  material  #3  can 
be  found  in  Figure  B5-13  and  are  typical  of  most  materials.  An  exception  to  this  was  seen  with  material 
#8  which  shows  almost  no  variation  with  compartment  size  or  ventilation  conditions.  The  predicted 
upper  gas  layer  temperatures  for  material  #3  are  shown  in  Figure  B5-14. 


Comparison  of  Experimental  HRR  to  Model  Results 
Material  #3  -  FR  Polyester 


Figure  B5-13.  -  Effect  of  compartment  size  on  HRR  for  material  #3. 
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Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 
Material  #3  -  FR  Polyeater 


Figure  B5-14.  -  Effect  of  compartment  size  on  Tgll  for  material  #3. 


The  effect  that  ventilation  has  on  the  model  predictions  was  examined  by  increasing  the  size  of 
the  ventilation  opening  from  0.8m  by  2.0m  to  2.4m  by  2.4m.  This  is  analogous  to  completely  opening  the 
swinging  walls  that  form  the  door  in  the  full-scale  test  configuration.  The  simulations  showed  a  cooler 
upper  gas  layer  and  a  higher  interface  height  compared  to  the  baseline  case.  The  change  in  heat  release 
rate  for  material  #9  is  shown  in  Figure  B5-15  and  its  corresponding  gas  layer  temperatures  are  shown  in 
Figure  B5-16. 
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Comparison  of  Experimental  HRRto  Model  Results 
Material  #9  -FR  Acrylic 


Figure  B5-15.  -  Compartment  ventilation  effect  on  HRR  for  material  #9. 


Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 

Material  #9-FR  Acrylic 


5.8.2  Initial  temperature 

The  effect  of  varying  the  initial  temperature  of  the  wall  lining  material  was  also  examined.  All 
simulations  presented  thus  far  were  run  with  an  initial  temperature  of  293. 15K  for  both  the  wall 
temperature  and  the  ambient  air  temperature.  It  was  expected  that  there  would  be  little  difference 
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between  the  results  of  a  simulation  with  an  initial  temperature  of  293.15K  compared  to  270K  with  all 
other  inputs  held  constant.  However,  this  is  not  the  case  as  the  model  exhibits  a  significant  dependence 
on  initial  temperature  conditions.  When  the  initial  temperature  is  lowered  from  293.15  K  to  270K,  the 
ignition  is  delayed  by  approximately  5  seconds.  To  validate  this  behavior,  the  integral  governing  the 
temperature  of  a  semi-infinite  solid  exposed  to  a  net  flux  that  varies  with  time  given  by  Carlsaw  and 
Jaeger  (Carlsaw,  1959)  was  numerically  integrated  using  the  trapezoidal  rule  (Bradley,  1995)  and  Gauss- 
Seidel  iteration  (Lay,  1997)  for  convergence.  The  delayed  ignition  at  a  lowered  initial  temperature  as 
exhibited  by  the  model  is  consistent  with  the  results  of  this  integration. 

After  ignition,  the  initial  fire  growth  on  the  material  with  the  reduced  initial  temperature  is 
significantly  more  rapid  than  on  the  material  with  the  higher  initial  temperature.  Burnout  also  occurs 
more  rapidly  on  the  material  with  the  lowered  initial  temperature.  This  phenomenon  is  shown  in  Figures 
5-17  and  5-18  for  materials  #3  and  #9. 


Effect  of  initial  Temperature  on  Simulatiotn  Reeults  (HRR  ve.  Time) 
Material  #3  -  FR  Polyeater 

Comparison  of  Baseline  Cases  to  Experimental  Data 


Figure  B5-17.  -  Effect  of  initial  temperature  on  heat  release  rate  for  material  #8. 
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Effect  of  Initial  Temperature  on  Slmulatioln  Reeulte  (HRR  ve.  Time) 
Material  *9-FR  Acrylic 

Compariaon  of  Baaelina  Caaae  to  Experimental  Data 


Figure  B5-18.  -  Effect  of  initial  temperature  on  heat  release  rate  for  material  #9. 


At  this  time  it  is  known  that  the  lowered  initial  temperature  causes  an  increase  in  the  overall  heat 
flux  distribution  to  the  wall.  Figure  B5-19  shows  the  total  heat  flux  to  the  wall  versus  distance  from  the 
floor  at  10  seconds  after  ignition  for  material  #3.  It  can  be  seen  that  for  a  distance  greater  than 
approximately  0.5  m,  the  total  flux  to  the  wall  for  an  initial  temperature  of  270K  is  significantly  higher 
than  that  with  an  initial  temperature  of  293 . 1 5K.  This  is  also  true  for  other  times  and  for  other  materials. 
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Total  Heat  Flux  to  Wall  vt.  Distance  at  10  Seconds  After  Ignition 
Material  f3  -  FR  Polyester 


Figure  B5-19.  -  Effect  of  initial  temperature  on  heat  release  rate  for  material  #9. 


The  increased  flux  to  the  wall  causes  an  increase  in  the  heat  release  rate  because  the  mass  loss 
rate  of  the  wall  fire  is  a  function  the  incident  flux  at  each  node.  The  mechanism  that  is  responsible  for 
this  increase  in  heat  flux  to  the  wall  is  unknown  at  this  time.  It  is  possible  that  this  behavior  is  the  result 
of  a  programming  error  or  oversight  that  occurred  as  part  of  the  algorithm  implementation  process.  It  is 
also  possible  that  the  model  has  been  properly  coded  but  the  algorithm  itself  (Mitler,  1993,  1994)  contains 
some  particularity  that  has  remained  unnoticed  to  this  point  and  is  responsible  for  this  behavior. 

It  should  be  noted  that  the  stand-alone  program  SPREAD  (Mitler,  1993)  does  not  exhibit  the 
same  dependency  on  initial  temperature  that  the  enhanced  CFAST/Mitler  implementation  does.  In  fact,  a 
higher  initial  temperature  corresponds  to  slightly  more  rapid  fire  growth  in  the  stand-alone  program. 
However,  the  changes  and  additions  outlined  in  Mitler’ s  second  paper  (Mitler,  1994)  (namely  radiative 
transfer  between  flames  on  the  ceiling  and  the  vertical  walls)  were  never  implemented  into  SPREAD.  It 
is  known  that  the  enhanced  radiation  subroutine  described  in  Section  2.3  is  not  causing  the  anomalous 
behavior  because  the  same  phenomenon  (i.e.,  more  rapid  initial  fire  growth  at  a  lower  initial  temperature) 
is  seen  when  the  enhanced  radiation  subroutine  disabled.  This  same  behavior  is  also  observed  when  the 
CFAST/Mitler  algorithm  is  run  as  a  line  burner  against  the  sidewall.  The  implication  is  that  this 
anomalous  behavior  is  related  to  the  model  formulation  or  a  coding  error,  however  the  exact  problem  has 
not  been  identified  at  this  time. 
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6.0  CONCLUSIONS 


This  report  has  described  the  development  and  testing  of  a  comer  flame  spread  algorithm 
embedded  within  the  zone  model  CFAST  (Peacock,  1993,  1997).  This  section  presents  the  strengths  and 
weaknesses  of  the  model  and  makes  recommendations  for  future  work. 

6.1  Strengths  of  Model 

Certain  modules  of  Mitler’ s  (Mitler,  1993,  1994)  algorithm  in  its  current  form  are  particularly 
powerful.  The  model  is  well  suited  for  simulating  upward  spread  on  walls.  It  couples  all  of  the  major 
variables  and  converges  to  a  solution  for  the  pyrolysis  front  position  at  each  time  step.  The  entire  wall  is 
discretized  into  elemental  strips  so  that  ignition  can  occur  at  only  one  strip,  or  over  a  range  of  strips.  The 
algorithm  also  has  a  good  pyrolysis  submodel  because  the  pyrolysis  rate  is  tracked  independently  at  each 
ignited  element  based  on  a  total  incident  flux  at  that  element  and  is  grounded  in  bench-scale  empirical 
data.  The  degree  of  preheating  is  also  tracked  independently  at  each  element.  This  is  possible  because 
the  subroutine  that  calculates  the  enhanced  radiative  exchange  in  the  comer  allows  for  a  more  detailed 
flux  distribution  than  is  feasible  with  other  flame  spread  models.  The  algorithm  in  its  current  form 
considers  a  heat  flux  that  varies  with  height  and  calculates  the  convective  and  radiative  components  of 
this  flux  separately. 

An  advantage  of  the  current  form  of  the  Mitler  flame  spread  model  is  that  it  has  the  ability  to 
easily  track  species  concentrations  in  the  upper  and  lower  layers  because  it  has  been  embedded  within 
CFAST.  This  can  be  accomplished  by  specifying  the  chemical  composition  of  the  lining  material  in  the 
CFAST  object  database.  The  zone  model  then  uses  the  heat  release  rate  of  the  burning  lining  material  to 
calculate  the  species  concentrations  that  develop.  This  is  useful  if  it  ever  becomes  desirable  to  use  the 
model  to  estimate  the  toxicity  conditions  during  a  room  comer  test. 

With  respect  to  the  simulations  on  the  nine  materials  that  were  discussed  in  this  report,  the  model 
seems  to  perform  well  with  materials  that  show  significant  burning,  at  least  initially.  Most  simulations 
follow  the  slope  of  the  experimental  heat  release  rate  curve  after  ignition  until  initial  burnout  occurs,  at 
which  point  the  slope  is  significantly  decreased.  This  is  indicative  that  the  flame  spread  and  pyrolysis 
modules  are  working  well  but  there  may  be  problems  with  the  portion  of  the  algorithm  that  governs 
burnout.  It  is  hypothesized  that  this  is  a  contributing  factor  toward  the  model’s  tendency  to  under-predict 
heat  release  rate,  except  for  material  #8  where  the  behavior  is  dominated  by  the  high  heat  release  rate  of 
the  cone  calorimeter  curve. 
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6.2  Weaknesses  of  Model 


The  flamespread  model  of  Mitler  (Mitler,  1993,  1994)  in  its  current  form  does  not  consider  gas 
phase  phenomena.  It  is  based  solely  on  heat  transfer  in  the  condensed  phase.  Seven  of  the  nine  materials 
against  which  the  model  was  validated  were  classified  as  fire-retardant.  It  is  not  known  what  type  of  fire- 
retardants  was  used  in  the  manufacture  of  these  materials;  however,  it  is  assumed  that  these  fire-retardants 
function  by  hindering  combustion  in  the  gas  phase.  The  model  is  not  able  to  directly  account  for  the  gas 
phase  effects  of  these  retardants,  although  these  gas  phase  phenomena  are  indirectly  accounted  for 
through  the  cone  calorimeter  curves  provided  to  the  model  as  input.  One  weakness  of  the  model  in  its 
current  form  is  that  the  bench-scale  tests  may  not  translate  well  to  full-scale  behavior  as  a  result  of  the 
fire-retardants  acting  in  the  gas  phase  during  combustion  in  the  cone  calorimeter.  It  was  previously 
hypothesized  that  this  is  a  significant  factor  contributing  toward  the  model’s  tendency  to  under-predict 
flame  spread  on  fire-retardant  materials. 

A  weakness  with  this  algorithm,  and  flame  spread  models  in  general,  is  that  the  experimentally 
observed  shape  of  the  pyrolysis  region  cannot  be  described  by  a  simple  geometric  expression.  As  a  result 
of  this,  an  over-simplification  of  the  shape  of  the  pyrolysis  region  is  required.  The  general  shape  of  the 
pyrolysis  region  in  the  model  of  Mitler  (Mitler,  1993,  1994)  is  shown  in  Figure  B2-1,  although  it  should 
be  noted  that  this  figure  depicts  more  rapid  lateral  spread  in  the  hot  upper  gas  layer.  The  shape  of  the 
pyrolysis  region  observed  in  the  video  record  of  the  tests  conducted  by  Janssens  (Janssens,  1998)  is 
dissimilar  to  that  simulated  by  the  modified  flame  spread  algorithm.  The  shape  of  the  pyrolysis  region  is 
not  seen  to  be  discontinuous  at  the  interface,  but  rather  it  fans  outward  toward  the  ceiling.  Additionally, 
there  are  three  other  issues  that  need  to  be  considered. 

#1)  The  simulations  revealed  a  potential  problem  with  the  algorithm’s  calculation  of  burnout, 
especially  on  the  ceiling.  Burnout  occurred  extremely  quickly  in  some  cases, 
detrimentally  affecting  the  simulation  results.  This  effect  was  even  more  pronounced  on 
the  ceiling.  Therefore,  the  simulation  of  fire  spread  across  the  ceiling  is  a  weakness  of 
the  model  in  its  current  form. 

#2)  Once  multiple  pyrolysis  zones  exist  and  they  are  sufficiently  far  apart,  the  validity  of  the 
flame  height  and  heat  flux  calculations  is  questionable  because  the  algorithm  was  written 
to  simulate  only  one  pyrolysis  zone. 

#3)  The  sensitive  dependence  on  initial  temperature  conditions  is  a  current  weakness  of  the 

model.  One  would  expect  a  slight  change  in  the  flamespread  behavior  to  be  imparted  as  a 
result  of  a  change  in  the  initial  temperature  of  the  lining  material.  However,  it  was 
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recently  discovered  that  the  model  behavior  exhibits  a  significant  dependence  on  the 
initial  temperature  conditions. 

6.3  Recommendations  for  Future  Work 

It  is  felt  that  the  results  of  the  flame  spread  model  can  be  improved  significantly  with  additional 
work.  The  model  development  is  still  in  its  incipient  stages  because  a  considerable  portion  of  the  time 
spent  on  programming  to  date  has  involved  the  insertion  of  the  algorithm  into  CFAST  and  the  coding  of 
enhanced  radiation  network.  Several  areas  have  been  identified  where  additional  work  will  likely 
improve  the  results  of  the  model  considerably.  This  model,  in  a  more  refined  version,  has  the  potential  to 
predict  room  comer  test  data  much  more  accurately  than  is  indicated  by  the  results  presented  in  this 
report.  Several  tasks  have  been  identified  that  would  be  beneficial  for  future  research.  These  are  as 
follows: 

(1)  Study  the  effect  of  fire-retardants  and  gas-phase  phenomena  to  determine  a  better  way  of 
modeling  the  complex  behavior  for  which  these  effects  are  significant.  A  possible 
solution  for  this  problem  is  to  introduce  a  time-dependent  combustion  efficiency  to 
describe  the  changing  state  of  material  combustion  from  the  fire-retardant.  This  is 
feasible  because  the  model  in  its  current  form  already  utilizes  a  constant  user-specified 
combustion  efficiency. 

(2)  Compare  additional  ISO  9705  experimental  data  to  the  results  predicted  by  the  model  for 
other  common  non-fire-retardant  materials.  This  would  be  helpful  in  further  verifying 
the  capabilities  of  the  algorithm.  In  this  project  only  one  such  material  was  simulated, 
leaving  the  possibility  that  other  issues  concerning  “simple”  materials  may  still  exist. 

(3)  Re-examine  the  manner  in  which  the  mass  loss  rate  (and  effectively  the  burn-out  rate)  is 
calculated,  particularly  on  the  ceiling.  Currently,  a  weighting  scheme  carried  over  from 
the  original  algorithm  is  used  in  the  mass  loss  rate  calculations.  This  scheme  was 
developed  for  rectangular  strips  of  constant  area  and  may  be  invalid  when  used  with  the 
circular  ceiling  strips  because  their  area  increases  with  distance  from  the  comer. 

(4)  Ensure  the  validity  of  the  portion  of  the  algorithm  that  tracks  multiple  pyrolysis  fronts. 
The  algorithm  currently  calculates  several  parameters  based  on  the  height  of  a  single 
pyrolysis  zone.  The  validity  of  these  calculations  (particularly  flame  height  and  heat  flux 
to  the  walls/ceiling)  must  be  examined  once  multiple  pyrolysis  zones  are  established  at 
separate  heights. 
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(5)  Review  the  assumptions  made  in  the  use  of  the  heat  flux  map  provided  by  Dillon  (Dillon, 
1998).  The  area  of  high  flux  near  the  ceiling  was  removed  because  ignition  occurred  near 
the  ceiling,  while  experimental  observation  places  ignition  further  down  on  the  wall. 
Experiments  designed  to  determine  why  sustained  ignition  does  not  occur  at  the  area  of 
highest  flux  would  be  a  good  first  step  in  attempting  to  implement  a  better  heat  flux  map 
into  the  algorithm.  Additionally,  the  validity  of  translating  these  two-dimensional  heat  flux 
maps  to  one-dimension  requires  examination.  Using  a  scheme  that  tracks  a  two- 
dimensional  wall  temperature  and  heat  flux  distribution  is  a  way  to  incorporate  more  detail 
into  the  model. 

(6)  Verify  the  validity  of  using  a  flame  transmissivity  factor  to  attenuate  radiation  and  ensure 
that  energy  is  being  conserved.  The  possibility  exists  that  the  use  of  this  factor  does  not 
conserve  energy  within  the  radiation  network. 

(7)  Identify  the  mechanism  that  is  causing  the  anomalous  behavior  with  regards  to  the 
sensitive  dependence  on  initial  temperature  conditions.  Determine  whether  this  is 
a  fundamental  shortcoming  of  the  algorithm  itself  or  the  result  of  a  coding  error  or 
oversight  during  the  implementation  process. 
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8.0  BASELINE  HRR  RESULTS 


Companion  of  Experimental  HRR  to  Model  Results 
Material  #1  -FR  Phenolic 
Baaellne  vs.  Experimental 


Comparison  of  Experimental  HRR  to  Model  Results 
Material  #2 .  Fire  Restricting  Material 
Baseline  vs.  Experimental 
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Comparison  of  Experimental  HRR  to  Model  Results 
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Comparison  of  Experimental  HRRto  Model  Reaulta 
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Baseline  vs.  Experimental 
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Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 
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Baseline  vs.  Experimental 


Comparison  of  Thermocouple  Reading  and  Predicted  Upper  Layer  Gas  Temperature 

Mat  ;rial  #8- Polyester 
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10.0  BASELINE  SMOKE  PRODUCTION  RATE  RESULTS 


Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  #1  -  FR  Phenolic 
Baseline  vs.  Experimental 
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Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  13  -  FR  Polyester 
Baseline  vs.  Experimental 


Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  #4  -  FR  Vinylester 
Baseline  vs.  Experimental 
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Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  IS  -  PR  Epoxy 
Baseline  vs.  Experimental 
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Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
Material  #8  -  Polyester 
Baseline  vs.  Experimental 
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Comparison  of  Experimental  Smoke  Production  Rate  to  Model  Results 
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Baseline  vs.  Experimental 
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11.0  BASELINE  FLOOR  HEAT  FLUX  RESULTS 


Heat  Flux  to  the  Floor  ve.  Time 
Material  *1  -  FR  Phenolic 
Baeeline  va.  Experimental 


Heat  Flux  to  the  Floor  va.  Time 
Material  #2  •  Fire  Reatrlcting  Material 
Baaeline  va.  Experimental 
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1.0  OVERVIEW 


The  current  version  of  the  flame  spread  model  is  an  expanded  and  improved  version  of  the  earlier 
vertical  wall  flame  spread  model  [Williams  et  al.,  1997,  Beyler  et  al.,  1997],  Additional  features  include 
two  dimensional  flame  spread,  area  source  fire  exposures,  comer/ceiling  configurations,  and  hot  layer 
effects.  The  model  retains  the  ability  to  calculate  flame  spread  on  vertical  walls  that  are  not  influenced  by 
a  comer  or  ceiling.  All  new  features  incorporated  into  the  current  flame  spread  model  are  summarized  in 
this  report.  A  model  verification  section  is  also  included.  The  model  results  when  used  to  predict  the 
USCG  ISO  9705  full-scale  room  fire  tests  are  presented  and  compared  with  the  test  data.  The  ISO  9705 
room  fire  tests  were  performed  on  a  variety  of  composite  materials  for  the  U.  S.  Coast  Guard.  The  reader 
is  referred  to  Williams  et  al.  [1997]  and  Beyler  et  al.  [1997]  for  details  regarding  the  mechanics  of 
previously  incorporated  features. 

The  flame  spread  model  retains  the  element  and  node  concept  and  the  surface-heating  algorithm 
that  was  part  of  the  original  version.  A  node  has  a  specific  coordinate  relative  to  the  base  of  the  comer 
and  an  element  is  a  region  bounded  by  four  nodes,  one  at  each  comer.  This  version  of  the  flame  spread 
model  divides  the  comer-ceiling  configuration  into  uniformly  spaced  horizontal  and  vertical  nodes. 
Symmetry  requires  that  the  ceiling  be  square  and  that  it  be  discretized  in  the  same  way  as  the  horizontal 
wall  dimension.  The  flame  spread  model  calculates  the  temperature  and  burning  condition  at  each 
element.  The  documentation  for  the  original  flame  spread  model  explains  in  detail  how  the  elements  and 
nodes  are  treated  and  the  specifics  of  the  surface  heating  equations  [Williams  et  al.,  1997].  The  focus  of 
this  document  is  to  present  the  new  features  of  the  model  and  the  results  of  the  model  predictions  of  the 
USCG  ISO  9705  tests. 

This  study  was  motivated  by  the  need  to  evaluate  the  performance  of  composite  materials.  A 
means  of  predicting  the  heat  release  rate  and  the  flame  spread  potential  of  materials  with  known 
properties  was  sought.  A  computer  flame  spread  model  seemed  ideally  suited  for  this  goal,  offering  the 
ability  to  calculate  conditions  on  an  elemental  basis  where  simple  hand  calculations  fall  short.  The  comer 
version  of  the  flame-spread  model  was  developed  for  the  U.  S.  Navy  to  predict  the  performance  of 
composite  materials  in  an  open  comer  test  configuration  [Lattimer  et  al.,  1999].  When  this  code  was 
adapted  for  use  with  the  U.  S.  Coast  Guard  ISO  9705  tests,  it  was  apparent  that  the  effect  of  the  hot  layer 
had  to  be  included.  Hence,  an  expanded  version  of  the  comer  flame  spread  model  was  developed. 
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2.0  SUMMARY  OF  NEW  FLAME  SPREAD  MODEL  CALCULATIONS 
2.1  Overview  of  the  Corner  Flame  Spread  Model  Calculations 

The  comer  flame  spread  model  uses  a  two-dimensional  flame  spread  algorithm  to  calculate  the 
total  heat  release  rate  of  a  comer  fire  as  a  function  of  time.  The  model  assumes  that  the  comer  is 
symmetrical  such  that  each  wall  segment  is  the  same  size  and  the  ceiling  is  square.  The  fire  source, 
which  can  either  be  a  line  fire  or  an  area  fire,  must  always  be  located  in  the  comer.  Figure  2-1  shows  the 
comer  configuration  after  the  wall  and  ceiling  have  ignited  for  an  area  source  fire  and  Figure  2-2  shows 
the  same  comer  configuration  for  a  line  burner  source  fire. 

Three  distinct  burning  regions  were  identified  in  comer  fires  during  the  full-scale  open  comer  fire 
tests  [Lattimer  et  al.,  1999].  These  regions,  also  shown  in  Figures  2-1  and  2-2,  are  (1)  the  area  on  the  wall 
ignited  by  the  fire  source,  (2)  the  ceiling,  and  (3)  the  upper  part  of  the  wall  ignited  by  the  ceiling  jet.  Heat 
flux  correlations  were  obtained  from  the  full-scale  open  comer  fire  test  data  for  each  of  the  three  regions. 
The  heat  flux  correlations  generally  depend  on  the  type  of  source  fire,  the  target  position  relative  to  the 
heat  source,  and  the  total  heat  release  rate  of  the  source  fire  and  the  walls. 

Two  of  the  three  regions  support  lateral  flame  spread  in  the  model.  The  portion  of  the  wall 
ignited  by  the  source  fire  laterally  spreads  away  from  the  comer  and  the  portion  of  the  wall  ignited  by  the 
ceiling  jet  laterally  spreads  downward.  The  ceiling  is  not  assumed  to  support  lateral  flame  spread  in  this 
model.  The  only  mechanism  that  a  ceiling  element  can  ignite  is  by  surface  heating  in  response  to  the 
incident  heat  flux. 

An  entirely  new  feature  to  the  flame  spread  program  is  the  ability  to  calculate  the  heat  transfer  to 
the  wall  and  ceiling  from  a  hot  gas  layer.  This  is  particularly  useful  for  predicting  flame  spread  in  small 
compartments.  When  reviewing  the  ISO  9705  test  data,  it  was  concluded  that  the  hot  gas  layer  has  a 
significant  effect  on  the  flame  spread  and  total  heat  release  rate  in  the  compartment.  The  current  method 
for  predicting  the  hot  layer  temperature  used  by  the  flame  spread  model  is  based  on  an  empirical 
correlation  that  was  modified  for  a  comer  configuration. 

An  important  task  of  the  flame  spread  model  is  to  keep  track  of  the  comer  fire  burning  stage. 
There  are  three  possible  burning  stages  during  the  comer  fire  burning  process,  though  all  three  do  not 
need  to  be  experienced  in  any  given  simulation.  The  stages  are  a  means  for  addressing  the  appropriate 
heat  release  rate  to  use  in  the  heat  flux  functions.  The  stages  depend  on  which  portions  of  the  wall  and 
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ceiling  are  burning  or  are  burned  out,  and  on  the  status  of  the  source  fire.  The  first  stage  assumes  that 
there  is  a  source  fire  and  that  the  wall  heat  release  rate  can  be  added  to  the  source  heat  release  rate  when 
evaluating  the  flame  height.  Normally  the  first  stage  describes  the  initial  burning  in  the  comer.  The 
second  stage  occurs  when  either  (1)  the  source  is  shut  off  or  (2)  the  wall  and  ceiling  have  burned  out 
beyond  the  area  exposed  by  the  source  fire.  This  stage  separates  the  source  heat  release  rate  from  the  wall 
heat  release  rate  when  calculating  the  flame  height.  The  third  stage  occurs  when  the  ceiling  has  spread 
beyond  the  wall  fire  and  has  burned  out  in  the  area  where  the  wall  fire  is  located.  This  stage  separates  the 
source  and/or  wall  heat  release  rate  from  the  ceiling  heat  release  rate  and  calculates  the  flame 
height/lengths  separately.  The  next  sections  describe  the  heat  flux  correlations  assuming  that  the  comer  is 
burning  under  the  first  stage  criteria.  The  discussion  is  readily  extended  to  the  second  and  third  stage 
burning  conditions.  The  model  automatically  determines  which  stage  to  use  at  each  time  step.  In  most 
cases,  the  first  stage  will  persist  until  the  end  of  the  simulation. 

2.2  Heat  Flux  Correlations 

The  heat  flux  incident  on  a  cell  determines  when  a  cell  will  ignite  and  how  much  heat  will  be 
released  from  that  cell  once  ignition  occurs.  Empirical  correlations  are  used  by  the  model  to  determine 
the  heat  flux  incident  on  each  cell  in  the  simulation  at  each  time  step.  Due  to  the  three-dimensional  nature 
of  the  comer  fire,  several  heat  flux  correlations  are  necessary  to  describe  the  heat  fluxes  on  the  wall  and 
the  ceiling.  As  previously  described,  three  distinct  regions  were  assumed  in  a  comer  fire:  the  walls,  the 
ceiling,  and  the  wall-ceiling  boundary.  One  or  more  heat  flux  correlations  is  necessary  for  each  of  these 
regions  to  accurately  describe  the  thermal  environment.  The  type  of  fire  in  the  comer  also  has  an  impact 
on  the  heat  fluxes  in  the  comer.  Two  types  of  source  fire  scenarios  are  currently  included  in  the  model: 
the  area  source  fire  and  the  line  source  fire.  The  model  treatment  of  both  situations  is  described  in  this 
section. 

2.2. 1  Area  Source  Fire  Scenario 

2.2. 1.1  Walls 

Initially  during  the  fire  before  the  walls  ignite,  the  incident  heat  flux  to  the  walls  is  dependent 
only  the  area  source  fire.  The  incident  heat  flux  to  the  wall  from  an  area  source  fire  located  in  the  comer 
was  found  to  vary  with  the  target  elevation  above  the  source  fire  base  and  the  lateral  distance  from  the 
comer.  The  peak  heat  flux  was  determined  to  be  a  function  of  the  source  diameter  and  is  given  by  the 
following  equation  [Lattimer  et  al.,  1999]: 
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ip,*  =  i2ori-exPr-4D;; 


(2-1) 


where  qpeak  "  is  the  peak  heat  flux  in  the  comer  (kW/m2)  and  D  is  the  diameter  of  the  source  (m)  for 

square  or  rectangular  source  fires  the  length  of  the  side  facing  the  wall  should  be  used  for  the  diameter. 
The  peak  heat  flux  occurs  in  the  comer  where  the  two  walls  join.  The  peak  heat  flux  varies  with  elevation 
above  the  base  of  the  source  fire  according  to  the  following  equations  [Lattimer  et  al.,  1999]: 
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where  q  "  (y)  is  the  maximum  heat  flux  from  an  area  source  fire  to  the  comer  where  the  two  walls 

meet  (kW/m2),  y  is  the  target  elevation  above  the  base  of  the  source  fire  (m),  and  Lftip  is  the  flame  tip 
length  measured  from  the  fire  base  (m). 


The  flame  tip  length  is  a  function  of  the  dynamic  diameter  of  the  fire,  d(t),  which  is  defined  in 
this  context  as  the  width  of  burning  wall  0.9  m  above  the  base  of  the  source  fire.  The  flame  spread  at  this 
elevation  was  measured  in  the  tests  and  was  the  assumed  diameter  when  correlating  the  data.  If  this  width 
is  less  than  the  diameter  of  the  source  fire,  then  the  source  fire  diameter  is  used.  The  flame  length  for  an 
area  source  was  determined  using: 

Lftip  =  5 .9d{t)Qj/2  (2-3) 

.  * 

where  d(t)  is  the  dynamic  diameter  (m)  and  Q^\s&  dimensionless  heat  release  rate  parameter.  The 
dimensionless  heat  release  rate  parameter  is  calculated  using  the  following  equation: 

Q(t) 

O .  =  -  —-7JT  (2-4) 

1090  d(tfn 
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where  Q(t)  is  the  total  heat  release  rate  of  the  source  fire  and  the  walls.  Note  that  the  heat  release  from 
the  ceiling  is  excluded  from  the  flame  length  calculation.  This  was  done  because  all  of  the  heat  flux 
correlations  were  developed  with  fires  located  below  the  ceiling. 

The  horizontal  heat  flux  distribution  from  the  source  was  found  to  be  a  function  of  the  distance 
from  the  comer  and  the  maximum  vertical  heat  flux  distribution: 

q" (x, y)  =  q"max  (y) exp ( -  7.5  *2  )  (2-5) 

where  q"  (x,  y)  is  the  horizontal  and  vertical  heat  flux  distribution  on  the  wall  (kW/m2)  and  x  is  the 
distance  from  the  comer  (m).  This  relation  was  applied  at  all  values  of  x. 

As  material  in  the  comer  ignites,  the  heat  flux  to  the  surface  becomes  some  combination  of  the 
heat  fluxes  generated  by  the  area  source  fire  and  those  generated  by  the  wall  fire.  Flat  wall  fires  have 
been  shown  to  produce  heat  fluxes  similar  to  those  produced  by  line  fires  against  a  flat  wall  [Quintiere  et 
al.,  1986].  This  is  attributed  to  the  flames  produced  in  both  cases  being  relatively  thin  flames  that 
buoyantly  rise  up  the  surface.  For  the  comer  configuration,  the  flames  from  burning  walls  were  assumed 
to  produce  heat  flux  levels  similar  to  those  measured  with  an  “L”  shaped  line  burner  in  the  comer. 

Cells  that  ignite  between  the  comer  and  the  edge  of  the  burner  had  an  incident  flux  that  was  equal 
to  either  result  from  Equation  2-3  or  an  incident  flux  of  50  kW/m2,  whichever  is  larger.  The  50  kW/m2 
was  determined  from  an  average  of  the  heat  fluxes  measured  over  the  width  of  the  flaming  region  in  the 
“L”  shaped  burner  tests.  The  50  kW/m2  minimum  heat  flux  was  applied  because  in  some  cases  the 
correlation  may  underpredict  the  heat  flux  to  the  wall  as  it  is  a  best-fit  function.  Because  50kW/m2  was 
the  measured  average  in  this  region,  this  limit  was  placed  on  the  correlation.  Attempts  to  use  the  line 
burner  correlation  instead  of  the  prescribed  minimum  heat  flux  resulted  in  poor  predictions.  It  is  not  clear 
why  this  occurs  at  this  time. 

As  cells  begin  to  ignite  outside  of  the  burner,  the  wall  fire  is  assumed  to  have  a  horizontal  heat 
flux  distribution, 

q(x,y)  =  6{)exy[-\.0(x  /  d(t))2],  (2-6) 

where  d(t)  is  the  distance  between  the  comer  and  the  horizontal  flaming  pyrolysis  front.  This  relation  is  a 
simplified  version  of  the  horizontal  relation  developed  from  the  “L”  line  burner  tests.  When  cells  are 
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burning  outside  the  edge  of  the  burner,  the  heat  flux  incident  upon  burning  cells  outside  the  burner  is 
taken  as  the  highest  of  the  three  values:  50  kW/m2,  the  result  from  Equation  2-5  ,  or  the  result  from 
Equation  2-6.  For  cells  that  are  not  ignited  but  the  wall  is  ignited  outside  of  the  burner  edge,  the  incident 
heat  flux  is  either  the  result  from  Equation  2-5  or  the  result  from  Equation  2-6,  whichever  is  largest.  This 
approach  allows  the  model  to  decide  when  the  wall  fire  begins  to  dominate  the  preheating  process. 


2.2. 1.2  Wall-Ceiling  Boundary 

The  heat  flux  along  the  wall-ceiling  boundary  is  physically  due  to  the  hot  ceiling  jet  traveling  out 
from  the  comer.  As  a  result,  the  heat  flux  decays  with  distance  from  the  comer  as  follows: 


This  relation  is  applied  along  the  wall-ceiling  boundary  outside  of  the  burning  region  on  the  wall. 
Along  this  region,  the  heat  flux  was  assumed  to  be  constant  over  the  depth  of  the  ceiling  jet,  which  was 
taken  to  be  5=  0.1  H  [Alpert,  1975],  Cells  over  the  ceiling  jet  depth  that  have  ignited  were  assumed  to  have 
an  incident  heat  flux  of  either  the  result  in  Equation  2-7  or  30  kW/m2,  whichever  is  greater.  Refer  to  Section 
2.2. 1 . 1  for  a  discussion  of  correlation  limits  such  as  this. 


After  the  cells  within  the  ceiling  jet  ignite,  cells  below  the  ceiling  jet  begin  getting  preheated. 
The  heat  flux  distribution  below  the  ceiling  jet  is  described  through  the  following  relation: 

q'(y)=Wev{-U[{H-y)-de(t)]}  (2-8) 

where  dcj(t)  is  distance  between  the  wall-ceiling  boundary  flame  front  and  the  ceiling.  This  relation  is 
based  on  results  from  the  “L”  shaped  burner  tests. 


2. 2. 1.3  Ceiling 

The  heat  flux  distribution  on  the  ceiling  was  found  to  be  a  function  of  the  radial  distance  from  the 
comer.  This  was  found  to  be  the  case  whether  or  not  the  ceiling  was  burning.  The  following  equations 
describe  the  heat  flux  distribution  on  the  ceiling: 
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where  qc"  (r)  is  the  incident  heat  flux  on  the  ceiling  (kW/m2)  at  a  distance  r  (m)  from  the  comer,  H  is  the 

height  of  the  ceiling  above  the  floor  (m),  and  Lf  tip  is  the  height  of  the  flame  tip  (m)  calculated  using 
Equations  2-3  and  2-4.  This  is  the  same  correlation  used  in  the  wall-ceiling  boundary  region,  except  the  x 
is  replaced  with  r.  The  radial  coordinate  of  each  ceiling  element  is  calculated  from  the  Cartesian 
coordinate  of  the  element  center: 


where  xc  is  the  x-coordinate  of  the  element  center  relative  to  the  comer  (m)  and  zc  is  the  z-coordinate  of 
the  element  center  relative  to  the  comer  (m). 


Burning  ceilings  have  been  shown  to  produce  heat  fluxes  of  20  kW/m2  [Hasemi  et  al,  1995].  For 
this  reason,  ignited  cells  on  the  ceiling  have  a  heat  flux  of  either  that  determined  in  Equation  2-9  or  20 
kW/m2,  whichever  is  greater.  Refer  to  Section  2.2.1 .1  for  a  discussion  on  correlation  limits  such  as  this. 


2.2.2  Line  Source  Fire  Scenario 

The  line  source  fire  scenario  is  treated  very  similar  to  that  of  the  area  source  fire  scenario  with  the 
exception  of  the  conditions  produced  by  the  line  source  fire  itself.  As  previously  described,  line  fire 
sources  result  in  flat,  thin  flames  that  buoyantly  rise  along  a  surface.  Line  fires  against  a  flat  wall  have 
been  shown  to  result  in  lower  heat  fluxes  compared  with  those  produced  by  the  area  fire  sources  in  a 
similar  configuration  [Quintiere  et  al,  1986;  Back  et  al.,  1994],  This  effect  is  believed  to  be  due  to  line 
fires  producing  thinner  flames  than  area  fires  [Quintiere  et  al.,  1992],  Results  from  tests  performed  in  this 
study  and  elsewhere  [Kokkala,  1993;  Hasemi  et  al,  1996]  have  shown  that  “L”  shaped  line  fires  placed  in 
a  comer  produce  lower  heat  fluxes  than  area  fires. 

The  peak  heat  flux  to  the  wall  surface  is  a  function  of  the  elevation  above  the  base  of  the  source 
only  and  is  not  a  function  of  the  source  heat  release  rate  or  the  width  of  the  base.  The  following  equation 
describes  the  maximum  heat  flux  to  the  wall  in  the  comer: 
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«max=70  ^—<0.5 

LfJip 

— — —  >  0.5  (2-11) 

Lf,tip 

where  y  is  the  target  elevation  above  the  base  of  the  source  fire  (m)  and  Lftip  is  the  flame  tip  length  (m). 
The  flame  tip  was  determined  using  Equations  2-3  and  2-4,  which  are  the  same  equations  used  in  the  area 
source  fire  scenario.  The  horizontal  heat  flux  distribution  is  given  by  the  following  equation: 

q'(x,y)=q'max(y)exp[-l.0(x/d)2]  (2-12) 

where  all  terms  have  been  defined.  Once  the  wall  begins  to  ignite,  the  heat  flux  to  the  wall  is  determined 
using  the  same  procedures  used  in  the  area  source  fire  scenario,  except  the  horizontal  heat  flux 
distribution  for  the  area  source  fire  is  replaced  by  the  one  for  the  line  source  fire  given  in  Equation  2-12. 

The  heat  flux  to  the  wall  ceiling  boundary  and  the  ceiling  were  determined  in  the  noncombustible 
boundary  testing  to  be  insensitive  to  the  type  of  source  fire  in  the  comer.  As  a  result,  the  model  uses  the 
same  procedures  outlined  in  the  area  source  fire  section  to  determine  the  heat  fluxes  in  these  regions  with 
a  line  source  fire  in  the  comer. 


2.3  Opposed  Flow  Flame  Spread 

Opposed  flow  flame  spread  is  the  mechanism  by  which  flame  spreads  laterally  along  the  wall 
from  the  comer  and  down  the  wall  below  the  ceiling.  The  opposed  flow  spread  velocity  has  been  shown 
to  be  a  function  of  the  thermal  properties,  flame  spread  properties,  and  the  surface  temperature  of  the 
elements  that  are  adjacent  to  a  burning  zone  [Quintiere  and  Harkleroad,  1985].  Thermal  and  flame  spread 
properties  for  a  material  are  determined  from  the  LIFT  apparatus  using  the  procedure  outlined  in  ASTM 
El 321  [1997].  To  model  the  opposed  flow  flame  spread,  flame  front  spread  velocity  is  calculated  using 
the  following  relation  [Quintiere  and  Harkleroad,  1985]: 


v(i,j,t) 


Q 

(kpc  )l(T  ig-T  sf 


(2-13) 


where  vfij.t)  is  the  lateral  flame  spread  velocity  across  element  i,j  at  time  t  (m/s),  0  is  the  flame  heating  parameter 
measured  in  the  LIFT  apparatus  (kW2/m3 ),  kAc  is  the  effective  thermal  inertia  and  is  also  measured  in  the  LIFT 
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apparatus  (kW2-s/m4-K2),  Tig  is  the  material  ignition  temperature  (K),  and  Ts  is  the  surface  temperature  (K).  If  Ts  is 
less  than  the  minimum  temperature  for  lateral  flame  spread,  Ts  mm,  determined  in  the  LIFT  apparatus,  then  no 
opposed  flow  flame  spread  is  calculated. 

The  three  opposed  flow  flame  material  properties,  0,  kAc,  Ts  min,  collectively  define  how  quickly 
a  material  will  propagate  flame  either  laterally  or  downward.  Materials  that  have  a  flame  heating 
parameter,  0,  equal  to  0  kW2/m3  do  not  spread  flame  laterally  or  downward,  while  materials  with  a  large 
flame  heating  parameter  tend  to  rapidly  spread  flame  laterally  or  downward.  Materials  that  spread  flame 
slowly  typically  have  a  minimum  flame  spread  temperature  close  to  their  ignition  temperature.  In  other 
words,  these  materials  need  to  be  heated  to  nearly  their  ignition  temperature  before  opposed  flow  flame 
spread  is  even  possible.  The  LIFT  kAc  value  is  derived  from  the  surface  heat  loss  coefficient  at  the 
ignition  temperature  and  the  ignition  time  data  [Quintiere  and  Harkleroad,  1985].  It  is  important  to  note 
that  this  kAc  parameter  used  for  opposed  flow  flame  spread  is  different  from  the  kAc  parameter  used  in 
the  surface-heating  algorithm,  and  the  two  parameters  are  not  interchangeable.  Consequently,  the  model 
requires  input  of  both  kAc  values. 

Opposed  flow  flame  spread  is  implemented  after  the  wall  has  ignited.  In  the  wall  region,  opposed 
flow  flame  spread  calculations  are  performed  on  cells  horizontally  adjacent  to  an  ignited  cell.  In  the  wall¬ 
ceiling  interface  region,  opposed  flow  flamespread  calculations  are  performed  on  cells  vertically  adjacent 
to  an  ignited  cell.  Elements  outside  of  the  burning  region  can  either  be  ignited  from  surface  heating  or 
through  the  opposed  flow  flame  spread.  The  maximum  ignition  time  of  an  adjacent  element  is  calculated 
using  the  opposed  flow  flame  spread  algorithm  using  the  flame  spread  velocity  and  the  dimensions  of  the 
element.  This  ignition  time  is  stored  and  updated  after  each  time  step.  If  the  element  has  not  ignited  and 
the  maximum  ignition  time  has  been  reached,  then  the  element  is  ignited  and  the  surface  temperature 
raised  to  the  ignition  temperature.  If  a  particular  element  happens  to  be  adjacent  to  both  regions,  the 
shortest  element  dimension  determines  which  direction  the  flame  spread  occurs. 

2.4  Hot  Layer  Effects 

The  effect  of  the  hot  layer  on  the  heating  and  subsequent  ignition  of  the  wall  and  ceiling  materials 
was  found  to  be  very  significant  after  reviewing  the  USCG  ISO  9705  fire  test  data.  This  effect  was  most 
pronounced  when  the  hot  gas  layer  temperature  exceeded  the  minimum  flame  spread  temperature,  Ts,mi„. 

If  the  hot  gas  layer  temperature  exceeded  the  ignition  temperature  of  the  wall  and  ceiling  materials  then 
flashover  was  found  to  be  imminent.  As  a  direct  result  of  these  findings,  it  was  apparent  that  an  accurate 
prediction  of  the  hot  layer  temperature  is  a  critical  aspect  of  the  flame  spread  model. 
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The  most  readily  integrated  approach  to  calculating  the  hot  gas  layer  temperature  in  a 
compartment  is  the  method  of  McCaffrey,  Quintiere,  and  Harkleroad  (the  MQH  Method)  [McCaffrey  et 
al.,  1981].  The  MQH  method  is  an  empirical  relation  that  uses  dimensionless  variables  to  correlate  the 
average  hot  gas  layer  temperature  in  a  compartment.  The  MQH  calculation  assumes  that  the  hot  gas  layer 
temperature  is  calculated  using  the  following  equation: 


Tr  (t)  =  1.63  JT« 


Q(t) 


\2/3 


p  oo  A.O  §  Ho 


hk  At 


\-,/3 


C  p  P  oo  Ac 


yfglfo 


+  T„ 


(2-14) 


where  Tr(t)  is  the  room  temperature  (K)  at  time  t  (s),  T4  is  the  initial/ambient  room  temperature  (K),  cp  is 
the  heat  capacity  of  air  (1.0  kJ/kg),  p4  is  the  density  of  air  at  ambient  temperature  (1.2  kg/m3),  Aa  is  the 
total  area  of  the  openings  (m2),  g  is  the  acceleration  of  gravity  (9.8  m/s2),  H0  is  the  height  of  the  opening 
(m),  hk  is  the  effective  heat  transfer  coefficient  (kW/m2-K),  and  Ar  is  the  total  wall  and  ceiling  surface 
area  of  the  compartment  (m2). 


The  effective  heat  transfer  coefficient  can  be  calculated  from  the  following  equations  [McCaffrey 
et  al.,  1981]: 


hk  = 


(2-15) 


where  k  is  the  wall  boundary  thermal  conductivity  (kW/m-K),  p  is  the  wall  boundary  density  (kW/m3),  cp 
is  the  wall  heat  capacity  (kJ/kg-K).  If  there  is  more  than  one  wall  material,  as  is  the  case  when  a 
composite  is  attached  to  a  substrate,  then  a  weighted  average  should  be  used  for  the  thermal  conductivity, 
density,  and  heat  capacity.  The  weighted  average  can  be  found  using  the  following  equation: 

TPavg  =  4  *  X  5‘TP‘  (2'16) 

0  £7 


where  TP  is  the  thermal  property,  i  is  the  ith  material,  *  is  the  wall  thickness,  (m),  and  n  is  the  total 
number  of  materials. 


Karlsson  and  Magnusson  [1991]  have  modified  the  MQH  method  to  account  for  comer  effects. 
They  found  that  using  a  connection  factor  of  1 .34  resulted  in  better  temperature  predictions  for  smoke 
layers.  This  approach  is  used  in  one  flame  speed  model. 
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2.5  Heat  of  Gasification 


The  flame-spread  model  requires  an  effective  heat  of  gasification  to  determine  the  mass  loss  rate  per 
unit  area  from  an  ignited  cell.  An  average  effective  heat  of  gasification  was  developed  using  cone 
calorimeter  data  at  three  different  irradiance  levels  where  ignition  occurred. 

The  procedure  for  determining  the  average  effective  heat  of  gasification  was  similar  to  that 
proposed  in  Quintiere  [1993].  Heat  of  gasification  was  defined  through  the  following  relation: 

■  w 

A  h.  =  %  (2-19) 

m 

where  A hg  is  the  heat  of  gasification  (kJ/kg),  qnet "  is  the  heat  flux  to  the  surface  of  the  material  (kW/m2), 

and  m"  is  the  mass  loss  rate  per  unit  surface  area  (kg/m2).  The  heat  of  gasification  was  used  in  the  model 
to  determine  the  heat  released  by  the  material.  As  such,  the  mass  loss  rate  per  unit  area  was  determined 
using  the  following  equation: 


O" 

m"  =  (2-20) 

A  He 

where  Q"  is  the  test  average  heat  release  rate  per  unit  area  (kW/m2)  and  AHC  is  the  effective  heat  of 
combustion  (kJ/kg).  The  net  heat  flux  to  the  surface  of  a  material  burning  in  a  cone  calorimeter  test  was 
determined  from  a  heat  balance  at  the  material  surface: 

€e,  =  <i,mp  +  4/  -  floss  (2-21) 

where  q'mp  is  the  imposed  irradiance  from  the  cone  heater  (kW/m2),  qf  is  the  heat  flux  from  the  flame  to 

the  target  surface  (kW/m2),  and  qloss  is  the  heat  loss  from  the  target  surface  (kW/m2).  Heat  losses  are  due 

to  both  convection  and  radiation  and  were  calculated  using  the  following  equation  [Quintiere  and 
Harkleroad,  1985]: 


q,„ss "  =  he (Tig -TJ  +  ea( T%  -  Tt)  (2-22) 

where  hc  is  the  convective  heat  transfer  coefficient  (0.015  kW/m2-K),  e  is  the  surface  emissivity  (0.95),  o 
is  the  Stefan-Boltzman  constant  (5.67E-1 1  kW/m2-K4),  Tig  is  the  ignition  temperature,  and  T4  is  the 
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ambient  temperature  (298  K).  Because  the  ignition  temperature  is  known,  the  heat  losses  from  the  surface 
could  be  calculated.  Inserting  Equations  2-20  and  2-21  into  Equation  2-19  results  in  the  following 
relation: 


O' 

tip  -  floss  =  -  qf  (2-23) 

where  all  terms  have  been  defined.  By  using  the  cone  calorimeter  data  at  different  irradiance  levels,  a 
plot  of  (  q'mp  —  q'oss )  versus  Q"  /A  Hc  can  be  generated  that  should  yield  a  straight  line.  The  slope  of  this 
line  is  the  average  effective  heat  of  gasification  and  the  j-intercept  is  the  negative  of  the  flame  heat  flux. 

The  results  of  using  this  approach  to  calculate  the  average  effective  heat  of  gasification  for  the 
materials  in  this  study  are  shown  in  Table  2-2.  For  most  of  the  composites,  this  analysis  resulted  in  flame 
heat  flux  values  much  higher  than  expected  for  the  small-scale  apparatus,  where  flame  heat  flux  are 
expected  to  range  from  20-30  kW/m2.  In  addition,  the  analysis  resulted  in  a  negative  average  effective 
heat  of  gasification  and  a  negative  flame  heat  flux  for  the  FR  Epoxy.  To  determine  whether  these 
inconsistent  values  were  due  to  the  calculation  method,  a  variety  of  other  more  common  materials  were 
evaluated. 

Cone  calorimeter  and  ignition  data  for  a  variety  of  common  plastics  were  taken  from  Hirschler 
[1992]  and  Tewarson  [1995].  The  data  for  the  plywood  was  taken  from  Lattimer  et  al.,  [1999].  The  data 
and  the  calculated  average  heat  of  gasification  are  also  shown  in  Table  2-2.  The  average  effective  heat  of 
gasification  ranged  between  2.0  and  4.6  MJ/kg  and  the  flame  heat  flux  ranged  between  17  and  3 1  kW/m2. 
The  flame  heat  flux  values  for  these  materials  are  in  the  expected  range,  verifying  that  the  approach  used 
to  determine  the  heat  of  gasification  is  appropriate. 

The  reason  for  the  high  flame  heat  flux  values  calculated  for  the  composite  materials  is  not 
known.  An  overestimate  in  the  flame  heat  flux  value  results  in  an  average  effective  heat  of  gasification 
that  is  higher  than  expected,  making  the  results  less  conservative.  More  conservative  values  for  the 
average  effective  heat  of  gasification  of  the  composite  materials  were  determined  by  assuming  the  flame 
heat  flux  is  20  kW/m2,  which  is  the  average  flame  heat  flux  determined  for  the  plastics  and  wood  in  Table 
2-2.  The  results  of  the  data  analysis  with  the  flame  heat  flux  set  to  20  kW/m2  are  shown  in  Table  2-3. 

The  heat  of  gasification  was  determined  to  range  between  4.2  and  20.5  MJ/kg.  These  heat  of  gasification 
values  determined  with  a  flame  heat  flux  of  20  kW/m2  were  used  in  flame  spread  model  calculations 
presented  in  this  analysis. 
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Table  2-2.  Data  Needed  to  Determine  the  Heat  of  Gasification  Using  the  Best-Fit  Method. 


Material 

Tig( K) 

A Hc 
(MJ/kg) 

Cone  Test  Average  Q"  (kW/m2) 

qf" 

(kW/m2) 

Ahg 

(MJ/kg) 

25 

50 

_ 

75 

100 

|  U.  S.  Coast  Guard  Materials  | 

j  1-FR  Phenolic 

876 

8.22 

N.I. 

29 

49 

54 

17 

KM  | 

«aa— 

905 

9.63 

N.I. 

N.I. 

24 

35 

24 

21.9  j 

|  3-FR  Polyester 

657 

11.26 

61 

_  .. .. 

80 

100 

- 

68 

14.4  j 

1  4-FR  Vinyl  Ester 

650 

13.4 

80 

90 

120 

- 

76 

15.5 

|  5-FR  Epoxy 

726 

8.7 

N.I. 

53 

50 

42 

-256 

-37.0 

916  • 

7.67 

N.I. 

N.I. 

30 

35 

121 

38.4 

|  7-Textile  Wall 

Covering 

756 

9.08 

22 

45 

50 

-- 

35 

14.3 

I  8-Polyester 

610 

21.6 

225 

300 

390 

-- 

62 

17.5  j 

|  9- FR  Acrylic 

658 

12.33 

50 

70 

85 

- 

62 

17.5  1 

|  Other  Materials  | 

I  Plywood4 

630 

11.74 

87 

135 

21 

- 

21 

4.7  | 

|  PMMA3 

598 

24.23 

380 

660 

880 

- 

31 

2.4  | 

658 

42.6 

550 

766 

1,283 

- 

29 

2.8  j 

598 

43.4 

385 

36 

1,260 

- 

17 

2.5 

|  Nylon3 

598 

27.9 

369 

629 

1,113 

-- 

19 

1.9  j 

Phenolic-resin-impregnate  glass  fiber  core  with  malamine  facing,  overall  thickness  0.01 18  m 
2Water  based  intumescent  coating  0.5  mm  thick 

3Cone  data  from  Hirschler  [1992]  and  ignition  data  from  Tewarson  [1995] 

4Cone  and  ignition  data  from  Lattimer  et  al.,  [1999] 

N.I.  -  No  ignition 


Table  2-3.  Calculated  Heat  of  Gasification  Assuming  Flame  Heat  Flux  of  20  kW/m2. 


Material 

Tig  (K) 

qf"  (kW/m2) 

Ahg  (MJ/kg) 

1-FR  Phenolic 

876 

20 

n.o 

2-Fire  Restricting 

Material1 

905 

20 

20.5  | 

3-FR  Polyester 

657 

20 

8.0  | 

4-FR  Vinyl  Ester 

650 

20 

7.9  j 

5-FR  Epoxy 

726 

20 

12.9  1 

6-Coated  FR  Epoxy2 

916 

20 

14.5  | 

7-Textile  Wall  Covering 

756 

20 

11.0 

8-Polyester 

610 

20 

4.2 

9-FR  Modified  Acrylic 

658 

20 

10.2 

2Water  based  intumescent  coating  0.5  mm  thick 


C-16 


Figure  2.1.  Corner  Configuration  with  Area  Source  Fire. 
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Figure  2.2.  Corner  Configuration  with  Line  Burner  Source  Fire. 
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(initially) 


3.0  MODEL  VALIDATION 


The  flame  spread  model  predictions  have  been  compared  to  full-scale  fire  test  data  to  validate  the 
calculation  algorithms.  Both  ISO  9705  room  fire  test  data  and  open  comer  fire  test  data  were  used  in  the 
program  validation  process.  Three  materials  were  used  in  the  validation  process:  fire  retardant  phenolic 
composite  [Janssens  et  al.,  1998];  vinyl  ester  composite  with  thirty  and  forty  percent  resin  [Lattimer  et  al, 
1999];  and  Douglas  fir  plywood  [Lattimer  et  al.,  1999].  These  materials  were  selected  because  the 
material  properties  represented  a  range  of  flame  spread  potential.  Also,  to  show  the  effects  of  the  hot  gas 
layer,  the  ISO  9705  fire  tests  were  modeled  using  the  modified  MQH  correlation.  The  predicted  heat 
release  rate  and  the  temperature  of  the  hot  gas  layer,  when  applicable,  were  compared  to  measured  test 
data. 

The  thermal  properties  of  the  three  materials  are  shown  in  Table  3-1  and  the  flame-spread 
properties  are  shown  in  Table  3-2.  Plywood  and  the  thirty-  percent  resin  vinyl  ester  are  not  listed  in  Table 
3-1  because  they  were  not  tested  in  the  ISO  9705  room  and  were  therefore  not  influenced  by  a  hot  gas 
layer.  The  substrate  is  listed  in  Table  3-1  because  it  has  a  significant  impact  on  the  heat  loss  through  the 
compartment  boundaries.  Two  different  types  of  vinyl  ester  were  used  each  with  slightly  different  flame 
spread  properties.  The  material  with  the  greatest  flame  spread  potential  was  the  plywood,  follot^ed"by  the 
vinyl  ester,  and  then  the  fire  retardant  phenolic  resin.  The  minimum  flame  spread  temperature  for  the 
forty-percent  vinyl  ester  was  assumed  to  be  50K  less  than  the  ignition  temperature  because  the  reported 
value  was  greater  than  the  ignition  temperature.  A  flame  spread  temperature  50K  less  than  the  ignition 
temperature  is  the  minimum  observed  among  the  materials  that  had  a  flame  spread  parameter  greater  than 
zero  and  is  consistent  with  the  value  for  the  thirty-percent  vinyl  ester. 


Table  3-1.  Thermal  Material  Properties  Used  to  Calculate  Hot  Layer  Temperature. 


Material 

Thickness  (m) 

Density  (kg/m3) 

Conductivity 

(kW/m-K) 

Heat  Capacity 

(kJ/kg-K) 

F.R.  Phenolic 

0.0038 

1,750 

0.00035 

1.35 

Vinyl  Ester  (30%  Res) 

0.0048 

1,630 

0.00035 

1.26 

Calcium  Silicate 
(Substrate) 

0.02 

700 

0.0001 1 

1.1 
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Table  3-2.  Flame  Spread  Properties  of  Materials  Used  for  Model  Validation. 


Material  Property 

Material 

Vinyl  Ester  (40% 
Resin)1 

Vinyl  Ester 

(30%  Resin)1 

Fire  Retardant 
Phenolic2 

Plywood1 

Ignition  Temp  (K) 

643 

675 

876 

623 

T,ml„(K)  • 

5933 

630 

876 

448 

<J>  (kW/m3) 

6.26 

1.8 

0 

9.85 

kAc  (LIFT)  (kW2-s/m4-K2) 

1.89 

2.77 

N/A 

0.84 

k Ac  (Cone)  (kW2-s/m4-K2) 

0.624 

0.456 

0.284 

0.766 

Fraction  Burned 

0.36 

0.3 

0.3 

0.82 

Thickness 

0.0048 

0.0125 

0.0038 

0.0092 

Density  (kg/m3) 

1,630 

1,920 

1,750 

770 

AHC  (kJ/kg) 

13,430 

10,740 

8,200 

11,740 

Latent  Heat  of  Vaporization 
(MJ/kg) 

7.94 

— .  _ 

7.05 

10.9 

4.65 

^Lattimer  et  al.,[  1 999] 

2Janssens  et  al.,  [1998] 

3TS  was  assumed  to  be  50  K  less  than  the  ignition  temperature 

3.1  Validation  Test  1:  Vinyl  Ester  (40  %  Resin)  in  ISO  9705  Room 

The  first  validation  comparison  used  the  test  data  from  the  vinyl  ester  test  in  the  ISO  9705 
compartment  [Lattimer  et  al.,  1999].  The  ISO  9705  room  is  2.4-m  long,  3.6-m  wide,  and  2.4  m  high. 
There  is  a  single  opening  that  is  0.8  m  wide  and  2.0  m  tall.  The  vinyl  ester  was  attached  to  a  calcium 
silicate  board  substrate  material  that  was  0.02  m  thick.  The  source  fire  was  a  0.17-mby  17-m  area  burner. 
The  heat  release  rate  of  the  source  fire  was  100  kW  for  the  first  ten  minutes  and  then  300  kW  for  the  last 
ten  minutes,  as  prescribed  by  ISO  9705  [1990]  and  Resolution  MSC.40(64)  [1994].  The  goal  of  these 
tests  was  to  determine  if  the  composite  material  lining  the  walls  and  ceiling  could  cause  flashover  in  the 
space.  Because  this  test  was  conducted  in  a  compartment,  the  hot  gas  layer  was  an  important  aspect  of 
the  test. 


Figure  3-1  shows  the  predicted  and  actual  heat  release  rate  for  this  test  configuration.  The  actual 
heat  release  rate  and  the  predicted  heat  release  rate  suggest  that  flashover  occurred  before  400  seconds. 
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The  model  predicted  a  more  gradual  rise  in  the  heat  release  rate  until  about  350  seconds,  at  which  time  a 
rapid  growth  in  heat  release  occurred  and  flashover  occurred.  This  rapid  heat  release  rate  increase 
corresponds  to  the  time  the  hot  layer  exceeded  the  ignition  temperature  of  the  vinyl  ester  lining  the  walls 
and  ceiling.  When  the  hot  gas  layer  exceeds  the  ignition  temperature  of  the  lining  material,  ignition  of  all 
exposed  surfaces  is  imminent  and  rapid.  The  actual  test  data  showed  a  faster  heat  release  rate  growth  and 
a  flashover  time  at  about  300  seconds.  Comparison  of  the  hot  layer  temperatures  shows  good  agreement 
up  to  about  250  seconds,  about  the  time  that  flashover  occurred  in  the  ISO  9705  compartment.  Although 
the  model  predictions  for  this  fire  and  material  are  not  unacceptable,  it  is  expected  that  they  could  be 
improved  if  the  flame-spread  data  were  regenerated.  The  minimum  temperature  for  flame  spread  had  to 
be  assumed,  as  noted  in  Table  3-2,  because  a  value  that  was  greater  than  the  ignition  temperature  was 
reported.  As  will  be  seen  in  the  sensitivity  analysis,  the  results  can  deviate  considerably  as  the 
uncertainty  in  the  flamespread  properties  increases. 

3.2  Validation  Test  2:  Vinyl  Ester  (30  %  Resin)  in  Full  Scale  Open  Corner  Test 

A  different  vinyl  ester  material  was  tested  in  an  open  comer  configuration  [Lattimer  et  al.,  1999]. 
This  vinyl  ester  had  less  resin  that  the  material  tested  in  the  ISO  compartment,  but  otherwise  had  similar 
material  properties.  The  open  comer  test  does  not  result  in  a  hot  layer  build-up,  thus  it  is  possible  to 
assess  the  impact  of  the  hot  layer  on  the  fire  development  by  comparing  these  results  with  those  presented 
in  Section  3-1 .  The  open  comer  used  in  this  test  was  1 .83  m  wide  and  2.28  m  tall.  The  ceiling  was 
1.83  mby  1.83  m. 

Figure  3-3  shows  the  predicted  and  actual  heat  release  rate  for  this  test.  The  open  comer  tests  do 
not  produce  a  hot  layer  thus  this  feature  was  turned  off  in  the  model.  The  actual  and  predicted  heat 
release  rate  are  in  agreement  for  this  test. 

3.3  Validation  Test  3:  Phenolic  Resin  Composite  in  ISO  9705  Room 

This  fire  test  was  conducted  with  a  fire  retardant  phenolic  resin  composite  in  the  ISO  9705  room 
[Lattimer  et  al.,  1999].  The  dimensions  of  the  room  and  opening  as  well  as  the  size  of  the  opening  were 
described  in  Section  3.1.  This  material  was  selected  for  validation  because  it  did  not  result  in  flashover 
despite  the  fact  that  hot  layer  was  present  in  the  space.  Figure  3-4  shows  the  actual  and  predicted  heat 
release  rate  and  Figure  3-5  shows  the  actual  and  predicted  hot  gas  layer  temperature.  In  both  cases  the 
agreement  between  the  actual  and  predicted  data  is  excellent.  The  temperature  predictions  are  within 
50  K  of  the  actual  temperature  readings  for  the  entire  20-minute  period. 
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Time  (seconds) 


Figure  3-2.  Predicted  and  Actual  Smoke  Layer  Temperature  for 
Vinyl  Ester  ISO  9705  Fire  Test 


Heat  Release  Rate  (kW) 


Figure  3-3.  Predicted  and  Actual  Heat  Release  Rate  for 
Vinyl  Ester  Open  Corner  Test 
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Heat  Release  Rate  (kW) 


Figure  3-4.  Predicted  and  Actual  Heat  Release  Rate  for  Fire  Retardant 
Phenolic  Resin  Composite  ISO  9705  Fire  Test 
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Figure  3-5.  Predicted  and  Actual  Smoke  Layer  Temperature  for  Fire  Retardant 
Phenolic  Resin  Composite  ISO  9705  Fire  Test 


3.4  Validation  Test  4:  Plywood  in  Full  Scale  Open  Corner  Test 

The  final  validation  comparison  was  done  using  the  plywood  test  data  from  the  open  comer  fire 
test  [Lattimer  et  al.,  1999].  No  hot  layer  effects  were  present.  Figure  3-6  shows  the  predicted  and  actual 
heat  release  rate  for  this  test  configuration.  The  heat  release  rate  predictions  are  in  good  agreement  up 
until  about  10  minutes,  when  the  source  fire  heat  release  rate  increases.  At  this  point,  the  model  under¬ 
predicts  the  heat  release  rate  by  about  twenty-percent.  This  may  be  due  to  the  plywood  detaching  from 
the  ceiling  and  igniting  the  wall.  This  type  of  phenomena  leads  to  a  rapid  increase  in  the  heat  release  rate 
as  is  evidenced  by  Figure  3-6  and  cannot  be  predicted  in  the  flamespread  model. 

3.5  Evaluation  of  the  Effect  of  the  Hot  Layer 

The  hue  effect  of  the  hot  gas  layer  was  evaluated  by  modeling  the  forty  percent  resin  vinyl  ester 
material  tested  in  the  ISO  9705  room  without  a  hot  layer.  This  has  been  done  and  is  shown  in  Figure  3-7. 
It  can  be  seen  that  the  vinyl  ester  material  does  not  readily  bum  and  would  not  lead  to  flashover 
conditions  if  there  were  no  hot  layer.  However,  if  the  effects  of  the  hot  layer  are  included,  this  material 
can  lead  to  flashover  within  eight  minutes. 
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Heat  Release  Rate  (kW) 


Figure  3-6.  Predicted  and  Actual  Heat  Release  Rate  of  Plywood 
Open  Corner  Fire  Test 
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Figure  3-7.  Comparison  of  Predicted  Vinyl  Ester  (40  %  Resin)  Heat  Release 
Rates:  Layer  and  No  Layer  Calculations 


4.0  MODEL  RESULTS 


The  flamespread  model  was  used  to  calculate  the  heat  release  rate,  compartment  temperature, 
smoke  production  rate,  and  floor  heat  flux  in  the  nine  ISO  9705  fire  tests.  The  heat  release  rate  and 
temperature  data  were  compared  directly  to  the  measured  fire  test  data  and  are  presented  in  Sections  4.2 
and  4.3,  respectively.  The  smoke  production  rate  and  the  floor  heat  flux,  both  of  which  were  calculated 
from  the  fire  test  data,  are  presented  in  Section  4.4  and  4.5,  respectively.  Each  ISO  9705  fire  test  used  a 
different  wall  lining  material.  The  source  fire  was  a  0.17  m2  area  burner.  The  heat  release  rate  was  100 
kW  for  the  first  ten  minutes  and  was  increased  to  300  kW  for  the  next  ten  minutes,  as  specified  in  the 
IMO  Resolution  [Resolution  MSV.41(64),  1994]. 

4.1  Material  Properties 

The  flame-spread  properties  for  the  nine  test  materials  are  shown  in  Tables  4-la,  4-lb,  and  4-lc. 
There  is  some  overlap  between  these  tables  and  Table  3-1 .  The  information  is  repeated  in  this  section  for 
completeness.  Five  of  the  material  listed  the  tables  do  not  spread  fire  laterally.  The  only  means  that  these 
materials  are  predicted  to  ignite  in  the  flame  spread  model  is  when  the  incident  heat  flux  raises  the  surface 
temperature  to  the  ignition  temperature. 


Table  4-la.  Flame  Spread  Properties  of  Materials  1-3  Tested  in  ISO  9705  Room. 


Material  Property 

Material  j 

Fire  Retardant 
Phenolic 

Fire  Restricting 
Material 

Fire  Retardant  Polyester 

Test/Material  Number 

1 

2 

3 

Ignition  Temp  (K) 

876 

905 

648 

T,  m,„  (K) 

876 

905 

598 

!  A/(kW/m3) 

0 

0 

19.5 

kAc  (LIFT)  (kW2  -s/m4-K2) 

N/A1 

N/A1 

1.65 

kAc  (CONE)  (kW2  -s/m4-K2) 

0.77 

0.6 

0.9 

Fraction  Burned 

0.3 

0.37 

0.33 

Thickness  (m) 

0.0038 

0.0118 

0.0052 

Density  (kg/m3) 

1,750 

240 

1,650 

)HC  (kJ/kg) 

8,200 

9.630 

11,280 

Latent  Heat  of  Vaporization 
(MJ/kg) 

10.9 

20.5 

7.97 

Not  available';  material  did  not  ignite 
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Table  4-lb.  Flame  Spread  Properties  of  Materials  4-6  Tested  in  ISO  9705  Room. 


Material  Property 

Material 

Vinyl  Ester 
(40%  Resin) 

Fire  Retardant  Epoxy 

Fire  Retardant  Epoxy 
(Coated) 

Teat/Material  Number 

4 

5 

6 

Ignition  Temp  (K) 

643 

726 

916 

Ts,  min  (K) 

593 1 

726 

916 

<P  (kW/m3) 

6.26 

0 

0 

kpc  (LIFT)  (kW2-s/m4-K2) 

1.89 

1.73 

8.0 

kpc  (Cone)  (kW2-s/m4-K2) 

1.16 

2.0 

1.5 

Fraction  Burned 

0.36 

0.19 

0.17 

Thickness  (m) 

0.0048 

0.0039 

0.0039 

Density  (kg/m1) 

1,630 

1,910 

. 

1,910 

AHC  (kJ/kg) 

13,430 

8,700 

7,670 

Latent  Heat  of  Vaporization  (MJ/kg) 

7.94 

12.88 

14.4 

Ts,mjn  was  assumed  to  be  50  K  less  than  the  ignition  temperature 


Table  4-lc.  Flame  Spread  Properties  of  Materials  7-9  Tested  in  ISO  9705  Room. 


Material  Property 


Material 


Textile  Wall  Covering 


Polyester 


Fire  Retardant  Acrylic 


Test/Material  Number 


8 


Ignition  Temp  (K) 


756 


610 


631 


i  (K) 


756 


507 


580 


O  (kW/m  ) 


19.1 


32.8 


kpc  (LIFT)  (kW  -s/m  -K  ) 


0.27 


0.74 


1.72 


kpc  (Cone)  (kW  -s/m  -K  ) 


0.68 


0.65 


2.0 


Fraction  Burned 


0.08 


0.58 


0.34 


Thickness  (m) 


Density  (kg/m3) 


0.00641 


0.0041 


00052 


1,760 


1,390 


1,880 


AHC  (kJ/kg) 


,080 


21,600 


12,280 


Latent  Heat  of 
Vaporization  (MJ/kg) 


Includes  backing  material 


11.02 


4.23 


10.08 


C-31 


The  thermal  material  properties  for  the  nine  test  materials  are  shown  in  Table  4-2.  The  density 
and  thickness  are  also  listed  in  Tables  4- la  through  4-lc  because  they  are  also  used  in  the  flamespread 
calculation. 


Table  4-2.  Thermal  Properties  of  Materials  Tested  in  ISO  9705  Room. 


1 - 

Material  Number 

Thickness  (m) 

Density  (kg/m3) 

Conductivity1 

(kW/m-K) 

Heat  Capacity' 
(kJ/kg-K) 

' 

0.0038 

1,750 

0.00035 

1.35 

2 

0.0118 

240 

0.00043 

1.0  i 

3 

0.0052 

1,650 

0.00035 

L26  j 

4 

0.0048 

1,630 

0.00035 

i.26  | 

5 

0.0039 

1,910 

0.00045 

1.1 

6 

0.0039 

1,910 

0.00045 

1.1 

7i 

0.0064 

1,760 

0.00069 

0.84 

8 

0.0041 

1,390 

0.0001 

1.26 

9 

0.0052 

1,880 

0.00045 

1.1 

Substrate 

0.02 

700 

0.0001 1 

1.1 

1  Conductivity  and  heat  capacity  properties  obtained  from  Laramee  [1987] 

2  Includes  backing  material 


4.2  Heat  Release  Rate  Comparisons 

Figures  4-1  through  4-9  show  the  calculated  and  predicted  heat  release  rate  for  all  nine-test 
materials.  The  results  can  be  grouped  into  three  categories:  those  materials  that  caused  flashover 
conditions  in  the  room  during  the  100  kW  source  fire  exposure,  those  that  caused  flashover  during  the  300 
kW  source  fire  exposure,  and  those  that  did  not  result  in  flashover  at  all.  The  flame-spread  model 
predicted  the  heat  release  rate  most  satisfactorily  in  seven  of  the  nine  cases.  These  seven  simulations 
included  materials  from  each  of  the  three  groups  listed  above. 

4.3  Hot  Layer  Temperature  Comparisons 

The  actual  and  predicted  hot  gas  layer  temperatures  for  the  nine  test  materials  are  shown  in 
Figures  4-10  through  4-18.  The  average  ceiling  jet  layer  temperature  data  and  the  hot  layer  temperature 
measured  at  the  door  are  shown  in  the  figures.  The  average  ceiling  jet  temperature  is  not  representative  of 
the  hot  layer  temperature.  It  was  assumed  then  that  the  hot  gas  temperatures  measured  near  the  door 
would  be  more  representatives  of  the  average  hot  layer  temperature  in  the  compartment.  When  the  hot 
layer  temperature  calculation  was  implemented  in  the  flame  spread  model,  it  was  intended  that  the 
predicted  temperatures  match  the  temperatures  measured  at  the  door. 
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Total  Heat  Release  Rate  (kW) 


Figure  4-2.  Actual  and  Predicted  Heat  Release  Rate  for  Material  #  2 
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Total  Heat  Release  Rate  (kW) 


Figure  4-5.  Actual  and  Predicted  Heat  Release  Rate  for  Material  #  5 
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Figure  4-8.  Actual  and  Predicted  Heat  Release  Rate  for  Matierial  #  8 
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Figure  4-10.  Actual  and  Predicted  Hot  Layer  Temperatures  for  Material  #  1 


Temperautre  (celsius) 


Figure  4-1 1 .  Actual  and  Predicted  Hot  Layer  Temperature  for  Material  #  2 
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Figure  4-12.  Actual  and  Predicted  Hot  Layer  Temperature  for  Material  #  3 
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600 


Figure  4-13.  Actual  and  Predicted  Hot  Layer  Temperatures  for  Material  #  4 
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TemDerature  ( Celsius) 


Figure  4-14.  MuaJ  and  Predcted  Ftat  Layer  Terrperatures  for  Material  #  5 
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Temperature  (celsius) 


Figure  4-15.  Actual  and  Predicted  Hot  Layer  Temperatures  for  Material  #  6 
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Temperature  (celsius) 


Figure  4-1 6.  Actual  and  Predicted  Hot  Layer  Temperatures  for  Material  #  7 
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Figure  4-17.  Acutal  and  Predicted  Hot  Layer  Temperatures  for  Material  #  8 
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Figure  4-18.  Actual  and  Predicted  Hot  Layer  Temperatures  for  Material  #  9 
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4.4  Heat  Flux  to  the  Floor 

The  heat  flux  to  the  floor  of  the  ISO  9705  compartment  was  computed  from  the  hot  layer  data 
using  the  following  equation: 

q  =  aTl  (4-1) 

where  q  is  the  heat  flux  to  the  floor  (kW/m2)  and  TL  is  the  temperature  of  the  hot  layer  (K).  Figures 
4-19  through  4-27  show  the  predicted  floor  heat  flux  as  a  function  of  time  for  each  of  the  nine  test 
materials. 

The  comparisons  of  the  predicted  and  measured  heat  fluxes  are  consistent  with  the  comparisons 
of  the  predicted  and  measured  hot  layer  temperatures  as  would  be  expected  from  Equation  4-1 .  As  such 
the  heat  flux  prediction  errors  are  directly  traceable  to  the  temperature  prediction  errors.  This  indicates 
that  the  assumption  that  the  layer  is  optically  thick  is  not  contributing  significantly  to  the  heat  flux 
prediction  errors.  It  is  apparent  in  Figure  4-27  that  the  radiometer  was  not  functioning  properly  for 
Test  #9. 


4.5  Smoke  Production  Rate 


The  smoke  production  rate  was  calculated  using  the  following  equation: 


SPR  =  Of  mf  =  of 


(e(<)-e,(0) 

A H, 


(4-2) 


where  SPR  is  the  smoke  production  rate  (m2/s),  0f  is  the  specific  extinction  area  as  measured  in  the  cone 
calorimeter  and  is  the  average  of  all  tests  (m2/kg),  mf  is  the  total  mass  loss  rate  of  the  wall  and  ceiling 

materials  (kg/s); )  Hc  is  the  heat  of  combustion  as  measured  in  the  cone  calorimeter  and  averaged  over  all 
tests  (kJ/kg),  Q  (t)  is  the  total  calculated  heat  release  rate  (kW)  at  time  t  (s),  and  Qs(t)  is  the  heat  release 
rate  of  the  source  fire  (kW)  at  time  t. 


Figures  4-28  through  4-36  show  the  calculated  smoke  production  rate  as  a  function  of  time  for  all 
nine  test  materials.  The  predicted  smoke  production  rate  correlated  well  for  only  four  tests,  seen  in 
Figures  4-30,  4-3 1 ,  4-32,  and  4-33.  It  is  not  known  why  the  predicted  results  do  not  correlate  well  for  the 
other  five  test  materials. 
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Figure  4-19.  Actual  and  Predicted  Heat  Flux  for  Material  #1 
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Figure  4-20.  Actual  and  Predicted  Heat  Flux  for  Material  #2 
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Figure  4-21 .  Actual  and  Predicted 
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Figure  4-22.  Actual  and  Predicted  Heat  Flux  for  Material  #4 
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Figure  4-24.  Actual  and  Predicted  Heat  flux  for  Material  #6 
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Figure  4-25.  Actual  and  Predicted  Heat  Flux  for  Material  #7 
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Figure  4-26.  Actual  and  Predicted  Heat  Flux  for  Material  #8 
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Figure  4-28.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #1 
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Figure  4-29.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #2 
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Figure  4-30.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #3 
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Figure  4-33.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #6 
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Figure  4-34.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #7 
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Figure  4-36.  Actual  and  Predicted  Smoke  Production  Rate  for  Material  #9 
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5.0  SENSITIVITY  ANALYSIS 

The  sensitivity  of  the  model  flame  spread  predictions  to  the  potential  uncertainty  in  the  material 
properties  was  investigated.  Solution  convergence  was  also  established  by  changing  the  time  step  and 
element  size  and  comparing  the  results.  This  section  describes  the  uncertainty  in  the  material  properties 
and  presents  the  results  of  a  systematic  variation  of  these  parameters  on  the  predictions  of  the  model, 

This  information  is  needed  to  highlight  the  most  significant  material  properties  that  impact  the  model 
results.  In  addition  to  the  material  properties  and  solution  convergence,  the  effect  of  the  compartment 
dimensions  and  the  opening  area  were  evaluated.  This  was  done  to  show  how  the  layer  temperature  and 
possibly  the  performance  of  a  composite  material  can  change  with  room  size. 

5.1  Solution  Convergence 

The  solution  convergence  was  investigated  by  varying  the  number  of  elements  used  on  the  wall 
and  ceiling  and  by  varying  the  time  step.  Note  that  symmetry  requires  the  ceiling  have  the  same  number 
of  elements  in  the  X  and  Z  directions.  Figure  5-1  shows  the  effect  of  the  number  of  nodes  on  the  model 
predictions  for  the  fire  retardant  polyester  (Material  3)  tested  in  the  ISO  9705  room.  Several  conclusions 
can  be  drawn  from  this  figure.  First,  the  model  is  less  sensitive  to  the  number  of  vertical  (Y)  elements  as 
it  is  to  the  number  of  horizontal  (X)  elements.  This  is  not  unexpected,  because  a  coarse  horizontal  mesh 
could  impede  the  lateral  flame  spread.  Vertical  flame  spread  is  influenced  more  by  the  size  of  the  source 
fire.  If  the  source  fire  were  small,  then  the  model  could  be  more  sensitive  to  the  number  of  vertical 
elements.  It  can  also  be  concluded  from  Figure  5-1  that  the  100x100  mesh  sized  used  to  predict  the  ISO 
9705  comer  fires  was  sufficient.  A  finer  wall  mesh  (230x230)  resulted  in  a  slightly  shifted  heat  release 
rate  curve  after  six  minutes.  The  additional  computation  time  for  such  a  fine  mesh  is  clearly  not  justified. 

Figure  5-2  shows  the  effect  of  the  time  step  size  on  the  predictions  for  Material  3.  The  time  step 
was  varied  from  0.1  second  to  10  seconds.  The  time  step  used  in  this  analysis  was  1.0  second.  Time 
steps  greater  than  2  seconds  were  not  expected  to  result  in  reliable  predictions  because  the  temperature 
rise  of  heated  elements  is  too  rapid  for  the  differential  equation  solver  over  this  amount  of  time.  Figure 
5-2  shows  that  a  5  and  10-second  time  step  result  in  completely  different  predictions.  However,  a  0.1  and 
1  second  time  step  are  nearly  the  same.  It  could  be  concluded  from  this  that  a  1 -second  time  step  was 
sufficiently  small  for  this  type  of  exposure. 
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Figure  5-1 .  Effect  of  the  Number  of  X  and  Y  Nodes  on  the  Model  Predictions 

for  Material  #  3. 
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5.2  Material  Properties 


The  sensitivity  of  the  model  predictions  to  the  uncertainty  in  the  material  properties  was 
evaluated  estimating  what  the  potential  range  of  values  for  a  specific  material  property  could  be,  then 
running  the  model  with  the  maximum  and  minimum  value  in  this  range.  Two  different  types  of  materials 
were  used  to  assess  the  impact  of  a  specific  material  property.  One  material  that  was  used  caused  the 
ISO  9705  room  to  reach  flashover  when  exposed  to  the  100  kW  fire  source  (Material  3)  and  the  other 
material  caused  the  room  to  reach  flashover  only  when  exposed  to  the  300  kW  fire  source  (Material  9). 

Some  properties  were  found  to  have  little  of  no  impact  on  the  model  predictions.  These 
properties  are: 

1 .  Heat  of  combustion  of  the  composite  material; 

2.  The  density  of  the  substrate  material; 

3.  The  thermal  conductivity  of  the  substrate  material; 

4.  The  thermal  conductivity  of  the  composite  material;  and 

5 .  The  heat  capacity  of  the  substrate  material; 

The  substrate  refers  to  the  0.02-m  thick  calcium  silicone  wall  boards  that  the  nine  composite  materials 
were  attached  to  during  the  ISO  9705  tests. 

5.2.1  Substrate  Thickness 

The  effect  of  the  substrate  thickness  on  the  model  predictions  is  shown  in  Figures  5-3  and  5-4. 
The  actual  thickness  of  the  substrate  in  the  fire  tests  was  20  mm  for  both  materials.  The  thickness  was 
arbitrarily  varied  from  10  mm  to  30  mm  to  test  the  sensitivity  of  the  results  to  this  parameter.  Figure  5-3 
shows  that  the  predicted  time  to  flashover  changes  from  about  six  minutes  for  a  10-mm  thick  substrate  to 
about  9  minutes  for  a  30-mm  thick  substrate.  Conversely,  the  substrate  thickness  has  no  appreciable 
effect  on  Material  #9,  as  evidenced  by  Figure  54.  Because  the  thickness  of  Materials  #3  and  #9  are  the 
same,  the  difference  likely  is  a  result  of  a  change  in  the  effective  thermal  properties  of  the  wall. 
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5.2.2  Composite  Thickness 


The  effect  of  the  composite  thickness  on  the  model  predictions  is  shown  in  Figures  5-5  and  5-6.  The 
actual  composite  thickness  was  5.2  mm  for  both  Material  #3  and  Material  #9.  The  thickness  was 
arbitrarily  varied  from  2.6  mm  to  10.4  mm.  Figures  5-5  and  5-6  show  that  the  impact  of  the  composite 
thickness  was  more  pronounced  with  Material  #3,  as  was  the  case  with  the  substrate  thickness.  The 
composite  thickness  will  have  an  impact  on  the  heat  loss  to  the  wall,  which  will  in  turn  influence  the  hot 
layer  temperature.  The  composite  thickness  will  also  determine  the  fire  duration,  though  that  is  not  an 
issue  with  Materials  #3  and  #9. 

5.2.3  Percent  Mass  Burned 

The  effect  of  the  percent  mass  burned  on  the  model  predictions  is  shown  in  Figures  5-7  and  5-8. 
The  range  of  values  for  the  percent  mass  burned  were  obtained  from  the  maximum  and  minimum  values 
measured  in  the  cone  calorimeter  for  each  material.  The  average  percent  mass  lost  for  Material  #3  was 
33%;  the  minimum  was  10  and  the  maximum  was  50%.  The  average  percent  mass  lost  for  material  #3 
was  34%;  the  minimum  was  1 8  and  the  maximum  was  5 1%.  The  percent  mass  influences  the  burning 
duration.  Figure  5-7  shows  that  this  has  more  of  an  impact  on  the  results  for  Material  #3  than  Material 
#9.  Other  factors  contribute  to  the  burning  duration  and  may  reduce  the  overall  effect  on  the  heat  release 
rate  predictions  for  Material  #9. 

5.2.4  LEFT  Thermal  Inertia 

The  effect  of  the  LIFT  thermal  inertia  (kpc)  on  the  model  predictions  is  shown  in  Figures  5-9  and  5-10. 
The  results  obtained  by  using  the  minimum  or  maximum  reported  value  in  one  data  set  were  compared  to 
the  results  obtained  using  the  actual  value.  The  actual  value  for  Material  #3  is  1.65  kW2-s/m4-K2  and  the 
minimum  reported  value  was  1.54  kW2-s/m4-K2.  The  actual  value  for  Material  #9  is  1.72  kW2-s/m4-K2 
and  the  maximum  reported  value  was  1 .89  kW2-s/m4-K2.  Figures  5-9  and  5-10  show  that  this  parameter 
does  not  have  a  significant  influence  on  the  model  predictions.  The  results  were  included  in  this  section 
to  demonstrate  the  difference  between  the  LIFT  thermal  inertia  and  the  cone  calorimeter  thermal  inertia. 
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Figure  5-3.  Effect  of  the  Substrate  Thickness  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-6.  Effect  of  the  Composite  Thickness  on  the  Model  Predictions  for  Material  #  9 
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Figure  5-7.  Effect  of  the  Percent  Mass  Burned  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-8.  Effect  of  the  Percent  Mass  Burned  on  the  Model  Predictions  for  Material  #  9 
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Figure  5-9.  Effect  of  the  LIFT  Thermal  Inertia  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-10.  Effect  of  the  LIFT  Thermal  Inertia  on  the  Model  Predictions  for  Material  #  9 
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5.2.5 


LIFT  Flame  Heating  Parameter 


The  effect  of  the  LIFT  flame  heating  parameter  on  the  model  predictions  is  shown  in  Figures  5-11 
and  5-12.  The  maximum  flame  heating  parameter  value  was  obtained  by  using  a  best-fit  line  from  one 
test.  The  minimum  value  was  obtained  by  using  a  best-fit  line  through  one  test  also.  The  range  for  the 
flame  heating  parameter  was  extended  for  Material  #3  by  using  the  critical  heat  flux  to  obtain  a  best-fit 
line.  This  was  assumed  to  be  a  typical  error  when  determining  this  material  property  and  resulted  in  an 
extremely  low  value,  which  effectively  would  shut  off  lateral  flame  spread  in  the  model. 

Figures  5-11  and  5-12  show  that  the  results  are  relatively  insensitive  to  the  flame  heating  parameter, 
except  for  the  erroneously  low  value  used  with  material  #3.  This  means  that  the  results  are  more  sensitive 
to  the  presence  or  absence  of  lateral  flame  spread  than  the  actual  flame  spread  velocity. 


5.2.6  Minimum  Temperature  for  Flame  Spread 

The  effect  of  the  LIFT  flame  heating  parameter  on  the  model  predictions  is  shown  in  Figures  5-13 
and  5-14.  The  range  of  values  for  the  minimum  flame  spread  temperature  were  based  on  the  potential 
error  in  determining  the  location  where  the  flame  spread  stopped  on  the  test  sample.  The  actual  minimum 
flame  spread  temperature  for  Material  #3  is  580  K  and  the  maximum  value  using  this  approach  is  646  K. 
The  actual  minimum  flame  spread  temperature  for  Material  #9  is  598  K  and  the  maximum  value  is  630  K. 
Figure  5-13  and  5-14  show  that  Material  #3  is  more  sensitive  to  this  parameter  than  Material  #9. 


5.2.7  Heat  of  Gasification 


The  effect  of  the  heat  of  gasification  is  shown  in  Figures  5-15  and  5-16.  The  maximum  and 
minimum  values  were  obtained  by  using  the  highest  and  lowest  critical  heat  flux,  heat  of  combustion,  and 
heat  release  rate  values  among  all  the  tests.  Table  5-1  summarizes  the  maximum,  minimum,  and  actual 
values  used  for  this  parameter  evaluation. 


Table  5-1.  Summary  of  Values. 


Material 

Heat  of  Combustion  (kJ/kg) 

Heat  of  Gasification  (kJ/kg) 

Actual 

Maximum 

Minimum 

Actual 

Maximum 

Minimum 

3 

b mmm 

14,280 

EB33HHI 

15,280 

9,280 

9 

8,000 

iim 

3,920 

10,240 

22.900 

3,720 
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Figure  5-1 1 .  Effect  of  the  Flame  Heating  Parameter  on  the  Model 
Predictions  tor  Material  #  3. 
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Figure  5-12.  Effect  of  the  LIFT  Flame  Heating  Parameter  on  the  Model 
Predictions  for  Material  #  9 
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Figure  5-13.  Effect  of  the  Minimum  Flame  Spread  Temperature  on  the  Model 

Predictions  for  Material  #  3. 
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Figure  5-15.  Effect  of  the  Heat  of  Gasification  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-16.  Effect  of  the  Heat  of  Gasification  on  the  Model  Predictions  for  Material  #  9 


The  results  show  that  the  heat  of  gasification  has  a  significant  impact  on  the  results.  This  is  easily 
understood  considering  the  mass  loss  rate  and  the  heat  release  rate  all  are  a  function  of  this  parameter. 

5.2.8  Critical  Heat  Flux 

The  effect  of  the  critical  heat  flux  is  shown  in  Figures  5-17  and  5-18.  ASTM  E1321  requires  that 
the  accuracy  of  the  reported  value  be  within  2  kW/m2.  The  reported  critical  heat  flux  for  Material  #3  is 

14.9  kW/m2.  The  maximum  expected  value  is  therefore  16.9  kW/m2.  The  reported  critical  heat  flux  for 
Material  #9  is  15.7  kW/m2  and  the  expected  maximum  value  is  17.7  kW/m2.  Figures  5-17  and  5-18  show 
that  Material  #3  is  more  sensitive  to  the  critical  heat  flux  than  Material  #9. 

5.2.9  Ignition  Time 

The  effect  of  the  ignition  time  is  shown  in  Figures  5-19  and  5-20.  The  range  of  values  for  this 
parameter  were  obtained  using  the  Reproducibility  equation  determined  from  the  Inter-laboratory  trials  as 
listed  in  ASTM  El 321  [1997]: 


R  =  7.4  +  0.22tig  (5-1) 

where  R  is  the  range  (s)  and  tig  is  the  ignition  time  (s).  Equation  5-1  is  applied  to  the  measured  ignition 
time  at  all  heat  flux  levels  measured  in  the  cone  calorimeter.  The  ignition  time  is  actually  entered  into  the 
comer  flame  spread  model  by  varying  the  cone  thermal  conductivity,  ignition  temperature,  and  the  heat  of 
gasification.  Figures  5-19  and  5-20  show  that  Material  #3  is  more  sensitive  to  the  ignition  time  than 
Material  #9,  which  is  nearly  insensitive  to  this  parameter. 

5.2. 10  Lateral  Flame  Spread  Velocity 

The  effect  of  the  lateral  flame  spread  velocity  on  the  model  predictions  is  shown  in  Figures  5-21 
and  5-22.  The  flamespread  velocity  was  minimized  and  maximized  by  taking  the  maximum  or  minimum 
of  all  the  parameters  that  contribute  to  the  flame  spread  velocity.  These  include  the  LIFT  thermal  inertia, 
the  LIFT  flame  heating  parameter,  and  the  minimum  temperature  for  flame  spread.  The  effect  of  each  of 
these  parameters  has  already  been  investigated  and  the  ranges  have  already  been  presented.  The  flame 
spread  velocity  parameter  changes  them  all  at  once.  Figures  5-21  and  5-22  shows  that  both  Material  #3 
and  Material  #9  are  sensitive  to  the  flame  spread  velocity. 
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Figure  5-17.  Effect  of  the  Critical  Heat  Flux  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-18.  Effect  of  the  Critical  Heat  Flux  on  the  Model  Predictions  for  Material  #  9 
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Figure  5-19.  Effect  of  the  Time  to  Ignition  in  the  Cone  Calorimeter 
Model  Predictions  for  Material 
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Figure  5-22.  Effect  of  the  Lateral  Flame  Spread  on  the  Model  Predictions  for  Material  #  9 
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5.2.1 1  Cone  Calorimeter  Heat  Release  Rate 

The  effect  of  the  heat  release  rate  on  the  predictions  of  the  model  are  shown  in  Figures  5-23  and 
5-24.  The  range  of  values  for  the  heat  release  rate  was  obtained  using  the  Reproducibility  equation  from 
the  Inter-laboratory  trials  as  listed  in  ASTM  E 1 32 1 : 

R  =  25.5  +  0.1510'  (5-2) 

where  R  is  the  heat  release  rate  range  (kW/m2)  and  Q"  is  the  measured  heat  release  rate  per  unit  area  at  a 

given  exposure  flux  (kW/m2).  Figures  5-23  and  5-24  show  that  both  materials  are  very  sensitive  to  this 
parameter. 

5.3  Room  Dimensions 

The  sensitivity  of  the  model  to  the  dimensions  of  the  room  and  the  area  of  the  vent  opening  was 
performed  to  show  the  effect  of  the  hot  layer.  Changing  the  room  size  will  cause  the  hot  layer  to  increase 
or  decrease,  thereby  exposing  the  boundaries  to  different  temperatures.  Figure  5-25  shows  the  effect  of 
increasing  the  room  floor  area  on  the  heat  release  rate  predictions  and  Figure  5-26  shows  the 
corresponding  hot  layer  temperatures.  It  is  interesting  to  note  that  increasing  the  room  area  causes  the 
material  to  flashover  after  the  source  increases  to  300  kW. 

Figure  5-27  shows  the  effect  of  increasing  the  vent  width  on  the  heat  release  rate  predictions. 
Figure  5-28  shows  the  corresponding  hot  layer  temperature.  Because  increasing  the  vent  width  increases 
the  area  and  outflow,  the  hot  layer  cools,  and  the  effect  is  analogous  to  increasing  the  floor  area. 
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Figure  5-23.  Effect  of  the  Heat  Release  Rates  on  the  Model  Predictions  for  Material  #  3. 
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Figure  5-25.  Effect  of  the  Compartment  Floor  Area  on  the  Model  Predictions  for  Material  #  3 
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Figure  5-26.  Effect  of  Compartment  Floor  area  on  the  Model  Predictions  for  Material  #  3 
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Figure  5-27.  Effect  of  the  Opening  Width  on  the  Model  Predictions  for  Material  #  3 
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6.0  SUMMARY 


The  room/comer  flame  spread  model  was  developed  to  evaluate  the  performance  of  composite 
materials  under  actual  installation  configurations.  The  model  is  an  extension  of  a  previously  developed 
flame  spread  program  that  calculated  the  flame  spread  on  vertical  walls.  The  room/comer  model  retains 
the  element  and  node  concepts  as  well  as  the  surface  heating  algorithms.  New  features  include  lateral 
flame  spread,  comer  heat  flux  correlations,  hot  gas  layer  effects,  area  source  fires,  and  a  new  means  of 
calculating  the  heat  of  gasification. 

Full-scale  fire  tests  of  combustible  materials  in  room/comers  revealed  that  there  are  three  distinct 
burning  regions  in  a  comer.  These  regions  are  the  wall,  the  ceiling,  and  the  wall-ceiling  boundary. 
Heating  and  flame  spread  in  each  of  these  regions  is  calculated  in  the  flame  spread  model  using  heat  flux 
correlations  specifically  developed  for  each  region.  Both  area  source  fires  and  line  burner  source  fires 
can  be  modeled. 

The  effect  of  the  hot  gas  layer  was  found  to  have  a  significant  impact  on  the  flame  spread  by  pre¬ 
heating  the  compartment  boundary  materials.  Consequently,  the  effects  of  the  hot  layer  were  included  in 
the  comer  flame  spread  model  using  the  Method  of  McCaffery,  Quintiere,  and  Harkleroad  [McCaffrey  et 
al.,  1981]  modified  for  comer  configurations  [Karlsson  and  Magnusson,  1991]. 

The  room/comer  flame  spread  model  requires  material  property  data  that  is  obtained  from  the 
Cone  Calorimeter  data  and  from  the  LIFT  apparatus.  These  properties  include  the  material  thermal 
inertia,  the  ignition  temperature  of  the  material,  the  flame  heating  parameter,  the  heat  of  gasification,  and 
the  minimum  flame  spread  temperature.  The  thermal  inertia  obtained  from  the  Cone  Calorimeter  data  is 
calculated  using  the  heating  algorithm  in  the  flame  spread  model  and  is  a  best-fit  value  to  the  ignition 
data.  Other  input  parameters  include  the  thermal  material  properties  of  the  boundary  materials  for  use 
with  the  hot  layer  temperature  calculation,  the  dimensions  of  the  compartment,  and  the  dimensions  of  the 
compartment  opening. 

6.1  Calculation  Results 

The  calculation  results  of  the  flame  spread  model  were  compared  to  the  full-scale  test  data  for 
nine  composite  materials  tested  in  the  ISO  9705  Test.  The  model  performance  was  evaluated  by 
comparing  the  time  to  flashover  heat  release  rate,  hot  layer  temperature,  smoke  production  rate,  and  floor 
heat  flux  predictions.  The  results  are  summarized  in  Table  6-1  time  to  flashover  for  the  heat  release  rate, 
layer  temperature,  and  smoke  production  rate.  The  test  data  was  averaged  over  thirty  seconds  for  peak 


C-104 


heat  release  rate  and  sixty  seconds  for  bench  smoke  production  rate,  sliding  intervals  to  reduce  noise.  The 
model  predictions  were  not  averaged. 

Table  6-1  shows  that  the  model  predicts  the  flashover  time  well.  Flashover  was  not  predicted  for 
Materials  1,  2,  6,  and  7,  consistent  with  test  observations.  The  predicted  time  to  Shaded  regions  indicate 
tests  that  were  terminated  due  to  severe  fire  conditions.  As  such,  the  peak  heat  release  rate  and  smoke 
production  rate  comparisons  are  invalidated  by  the  test  termination.  Average  heat  and  smoke  release  rates 
are  averaged  up  to  the  time  of  predicted  flashover  to  be  consistent  with  the  tests.  *Material  No.  7  (textile 
wall  covering)  fell  off  the  wall  during  the  Room/Comer  test.  Flashover  for  Materials  3, 4,  5,  8,  and  9  was 
within  two  minutes  of  the  measured  time  to  flashover  in  all  cases. 

The  peak  heat  release  rate  predictions  are  also  shown  in  Table  6-1,  but  may  not  be  a  good 
parameter  to  assess  the  flame  spread  model  performance  for  materials  that  resulted  in  flashover.  This  is 
because  these  tests  were  terminated  before  the  intended  twenty-minute  test  period  due  to  the  severity  of 
the  fire.  Flence,  the  experimental  data  corresponds  to  the  period  just  before  the  test  was  terminated, 
whereas  the  model  results  reflect  the  peak  heat  release  under  the  assumption  that  the  test  were  run  to 
completion.  Tests  that  resulted  in  flashover  are  shaded  in  Table  6-1.  Predicted  peak  heat  release  rates  for 
materials  that  did  not  flashover  are  25-35  percent  of  the  measured  values. 
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_ Lining  and  Ceiling  Materials.  _ 
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The  average  heat  release  rates  for  the  model  shown  in  Table  6-1  were  averaged  only  up  to  the 
time  of  flashover,  such  that  this  is  a  better  indicator  of  the  model  performance.  Overall,  the 
preflashover  heat  release  rates  compare  well  with  the  test  data.  Predicted  average  heat  release  rates 
are  50-90  percent  of  the  measured  values. 

Smoke  production  is  generally  under-predicted  in  the  room/comer  flame  spread  model.  The 
model  predicted  the  smoke  production  to  be  10-25  percent  of  the  measured  value  for  three  of  the  four 
tests  that  did  not  result  in  flashover.  The  model  over-predicted  the  smoke  production  for  Material  7. 
Although  the  reason  for  this  is  not  know,  material  was  observed  to  fall  off  during  the  test,  possibly 
reducing  the  experimentally  observed  smoke  production. 

The  smoke  production  predictions  for  the  remaining  five  materials  varied  widely.  The  average 
smoke  production  rate  for  Materials  3  and  4  were  within  10  percent  of  the  experimental  value.  On  the 
other  hand,  the  smoke  production  for  Material  5  was  predicted  to  be  25  percent  the  measured  value 
and  for  Material  8  it  was  predicted  to  be  25  times  the  experimental  value. 

These  results  are  much  worse  that  the  heat  release  rate  predictions.  This  may  be  partly  due  to 
the  fact  that  the  smoke  production  calculation  uses  the  predicted  heat  release  rate  and  the  specific 
extinction  area  measure  in  the  Cone  Calorimeter.  As  such,  there  are  additional  sources  of  uncertainty 
in  these  predictions.  Nevertheless,  the  results  are  unsatisfactory  and  are  indicative  of  a  lack  of  insight 
into  smoke  generation  in  these  fires.  It  is  possible  that  for  the  materials  that  do  not  cause  flashover,  the 
smoke  generation  is  dominated  by  material  pyrolysis,  which  is  not  ignited.  This  phenomena  is  not 
included  in  the  model. 

6.2  Model  Sensitivity 

The  sensitivity  of  the  flame  spread  model  was  evaluated  by  changing  different  input 
parameters.  The  parameters  were  modified  according  to  the  expected  uncertainty  in  the  particular 
parameter.  Such  uncertainty  included  reported  error  ranges  in  the  data  and  potential  user  input  errors 
when  processing  the  Cone  Calorimeter  data  for  input  into  the  flame  spread  model.  In  addition  to  the 
material  properties,  the  convergence  was  verified  by  changing  the  number  of  elements  and  the  solution 
time  step  size. 
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It  was  demonstrated  that  the  flame  spread  model  was  not  sensitive  to  the  number  of  elements 
and  the  time  step  size  for  the  ISO  9705  room/comer  configuration  predictions.  Less  elements  or  a 
larger  time  step  could  result  in  a  non-convergent  solution. 

Two  composites  were  used  to  test  the  model  sensitivity  to  the  material  properties:  Material  3 
and  Material  9.  Material  3  represents  a  typical  composite  that  results  in  flashover  in  the  ISO  9705 
room/comer  tests  when  exposed  to  the  100  kW  fire  source.  Material  9  represents  a  composite  that 
does  not  flashover  the  ISO  9705  room/comer  test  when  exposed  to  the  100  kW  source  fire  but  does 
when  exposed  to  the  300  kW  source  fire.  The  model  was  shown  to  be  sensitive  to  the  ignition,  the 
lateral  flame  spread,  and  the  heat  release  rate  parameters.  Material  3  was  more  sensitive  than 
Material  9  in  all  cases.  The  effect  of  room  size  and  opening  area  were  also  noted  to  be  very 
significant,  primarily  by  changing  the  hot  layer  temperature. 


C-108 


7.0  REFERENCES 


Alpert,  R.L,  “Turblent  Ceiling  Jet  Induced  by  Large-Scale  Fires,”  Combustion  Science  and 
Technology,  Vol.  1 1.  Pp.  197-213,  (1975). 

Alpert,  R.  L.  (1972),  "Calculation  of  Response  Time  of  Ceiling-Mounted  Fire  Detectors,"  Fire 
Technology,  8(3),  National  Fire  Protection  Association,  Quincy,  MA,  1972. 

ASTM  E  1321-97A,  “Standard  Test  Method  for  Determining  Material  Ignition  and  Flame  Spread 
Properties,”  1997  Annual  Book  of  ASTM  Standards,  Vol  04.01:  Building  Seals  and  Sealants;  fire 
Standards;  Dimension  Stone,  American  Society  for  Testing  and  Materials,  Philadelphia,  PA,  pp.  620- 
635,(1997). 

Back,  G.,  Beyler,  C.,  DiNenno,  P.,  and  Tatem,  P.,  1994,  “Wall  Incident  Heat  Flux  Distribution 
Resulting  from  and  Adjacent  Fire,”  Fire  Safety  Science-Proceedings  from  the  Fourth  International 
Symposium,  Editor  T.  Kashiwagi,  International  Association  for  Fire  Safety  Science,  pp.  241-252. 

Beyler,  C.L.,  Hunt,  S.P.,  Iqbal,  N.,  and  Williams,  F.W.,  “A  Computer  Model  of  Upward  Flame  Spread 
on  Vertical  Surfaces,”  Fire  Safety  Science-Proceedings  of  the  Fifth  International  Symposium, 
International  Association  of  Fire  Safety  Science  (IAFSS),  Hasemi  and  Taukuba,  Editors,  Australia, 
March  307,  1997,  pp.  297-308. 

Deal,  S.  and  Beyler,  C.  (1990),  "Correlating  Preflashover  Room  Fire  Temperatures,"  Journal  of  Fire 
Protection  Engineering,  2(2),  1990. 

Hasemi,  T.,  Yokobayashi,  S.,  Wakamatsu,  T.,  and  Ptchelintsev,  A.,  “fire  Safety  of  Building 
Components  Exposed  to  a  Localized  Fire-Scope  and  Experiments  on  Ceiling/beam  System  Exposed  to 
a  Localized  Fire,”  Proceedings  from  ASIAFLAM,  pp.  351-361  (1995). 

Hasemi,  T.,  Yoshida,  M.,  Takashima,  S.,  Kikuchi,  R.,  and  Yokobayashi,  Y.,  “Flame  Length  and  Flame 
Heat  Transfer  Correlations  in  Comer-Wall  and  Comer-Wall  Ceiling  Configurations,”  Proceedings  of 
INTERFLAME' ’96,  Eds.  Franks  and  Grayson,  pp.  179-188  (1996). 

Hirschler,  M.M.  (1992),  "Plastics:  A.  Heat  Release  from  Plastic  Materials,"  Heat  Release  in  Fires, 
Elsevier,  London,  1992. 

ISO  9705  (1990),  "Fire  Tests  -  full-scale  Room  Test  for  Surface  Products,"  International  Organization 
for  Standards,  Geneva,  Switzerland,  1990. 

Janssens,  M.  L.  and  Garabedian,  A.,  and  Gray,  W.  (1998),  "Establishment  of  International  Standards 
Organization  Criteria  for  Fire  Restricting  Materials  Used  on  High  Speed  Craft,"  USCG  Report  No. 
CG-D-22-98,  1998. 

Karlsson,  B.  and  Magnusson,  S.  E.  (1991),  "Combustible  Wall  Lining  Materials:  Numerical 
Simulation  of  Room  Fire  Growth  and  the  Outline  of  a  Reliability  Based  Classification  Procedure," 
Proceedings  of  the  Third  International  Fire  Safety  Science  Symposium,  Elsevier  Applied  Science,. 
London,  1991. 


C-109 


Kokkala,  M,,  “Characteristics  of  a  Flame  in  an  Open  Comer  of  Walls,”  Proceedings  from  Interflam 
’93,  pp.  13-24,  (1993). 

Laramee,  R.  C.  (1987),  "Forms  and  Properties  of  Composite  Materials,"  Volume  1,  Engineered 
Materials  Handbook,  Composites,  American  Society  of  Materials,  Metals  Park,  Ohio,  1987. 

Lattimer,  B.Y.,  Hunt,  S.P.,  Sorathia,  U.,  Blum,  M.,  Graciek,  T.,  MacFarland,  M.,  Lee,  A.,  and  Long, 

G.  (1999),  "Development  of  a  Model  for  Predicting  Fire  Growth  in  a  Combustible  Comer," 
NSWCCD-64-TR-1999/XX,  U.S.  Navy,  Naval  Surface  Warfare  Center-Carderock  Division,  West 
Bethesda,  MD,  September  1999. 

McCaffrey,  B.,  Quintiere,  J.,  and  Harkleroad,  M.  (1981),  "Estimating  Room  Fire  Temperatures  and  the 
Likelihood  of  Flashover  using  Fire  Test  Data  Correlations,"  Fire  Technology ,  17(2),  1981. 

Mowrer,  F.  and  Williamson,  B.  (1987),  "Estimating  Room  Temperatures  from  Fires  along  Walls  and 
in  Comers,"  Fire  Technology,  23(2),  1987. 

Quintiere,  J.G.,  Harkleroad,  M.,  and  Hasemi,  T.,  “Wall  Flames  and  Implications  for  Upward  Flame 
Spread,”  Combustion  Science  and  Technology,  48  (1),  (1986). 

Quintiere,  J.G.  and  Cleary,  T.G.,  “Heat  Flux  from  Flames  to  Vertical  Surfaces,”  ASME-Heat  and  Mass 
Transfer  in  Fire  and  Combustion  Systems,  HTD-Vol.  223,  pp.  1 1 1-120,  (1992) 

Quintiere,  J.G.  (1993),  "A  Model  for  Fire  Growth  on  Materials  in  a  Room-Comer  Test,"  Fire  Safety 
Journal,  Vol.  20  1993. 

Quintiere,  J.  and  Harkleroad,  M.  (1985),  "New  Concepts  for  Measuring  Flame  Spread  Properties," 

Fire  Safety  Science  and  Engineering,  ASTM  STP  882,  American  Society  for  Testing  and  Materials, 
Philadelphia,  PA,  1985. 

Resolution  MSC.40(64)  (1994),  "Standard  for  Qualifying  Marine  Materials  for  High  Speed  Craft  as 
Fire-Restricting  Materials,"  International  Maritime  Organization,  London,  England,  1994. 

Tewarson,  A.  (1995),  "Generation  of  Heat  and  Chemical  Compounds  in  Fires,"  Section  3,  Chapter  4, 
The  SFPE  Handbook  of  Fire  Protection  Engineering,  National  Fire  Protection  Association,  Quincy, 
MA,  1995. 

Tu,  K.  M.  and  Quintiere,  J.  Q.  (1991),  "Wall  Flame  Heights  and  External  Radiation,"  Fire  Technology, 
27(3),  1991. 

Williams,  F.  W.,  Beyler,  C.  L.,  Hunt,  S.  P.,  Iqbal  (1997),  "Upward  Flame  Spread  on  Vertical 
Surfaces,"  NRL/MR/6 180-97-7908,  Navy  Technology  Center  for  Safety  and  Survivability,  Chemistry 
Division,  Washington,  DC,  1997. 


C-110 


